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Attempts  by  many  investigators  to  measure  the  electrical 
impedance  properties  of  animal  tissues  in  situ  at  frequencies  of 
physiological  interest  (0  to  200  kHz)  have  met  with  only  moderate 
success.   The  effects  of  electrode  polarization  have  been  reduced  or 
eliminated  by  using  separate  voltage  pickup  and  current  electrode 
pairs.   However,  the  practical  limitations  of  this  method  have  not 
been  clearly  defined,  nor  have  appropriate  analytical  techniques  been 
applied  to  the  analysis  and  design  of  the  electrode  and  instrumenta- 
tion systems.   These  deficiencies  and  the  inadequacy  of  existing 
instrumentation  techniques  have  hindered  progress  in  the  impedance 
art. 

In  the  dissertation  abstracted  here,  the  ensemble  consisting 
of  four  electrodes  in  an  electrolytic  medium  (or  biological  tissue) 
is  depicted  as  a  linear  four  terminal  device  and  described  by  means 
of  a  3  x  3  matrix  of  complex  impedance  parameters.   The  matrix  elements 
are  dissected  into  electrolyte  and  interface  components  which  are 
defined  in  terms  of  the  field  geometry,  and  the  impedance  characteris- 
tics of  the  electrolyte  and  the  electrode-electrolyte  interfaces. 


xiii 


Although  the  matrix  description  is  more  direct  and  often  simpler  to 
use,  an  equivalent  network  has  also  been  evolved.   The  matrix  elements 
of  two  different  electrode  configurations  are  calculated  and  used  to 
illustrate  major  analytical  results. 

The  immediate  objective  of  most  four-electrode  impedance  measure- 
ments is  to  identify  what  has  been  defined  as  the  first  element  in  the 
impedance  parameter  matrix.   This  element  is  estimated  experimentally 
by  determining  the  ratio  of  an  arbitrary  measure  of  the  nominal  voltage 
response  to  an  arbitrary  measure  of  the  nominal  exciting  current  under 
the  assumption  that  the  current  in  the  voltage-response  electrodes 
is  zero.   The  parameters  defined  in  the  impedance  matrix  permit  a 
comprehensive  analysis  to  be  made  of  the  errors  involved  in  this 
determination.   On  the  basis  of  such  an  analysis,  the  classical  method 
of  grounding  one  current  electrode  and  recording  differentially  from 
the  potential  electrodes  has  been  found  to  be  impractical  for  small 
electrode  systems  at  low  and  high  frequencies. 

To  overcome  many  of  the  limitations  of  the  classical  method,  a 
new  technique,  called  the  virtual-ground  system,  has  been  developed. 
In  this  system,  one  potential  electrode  is  maintained  close  to  ground 
potential  by  means  of  negative  feedback.   Three  classes  of  error — 
loading,  common  mode  rejection,  and  current  control — are  analyzed  for 
the  classical  and  virtual-ground  systems.   In  addition,  several  types 
of  noise  and  methods  of  noise  filtering  are  considered  in  each  system. 
The  results  of  the  error  and  noise  analyses  demonstrate  the  advantages 
of  the  virtual-ground  system  and  provide  a  set  of  design  formulae.   To 
complement  the  virtual-ground  system,  a  current  control  circuit  has 
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been  adopted  from  the  literature  and  analyzed  for  application  to  broad 
band  impedance  measurements . 

A  successfully  operating  system  using  the  virtual-ground  system 
and  the  adopted  current  control  circuit  has  been  built  at  modest  cost 
and  is  described  in  the  dissertation.   It  has  been  tested  with  small 
electrode  assemblies  in  saline  and  with  similar  assemblies  implanted 
in  the  cerebral  cortex  of  live  rats.   Several  experiments  are  reported 
to  verify  system  performance  and  indicate  the  quality  of  results 
already  obtained. 
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CHAPTER  1 


INTRODUCTION 


Problems  and  Contributions 

Attempts  by  many  investigators  to  measure  the  electrical 
impedance  properties  of  animal  tissues  in  situ  at  frequencies  of 
physiological  interest  (0  to  200  kHz)  have  met  with  only  moderate 
success.   The  effects  of  electrode  polarization  have  been  reduced  or 
eliminated  by  using  separate  voltage  pickup  and  current  electrode 
pairs.   However,  the  practical  limitations  of  this  method  have  not 
been  clearly  defined,  nor  have  appropriate  analytical  techniques  been 
applied  to  the  analysis  and  design  of  the  electrode  and  instrumenta- 
tion systems.   These  deficiencies  and  the  inadequacy  of  existing 
instrumentation  techniques  have  hindered  progress  in  the  impedance 
art. 

In  the  present  chapter,  definitions  and  assumptions  are 
explained  which  are  basic  to  the  study  of  biological  impedance 
phenomena;  the  nature  of  tissue  impedance  is  reviewed  according  to 
current  knowledge;  and  a  review  of  selected  literature  is  presented, 
which  reveals  various  measurement  techniques  and  common  sources  of 
error.   The  material  in  Chapters  2,3,4  and  5  represents  an  original 
application  of  engineering  principles  toward  the  general  fulfillment 
of  the  needs  expressed  in  the  opening  paragraph.   In  Chapter  2, 


special  parameters  are  defined  which  lead  to  an  unprecedented  under- 
standing of  the  electrode-tissue  system  and  provide  an  appropriate 
basis  for  the  analysis  and  design  discussions  that  follow.   In  Chapter 
3,  the  classical  method  for  transfer  impedance  measurements  is 
thoroughly  investigated,  resulting  in  useful  design  criteria  and 
sophistication  in  understanding  the  major  types  of  measurement  error. 
In  Chapter  4,  a  new  instrumentation  technique  is  proposed  and  analyzed 
in  depth.   The  new  technique  is  shown  to  offer  major  advantages  over 
existing  methods.   An  experimental  system  which  incorporates  this 
technique  and  several  other  special  advances  is  described  in  Chapter  5 
along  with  pertinent  experimental  results. 

Network  Classification 

To  make  impedance  measurements  more  intelligible  at  the  present 
state  of  the  art,  one  usually  invokes  the  restrictions  of  time 
invariance  and  linearity  on  the  medium  of  interest. 

Time  invariance  implies  that  the  impedance  properties  of  the 
medium  do  not  vary  during  the  critical  period  of  measurement.   The 
critical  period  may  just  include  one  to  several  repetitions  of  the 
applied  signal  waveform  or,  for  example,  it  may  include  a  whole 
sequence  of  separate  measurements  at  different  sinusoidal  frequencies. 
Time  invariance  does  not  exclude  the  measurement  of  impedance  changes 
provided  that  the  changes  are  slow  compared  to  the  time  required  for 
a  single  measurement.   Certain  impedance  characteristics  of  tissues 
in  situ  are  known  to  vary  spontaneously  and  under  various  stimuli  in 
periods  ranging  from  many  hours  to  several  seconds. 


A  system  is  linear  if  the  response  to  a  forcing  function 
comprised  of  several  components  is  equal  to  the  sum  of  the  responses 
which  would  be  obtained  from  each  component  acting  alone  (Dern  and 
Walsh,  1963).   Further,  if  the  amplitude  of  a  given  forcing  function 
is  doubled,  then  the  amplitude  of  the  response  function  should  also 
double  while  the  waveform  remains  the  same.   This  restriction  appears 
to  be  satisfied  in  practice  if  the  charge  and  current  density  in  the 
tissue  are  kept  sufficiently  low.   Hence,  currents  and  voltages  are 
typically  measured  in  microamperes  and  microvolts.   At  these  levels, 
several  different  sources  of  noise  may  be  significant.   This  topic  is 
considered  later  in  the  work. 

When  the  linearity  condition  is  satisfied,  it  is  almost  certain 
that  the  medium  will  also  be  bilateral  and  passive.   A  system  is 
bilateral  if  the  negative  of  the  response  to  a  given  forcing  function 
is  the  response  to  the  negative  of  that  forcing  function.   Most  tissues 
are  not  strictly  passive  in  the  network  sense  (Raisbeck,  1954)  because 
they  contain  electrochemical  sources  of  energy.   In  the  linear  range, 
however,  the  electrical  activity  of  these  sources  appears  to  be  indepen- 
dent of  the  applied  signal.   Since  measurement  techniques  attempt  to 
weed  out  the  effects  of  background  noise,  the  final  results  ideally 
correspond  to  those  of  a  passive  medium. 

Complex  Impedance 

If  a  sinusoidal  forcing  function  is  applied  to  a  time  invariant 
linear  bilateral  network,  every  waveform  in  the  network  will  also  be 
sinusoidal  with  the  same  frequency  once  the  initial  transients  die 


out,  i.e.,  when  the  steady  state  is  reached.   Only  the  phase  and 
amplitude  may  differ  from  place  to  place.   This  convenient  property 
leads  to  the  definition  of  complex  impedance, 
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where  | v|  and  | l|  are  the  peak  amplitudes  of  an  arbitrary  steady  state 
voltage-current  pair  and  9  is  the  algebraic  phase  difference  between 
them,  measured  with  respect  to  the  voltage  waveform. 

The  factor  e^e  is  known  as  the  polar  form  of  the  complex  number 
given  by 

eJ°  =  cos  6  +  j  sin  6  (1.2) 

where  e  is  the  natural  logarithm  base  and  j  is  the  imaginary  operator 
/^-l.   The  cosine  and  sine  terms  are  known  as  the  real  and  imaginary 
parts  of  the  complex  number,  respectively. 

If  V  and  I  occur  at  the  same  terminal  pair,  Z  is  called  a  driving 

o 

point  impedance.   When  they  occur  at  different  terminal  pairs,  Z  is 
known  as  a  transfer  impedance.   In  either  case,  it  is  an  inherent 
physical  parameter  of  the  network.   Sometimes,  it  is  more  convenient 
to  use  the  reciprocal  of  the  complex  impedance  called  the  complex 
admittance.   For  the  same  quantities  as  before,  the  complex  admittance 
is  defined  as 


'•■m--1*  <i-3) 


The  concepts  of  complex  impedance  and  complex  admittance  will  be  used 
frequently  in  the  present  work.   In  addition,  all  voltage  and  current 
signals  in  the  electrode-tissue  system  will  be  written  in  complex  form. 
The  complex  form  of  the  signal 


v(t)  =  |v|  Cos  (tut  +  9)  (1.4) 


is  defined  as 


V  E  Ivle^1 


(1.5) 


Explanations  of  the  use  of  complex  notation  may  be  found  in  basic  texts 
in  electrical  engineering  (e.g.,  Van  Valkenburg,  1955). 

The  theory  of  linear  time  invariant  systems  shows  that  if  the 
sinusoidal  steady  state  response  of  a  system  is  known  at  a  sufficient 
number  of  frequencies,  the  response  to  other  waveforms,  periodic  or 
transient,  can  be  determined  analytically  and  vice  versa.   This 
equivalency  is  especially  welcomed  when  spectral  information  is  desired 
and  the  time  for  measurements  is  limited.   Although  useful  results  have 
been  obtained  from  nonsinusoidal  methods  (Smith,  1967),  much  of  the 
available  information  is  lost  due  to  inadequate  resolution  and  band- 
width in  the  measuring  system.   Improved  data  processing  techniques 
using  digital  computation  are  certain  to  overcome  some  of  the  difficulty. 
The  remaining  problem  lies  in  the  extraction  and  amplification  of  the 
raw  analog  data  which  is  a  major  theme  in  the  present  work. 


Tissue  Impedance 

Animal  tissues  fall  into  four  major  classifications:   epithelial, 
connective,  muscular  and  nervous  (Ham,  1965);  and  each  of  these  may 
be  further  divided  into  a  number  of  specific  subclasses.   The  constitu- 
tion and  structure  of  tissues,  even  within  the  same  classification,  may 
vary  markedly.   Furthermore,  the  electrical  properties  of  relatively 
few  tissue  types,  such  as  muscular  and  nervous,  have  been  studied  in 
detail.  The  following  discussion  presents  a  brief  synopsis  of  the 
highlights  of  tissue  impedance  studies.   It  pertains  most  directly  to 
muscular  and  nervous  tissues  but  is  also  applicable,  with  some  modifi- 
cations,  to  other  tissues. 

In  tissue,  the  applied  current  may  travel  between  cells  in  the 
extracellular  fluid  and  across  cell  membranes  into  the  intracellular 
fluid  and  subcellular  structures.   At  frequencies  of  the  order  of  200 
kHz  and  below,  the  electrical  properties  of  the  fluids  are,  for  all 
practical  purposes,  independent  of  frequency  and  are  characterized  by 
resistivity  values  of  100  to  200  ohm-cm  (Schwan,  1963).   However, 
membranes  (cellular  and  subcellular)  which  are  less  than  200  Angstroms 
thick  are  associated  with  relatively  high  transverse  resistance  (1  to 
ip00  ohm-cm  )  and  high  transverse  capacitance  (^  1  uf/cm2). 

The  charging  and  discharging  of  adjacent  membrane- fluid  interfaces 
(similar  to  plates  of  a  fixed  capacitor)  is  one  of  several  mechanisms 
which  impart  frequency  dependent  characteristics  to  the  overall  tissue 
impedance.   The  characteristics  due  to  this  specific  mechanism  in 
tissue  have  been  termed  $  dispersion  by  Schwan  (Schwan,  1957,  1959, 
1963;  Schwan  and  Cole,  1960).   The  membraneous  nature  of  tissue  is  also 


responsible,  in  less  evident  ways,  for  other  dispersion  phenomena  at 
lower  frequencies  than  the  6  dispersion.   These  phenomena  have  been 
termed  a  dispersion  (in  tissues)  by  Schwan  and  have  been  associated 
with  several  different  mechanisms  which  are  still  under  debate.   Further 
discussions  of  a   and  0  dispersion  are  given  in  later  paragraphs. 

In  addition  to  microscopic  heterogeneities  measured  in  Angstroms 
and  microns,  tissue  properties  may  vary  grossly  in  distances  ranging 
from  one  to  many  millimeters.   Many  tissues  are  also  anisotropic  (Rush 
et  al,  1963;  Nicholson,  1965;  Ranck  and  Be  Merit,  1965).   These 
characteristics  and  the  close  boundaries  of  some  tissues  in  situ  may 
be  expected  to  influence  the  impedance  determined  with  different 
electrode  configurations  in  complex  ways.   Hence,  the  concept  of 
specific  impedance,  defined  for  a  one-centimeter  cube  of  tissue  in  a 
uniform  field,  is  an  idealization  which  requires  special  qualification 
in  each  case.   To  be  sure,  data  obtained  with  different  electrode 
configurations  in  the  same  tissue  may  not  simply  differ  by  a  constant 
geometrical  factor  and  a  comparison  of  results  is  at  best  approximate. 
In  practice,  one  tries  to  use  an  electrode  system  which  measures  average 
properties  of  a  large  population  of  cells,  but  avoids  the  effects  of 
gross  inhomogeneities.   There  are,  however,  applications  in  which 
careful  localization  is  not  an  important  factor. 

For  an  arbitrary  electrode  system,  let  the  Measured  complex 
impedance  be  written 

Z  =  R(f)  +  jX(f)  (1>6) 


where  R  and  X  designate  resistance  and  reactance,  respectively,  and 
both  are  functions  of  frequency.   It  has  been  reasonably  well  established 
that  two  fairly  distinct  frequency  dispersions  occur  in  many,  if  not 
all,  tissues  in  the  0  to  200  kHz  frequency  band  (Schwan,  1957;  Fatt , 
1963;  Ranck,  1963a).   Using  Schwan 's  terminology,  the  a  and  3  disper- 
sions are  characterized  approximately  by  circular  arcs  in  the  complex 
impedance  plane  as  indicated  in  Figure  1.1.   The  arcs  usually  overlap 
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Figure  1.1.   Frequency  Dispersion  in  Tissues 

and  have  centers  below  the  R  axis.   In  results  so  far  reported,  the 
overlap  region  generally  occurs  at  frequencies  below  2  kHz.   The  low 
frequency  for  maximum  reactance,  f  ,  may  vary  from  several  hundred 
Hz  to  below  20  Hz  depending  on  tissue  and  experimental  conditions. 


Likewise,  the  high  frequency  for  maximum  reactance,  f  ,  may  vary  from 
roughly  1  kHz  to  well  above  200  kHz. 

Scales  have  been  provided  in  Figure  1.1  to  illustrate  relative 
features  of  the  dispersion  phenomena.   It  is  evident  that  a  high  degree 
of  measurement  accuracy  is  required  to  resolve  points  on  the  a  disper- 
sion arc.   On  this  arc,  the  reactance  is  usually  less  than  one  fifth 
as  large  as  the  resistance  with  corresponding  phase  angles  of  less  than 
12  degrees.   In  addition,  the  resistance  usually  changes  by  less  than 
20  percent  over  the  arc.   The  accuracy  required  to  measure  a   dispersion 
is  difficult  to  achieve  in  vitro,  but  it  is  especially  difficult  to 
achieve  in  situ. 

In  vitro,  the  a  and  8  dispersions  are  both  sensitive  to  age 
and  environmental  factors.   The  a  dispersion  appears  to  be  more  closely 
related  to  in  vivo  properties  of  the  tissue  and  breaks  down  rapidly 
after  excision.   The  £  dispersion  may  not  change  significantly  for  many 
hours  after  excision.   It  disappears  only  after  cellular  structure 
completely  deteriorates.   A  detailed  discussion  of  B  and  g  dispersion 
is  beyond  the  scope  of  the  present  work.   However,  a  brief  introduction 
to  the  mechanisms  which  underlie  these  phenomena  is  of  interest  and  is 
given  in  the  following  sections. 

3  dispersion 

The  explanation  of  g  dispersion  follows  from  principles  of 
classical  electric  field  theory  first  presented  by  Maxwell  in  1873. 
The  theory  has  been  applied  to  suspensions  of  conducting  spherical, 
elipsoidal  and  cylindrical  particles  surrounded  by  poorly  conducting 
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shells  (membranes)  to  simulate  various  biological  materials.   In  brief, 
3  dispersion  is  a  structural  phenomenon  depending  on  the  distribution 
of  the  membranes  and  conducting  fluids,  the  resistivity  of  the  fluids, 
and  the  dielectric  properties  of  the  membrane  material. 

At  frequencies  just  above  the  3  arc,  the  membrane  impedance  is 
negligible  in  comparison  to  the  effective  resistance  of  the  cell 
interior,  and  the  tissue  behaves  like  a  suspension  of  conducting  bodies 
in  a  similar  conducting  medium.   As  the  3  arc  is  entered,  the  membrane 
impedance  grows  comparable  to  the  effective  resistance  of  the  cell 
interior.   Toward  the  low  frequency  end  of  the  arc,  the  membrane 
impedance  has  grown  much  larger  than  the  effective  resistance  of  the 
cell  interior  and  the  tissue  begins  to  approach  (but  does  not  necessarily 
attain)  a  condition  in  which  the  cells  appear  like  empty  holes  in  the 
extracellular  fluid. 

While  3  dispersion  is  primarily  associated  with  the  cell 
membrane  proper,  the  membraneous  structures  within  the  cell  may  also 
contribute  to  the  dispersion  phenomena  by  a  similar  mechanism.   However, 
these  effects  are  smaller  and  tend  to  occur  at  higher  frequencies  (1 
MHz  to  10  MHz).   Dispersions  due  to  proteins  and  other  macromolecular 
matter  in  the  fluids  also  occur  in  the  MHz  range,  but  these  mechanisms 
are  distinctly  different  from  that  of  the  3  dispersion. 

ct  dispersion 

The  additional  resistance  rise  associated  with  a  dispersion 
originates  from  one  or  more  mechanisms  which  cannot  follow  the  rapid 
electrical  changes  that  occur  at  frequencies  above  the  a  dispersion 
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range.   Hence,  this  dispersion,  as  3  dispersion,  is  characterized  by 
capacitive  properties  (negative  reactance).   But,  contrary  to  3  disper- 
sion, the  underlying  phenomenon  does  not  necessarily  involve  the  charge 
and  discharge  of  adjacent  faces  of  a  dielectric  material.   At  least 
five  different  mechanisms  have  been  proposed  to  explain  a  dispersion 
(Schwan,  1957;  Schwan  and  Cole,  1960;  Fatt,  1963;  Ranck,  1963b).   These 
are  summarized  in  the  following  paragraphs. 

Membrane  permeability. --The  term  permeability  refers  to  the 
physical  transport  of  ions  through  the  cell  membrane  and  is  envisioned 
to  occur  at  specific  sites  in  the  membrane.   If  the  net  ion  current  is 
regulated  by  the  cell  such  that  it  is  related  to  the  membrane  voltage 
through  a  differential  equation  in  time,  then  the  permeability  is 
electrically  equivalent  to  an  impedance  characterized  by  resistive  and 
capacitive  elements,  even  though  no  physical  capacitance  exists. 
Inductive  elements  are  also  possible.   Evidence  for  the  permeability 
mechanism  is  based  mainly  on  studies  of  large  nerve  axons  (Cole,  1955; 
Schmitt,  1955). 

Membrane  structure . — The  membrane  is  envisioned  as  a  double 
layered  structure  of  electrically  different  materials  (lipid  and  protein) 
This  structure  results  in  a  frequency  dispersion  of  the  total  membrane 
impedance,  which  is  characterized  by  a  bounded  increase  in  the  effective 
parallel  resistance  and  capacitance  of  the  membrane  as  the  frequency 
is  decreased.   Schwan  (1957)  has  predicted  that  this  dispersion  should 
occur  in  the  frequency  range  of  10  Hz  to  1  kHz,  which  is  in  agreement 
with  the  typical  range  for  a  dispersion. 
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Intracellular  channels .— Fatt  (1963)  has  developed  a  theory  of 
dispersion  based  on  the  structural  organization  of  the  muscle  fibre. 
He  considers  two  parallel  current  paths.   One  is  formed  by  the  bulk 
impedance  of  the  cell  membrane  in  series  with  the  bulk  intracellular 
fluid  and  is  primarily  responsible  for  0  dispersion.   The  other  path 
originates  at  specific  sites  (pores)  in  the  cell  membrane  and  continues 
through  a  system  of  intracellular  channels  (the  sarcotubular  system 
of  the  endoplasmic  reticulum).   It  is  primarily  responsible  for  a 
dispersion.   For  convenience,  let  the  parallel  paths  be  designated  by 
a  and  0  corresponding  to  their  respective  dispersion  phenomena. 

Each  path  is  represented  by  a  resistor  in  series  with  a  capacitor. 
Since  the  channels  making  up  the  a   path  collectively  occupy  a  minute 
cross  section  of  the  membrane  surface  and  cell  interior,  the  effective 
resistance  of  the  a  path  is  much  higher  (order  of  100  times)  than  that 
of  the  0  path.   On  the  other  hand,  a  special  mechanism  is  postulated 
for  the  channel  entrance  which  permits  the  effective  capacitance  of  the 
a  path  to  be  much  larger  (order  of  25  times)  than  that  of  the  9  path. 
Hence,  at  high  frequencies,  the  a   path  appears  as  a  high  resistance 
across  the  0  path  and  has  little  effect  on  the  0  dispersion.   At  low 
frequencies,  the  0  path  appears  as  a  large  reactance  across  the  a  path 
and  has  little  effect  on  the  a   dispersion.   The  a  dispersion  is  much 
smaller  than  the  0  dispersion  because  the  effective  impedance  of  the 
extracellular  fluid  that  shunts  the  cell  remains  constant,  while  the  a 
path  phenomenon  occurs  at  a  much  higher  impedance  level  than  that  of 
the  0  path. 

With  a  constant  resistor  and  capacitor  representing  the  two 
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paths,  the  a  and  3  arcs  are  centered  on  the  real  axis,  corresponding 
to  one  time  constant  each.   In  practice,  the  arcs  are  centered  below 
the  axis,  corresponding  to  a  distribution  of  time  constants.   To 
account  for  this  fact,  Fatt  postulates  and  investigates  the  existence 
of  distributed  electrical  coupling  between  the  channel  system  and  the 
bulk  intracellular  fluid. 

Cable  properties.— Ranck  (1963b)  has  developed  a  theory  of  a 
dispersion  based  on  the  cable-like  properties  of  cell  processes  in 
cortical  tissue.   He  postulates  that  at  low  frequencies  (below  5  kHz), 
the  signal  current  is  constrained  by  the  cell  membranes  to  flow 
primarily  along  the  cell  processes  (axons  and  dendrites).   The  tissue 
then  appears  as  a  system  of  cables  consisting  of  intracellular  and 
extracellular  fluid  conductors  separated  by  a  leaky  dielectric  material 
(the  cell  membrane).   Using  the  cable  equation  and  certain  assumptions 
to  account  for  geometrical  and  material  aspects  of  the  tissue,  Ranck 
shows,  in  quantitative  terms,  that  the  cable  mechanism  may  well  give 
rise  to  the  experimentally  observed  a  dispersion  (Ranck,  1963a). 

At  higher  frequencies  (above  5  kHz),  the  membrane  impedance  has 
dropped  sufficiently  so  that  the  signal  current  is  no  longer  constrained 
to  paths  along  the  cell  processes.   The  current  then  flows  primarily  in 
directions  transverse  to  the  cell  membranes.   The  mechanism  for  $  disper- 
sion occurs  according  to  the  principles  of  Maxwell  as  previously 
discussed. 


Relaxation  of  ionic  atmospheres . --This  mechanism,  suggested  by 
Schwan,  is  based  on  impedance  studies  of  suspensions  of  small  non- 
conducting particles  (polystyrene,  glass,  kaolin,  etc.)  in  electrolyti 
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media.   These  suspensions  give  rise  to  low  and  high  frequency  dispersion 
phenomena  which  are  similar  to  the  a  and  g  dispersion  of  tissue.   At 
low  frequencies,  the  particles  are  perfectly  insulating.   Hence,  the 
low  frequency  dispersion  can  only  arise  through  a  variation  in  the 
effective  impedance  of  the  current  paths  between  the  particles.   The 
mechanism  lies  in  the  displacement  from  equilibrium  position  of  the 
ionic  layers  which  are  known  to  collect  around  particles  in  electrolyte 
suspension.   The  displacement  is  accompanied  by  a  change  in  the  stored 
electrostatic  energy  of  the  layers  and  a  purturbation  of  the  field 
around  the  particle.   These  phenomena  are  frequency  dependent  and  are 
reflected  as  reactive  and  resistive  components  in  the  suspension 
impedance . 

Many  biologically  significant  phenomena  have  been  reported  to 
influence  the  impedance  properties  of  certain  tissues  in  vivo.   In 
view  of  these  findings  and  the  possible  mechanisms  that  underlie 
impedance  phenomena,  the  measurement  and  analysis  of  tissue  impedance 
have  become  an  important  research  tool  (e.g.,  Adey  et  al,  1962; 
Aladjalova,  1964;  Cole,  1962;  Geddes  and  Baker,  1967;  Nevis  and  Collins, 
1967;  Shalit  and  Mahler,  1966;  Van  Harreveld,  1966).   Impedance  measure- 
ments are  also  of  potential  interest  in  clinical  applications.   Unfor- 
tunately, inadequacies  and  uncertainties  in  experimental  technique  have 
hindered  their  development  for  widespread  use.   Some  of  the  problems  are 
resolved  in  the  present  work. 

Literature  Review 

Three  principal  classes  of  electrode  systems  have  evolved  to 
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meet  the  varied  demands  of  biological  impedance  studies.   In  each  case, 
two  or  more  electrodes  (the  current  electrodes)  are  required  to  supply- 
signal  current  to  the  tissue  under  study.   The  classes  are  distinguished 
according  to  the  method  of  observing  the  voltage  signal.   Thus,  in  two- 
electrode  systems,  the  voltage  is  observed  directly  at  the  current 
electrodes;  in  three-electrode  systems,  the  voltage  is  observed  between 
a  potential  electrode  and  one  of  the  current  electrodes;  in  four- 
electrode  systems,  the  voltage  is  observed  between  two  potential 
electrodes . 

The  major  reason  for  using  three-electrode  and  four- electrode 
systems  is  that  the  potential  electrodes  nominally  do  not  draw  current 
so  that  the  observed  voltage  response  is  not  influenced  by  interface 
impedance.   Higher  order  systems  also  permit  greater  freedom  in 
selecting  the  current  field  configuration,  which  sometimes  leads  to 
more  versatile  experiments  and  better  anatomical  localization  of 
observed  impedance  phenomena.   Furthermore,  they  facilitate  control 
of  the  signal  current  level  and  are  especially  convenient  for  monitoring 
ongoing  biological  activity  during  the  impedance  measurement. 

To  provide  further  insight  into  the  nature  and  limitations  of 
biological  impedance  measurements,  some  of  the  major  works  in  this 
field  are  reviewed  in  the  following  sections.   An  effort  is  made  to 
indicate  historical  background  and  reveal  common  relationships  between 
the  three  classes  of  electrode  systems. 

Two-Electrode  Measurements 

Early  studies  of  the  conductivity  of  biological  materials  may 
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be  traced  as  far  back  as  the  nineteenth  century.   However,  the  earliest 
major  work  to  influence  modern  tissue  impedance  research  was  that  of 
Fricke  and  Morse  (1925),  who  first  applied  field  theory  to  the  inter- 
pretation of  experimental  impedance  data.   Thus,  from  measurements  of 
suspensions  of  red  blood  cells,  the  authors  were  able  to  estimate  the 
resistive  and  capacitive  properties  of  the  cell  membranes.   The  suspen- 
sions were  contained  in  an  electrolytic  cell  and  measured  by  means  of 
a  wheatstone  bridge.   The  electrolytic  cell  was  shaped  like  an  hour- 
glass with  large  platinized  platinum  electrodes  sealed  into  either 
end.   Although  this  design  was  selected  to  minimize  electrode  polariza- 
tion effects,  the  authors  were  well  aware  that  the  low  frequency 
utility  of  the  cell  was  still  quite  limited.   The  wheatstone  bridge 
was  operable  from  800  Hz  to  4.5  MHz,  but  no  results  were  reported  for 
frequencies  below  3600  Hz. 

A  large  body  of  impedance  work  has  been  done  on  excised  tissues 
as  well  as  cell  suspensions,  using  methods  similar  to  that  of  Fricke  and 
Morse.   Various  improvements  have  appeared  with  respect  to  the  bridge 
system  (Cole  and  Curtis,  1937-  Schwan  and  Sittel,  1953).   In  addition, 
a  number  of  techniques  have  evolved  for  extending  the  low  frequency 
capability  of  the  electrolytic  cell  (Schwan,  1963).   The  three  principal 
methods  are  summarized  below. 

Distance  variation  method.— In  the  distance  variation  method 
(Fricke  and  Curtis,  1937),  measurements  are  made  at  two  different 
electrode  separations.   By  subtracting  the  smaller  from  the  larger  ■ 
result,  polarization  effects  are  eliminated  and  the  impedance  of  the 
difference  volume  is  obtained. 
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Substitution  method. — In  the  substitution  method  (Cole  and 
Curtis,  1937),  the  sample  with  unknown  impedance  properties  is  replaced 
by  a  sample  with  known  impedance  properties .   The  electrode  polariza- 
tion impedance  is  then  calculated  from  measurements  on  the  known  sample 
and  subtracted  from  measurements  on  the  unknown  sample  to  give  the 
desired  impedance.   Unfortunately,  the  method  is  not  completely  straight- 
forward because  polarization  impedance  is  influenced  by  the  presence 
of  cells  in  the  sample  (Schwan,  1963).   One  way  to  avoid  this  discrepancy 
is  to  connect  the  electrodes  indirectly  to  the  sample  through  a  fixed 
intermediate  electrolyte  (Schwan,  1954;  and  Fatt,  1963). 

Frequency  variation  method. — Theoretical  (e.g.,  Fricke,  1932)  and 
experimental  (e.g.,  Smith,  et  al.,  1967)  studies  have  shown  that  electrode 
polarization  impedance  may  be  approximated  as  a  negative  power  function 
of  frequency  with  a  constant  phase  angle,  i.e.,  Z  =  Mf~ae~Je,  where  M, 
a,  and  9  are  real  positive  constants  and  a  <  1,  9  <  90  degrees. 
Several  schemes  based  on  this  fact  have  been  given  to  distinguish 
between  electrode  polarization  and  Tissue  impedance  (Schwan,  1963;  Smith, 
et  al.  ,  1967).   In  essence,  the  polarization  impedance  is  first  measured 
at  low  frequencies  where  the  sample  impedance  is  negligible.   These 
data  are  used  to  estimate  the  polarization  impedance  at  higher 
frequencies  so  that  it  may  be  discounted  from  the  sample  measurements. 

In  all  of  the  techniques  for  eliminating  the  effects  of  electrode 
polarization,  the  sample  impedance  is  estimated  as  the  difference 
between  two  quantities  which  are  assumed  to  share  a  common  polarization 
impedance.   As  the  frequency  is  decreased,  the  common  polarization 
impedance  increases  with  respect  to  the  sample  impedance.   Hence,  a 
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small  relative  error  in  the  quantities  to  be  differenced  leads  to  an 
increasingly  larger  error  in  the  result.   In  practice,  a  satisfactory 
degree  of  accuracy  can  only  be  achieved  at  frequencies  where  the  common 
impedance  is  less  than  the  impedance  to  be  estimated. 

In  1950,  the  status  of  in  situ  impedance  measurements  was 
reviewed  (Benjamin,  et  al.,  1950)  and  many  of  the  technical  difficul- 
ties involved  in  such  measurements  were  discussed.   The  authors 
concluded  that  reliable  in  situ  data  were  not  yet  available  at  physio- 
logical frequencies  and  called  for  greater  effort  in  that  direction. 
For  obvious  reasons,  methods  of  correcting  for  electrode  polarization 
are  more  difficult  to  implement  in  situ  than  in  vitro.   Few,  if  any, 
two-electrode  studies  in  situ  are  beyond  serious  criticism. 

The  earliest  of  the  more  precise  in_  situ  studies  was  reported 
by  Schwan  and  Kay  (1956,  1957).   They  measured  resistive  and  capacitive 
properties  of  various  thoracic  tissues  in  anesthesized  dogs.   The 
electrode  system  consisted  of  a  cylindrical  probe  with  one  electrode 
forming  the  tip  and  one  electrode  forming  part  of  the  shank.   The 
electrode  material  was  platinized  platinum  and  the  surface  area  of  each 
electrode  was  0.3  cm  .   Measurements  were  performed  with  a  wheatstone 
bridge  at  frequencies  from  10  Hz  to  10  kHz. 

To  correct  for  electrode  polarization,  the  polarization  impedance 
was  measured  as  a  function  of  frequency  in  isotonic  saline.   These 
results  were  modified  by  a  masking  factor  to  estimate  the  polarization 
impedance  in  the  tissue  at  each  frequency.   The  masking  factor  was 
assumed  to  be  independent  of  frequency  and  was  determined  by  comparing 
the  in-saline  and  in-tissue  polarization  impedances  at  10  Hz.   Special 
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precautions  were  necessary  to  assure  that  the  fragile  electrode  surfaces 
were  not  altered  between  saline  and  tissue  measurements.  The  accuracy 
of  the  series  polarization  resistance  and  capacitance  calculated  for 
in-tissue  conditions  was  stated  as  about  ±10  percent.   The  error  in 
the  result  for  tissue  capacitance  was  given  as  ±10  percent  at  100  Hz 
and  ±100  percent  at  10  Hz.   In  other  words,  the  reactive  part  of  the 
polarization  impedance  was  equal  to  that  of  the  tissue  at  100  Hz  and 
was  ten  times  that  of  the  tissue  at  10  Hz.   For  the  tissue  resistance, 
the  error  was  given  as  ±3  percent  at  10  Hz.   The  higher  accuracy 
reflects  the  fact  that  at  10  Hz  the  tissue  resistance  is  much  larger 
(more  than  30  times)  than  the  tissue  reactance. 

Schwan  and  Kay  reported  that  the  resistance  of  certain  tissues 
rises  about  20  percent  beginning  at  a  frequency  below  1  kHz  and 
continuing  more  sharply  toward  lower  frequencies.   This  finding  and 
certain  of  the  values  found  for  tissue  resistance  have  since  been 
disputed  by  others  (Burger  and  van  Dongen,  1960-61;  Rush,  et  al,  1963), 
who,  on  the  basis  of  four-electrode  measurements,  claimed  the  resistivity 
to  be  independent  of  frequency  in  this  range.   Several  reasons  were 
proposed  by  Rush  et  al.  for  the  discrepancy.   Important  among  these 
were  erroneous  assumptions  in  the  correction  for  polarization  impedance. 
Further  evidence  (Schwan,  1963,  Figure  10;  Smith  et  al.,  1967)  indicates 
that  contrary  to  what  was  assumed,  the  masking  fao»tor  may  well  have  been 
dependent  on  frequency,  especially  below  500  Hz.   If  this  is  true,  then 
the  anomalous  change  in  slope  of  the  capacitance  curve  (Schwan  and  Kay, 
1957,  Figure  1)  at  frequencies  below  500  Hz  is  also  in  question.   It  has 
been  reasonably  well  established  that  a  small  variation  in  resistivity 
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does  occur  in  many  tissues  at  low  frequencies  (a  dispersion).   Therefore, 
in  view  of  the  many  uncertain  factors,  it  is  suspected  that  none  of 
the  conflicting  reports  are  completely  reliable. 

The  work  of  Schwan  and  Kay  illustrates  some  of  the  problems  and 
uncertainties  encountered  in  two-electrode  measurements  in  situ.   It 
probably  also  demonstrates  the  best  accuracy  that  can  be  achieved 
under  such  conditions-at  least  in  the  light  of  current  knowledge. 
Many  experiments  (especially  chronic  ones)  require  smaller  electrodes 
and  the  error  due  to  polarization  is  correspondingly  larger.   In  any 
case,  manual  and  computational  correction  procedures  are  time  consuming 
and  difficult  to  carry  out  with  the  degree  of  caution  necessary.   Further- 
more, the  validity  of  the  key  assumptions  on  which  the  accuracy  of  two- 
electrode  results  depend  is     difficult,  if  not  impossible,  to  verify 
in  specific  experiments. 

Three-Electrode  Measurements 

Three-electrode  systems  are  useful  where  there  is  extreme 
assymmetry  in  the  current  electrode  configuration.   They  have  been  used 
primarily  to  study  the  membrane  properties  of  individual  cells  in  situ 
and  in  vitro.   The  current  is  applied  between  a  microelectrode  (0.5  to 
5  micron  tip)  in  the  cell  body  and  a  large  reference  electrode  in  the 
external  medium.   It  passes  across  the  membrane  of  the  cell  body  and 
in  the  case  of  neurons,  also  through  the  dendritic  branches  (Rail,  1959). 
Almost  all  of  the  potential  difference  between  the  interior  of  the 
cell  body  and  the  reference  electrode  occurs  across  the  cell  membrane. 
The  reference  electrode  can  usually  be  designed  to  produce  a 
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negligible  polarization  component  in  the  observed  voltage.   However, 
the  impedance  of  the  microelectrode  is  quite  large  (megohms)  and 
cannot  be  neglected.   To  overcome  this  problem,  a  second  microelectrode 
may  be  added  in  the  cell  body  to  record  potential  (Fessard  and  Tauc, 
1956;  Combs  et.  al. ,  1955,  1959).   The  high  impedance  of  the  intra- 
cellular electrodes  makes  such  systems  especially  sensitive  to 
parasitic  coupling  between  electrodes.   Fessard  and  Tauc  used  completely 
separate  micropipettes  to  minimize  this  hazard.   For  smaller  cells, 
the  double-barrelled  micropipette  of  Combs  et.  al.  is  more  convenient, 
but  the  problem  of  parasitic  coupling  is  intensified. 

The  foregoing  methods  have  been  used  only  with  pulse  waveforms 
to  determine  membrane  time  constants  and  for  stimulating  and  observing 
active  responses  of  neurons.   A  systematic  application  of  sinusoidal 
techniques  would  be  desirable  both  from  the  standpoint  of  defining 
measurement  errors  and  obtaining  more  precise  information  about  the 
electrical  properties  of  the  cell  membrane. 

Three-electrode  systems  have  also  been  used  with  a  distance 
variation  technique  to  study  tissue  impedance  in  anesthetized  animals. 
(Van  Harreveld,  et.  al. ,  1963;  Nicholson,  1965).   The  current  field  is 
set  up  by  means  of  a  fixed  pair  of  current  electrodes  and  the  voltage 
is  measured  with  respect  to  one  current  electrode  as  a  function  of 
potential  electrode  position.   Regardless  of  the  particular  procedure 
and  field  configuration  used,  the  calculation  of  specific  impedance 
requires  the  evaluation  of  the  voltage  difference  between  probe  posi- 
tions.  Hence  this  method  is  subject  to  similar  limitations  with 
respect  to  the  ratio  of  common  to  difference  voltage  as  were  discussed 
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earlier  for  two-electrode  systems.   Another  limitation  of  this  technique 
is  that  the  measurements  at  each  electrode  position  are  made  at 
different  time  intervals. 


Four-Electrode  Measurements 


It  is  convenient  to  classify  four-electrode  measurements 
according  to  the  manner  in  which  the  potential  electrodes  are  used. 
In  this  section,  a  method  is  labeled  continuous  if  the  potential 
electrodes  take  part  in  the  impedance  measurement  simultaneously  and 
discontinuous  if  separate  measurements  are  required  at  each  electrode. 

Discontinuous  methods 

The  only  method  of  current  interest  which  deals  separately  with 
each  potential  electrode  is  a  bridge  balancing  scheme.   The  basic 
technique  was  originally  described  by  Shedlovsky  (1930a,  b)  for 
measuring  conductance  in  nonbiological  electrolytes.   Recently  it  was 
applied  to  biological  media  (Songster,  1967;  Hill,  unpublished)  and 
modified  to  include  capacitive  components  and  a  broader  frequency 
range  (5  Hz  to  200  kHz).   The  main  features  are  illustrated  in  Figure 
1.2.   The  electrode-tissue  system  is  represented  by  the  simplified 
equivalent  network,  Z  ,  Z  ,  and  Z   (see  Figure  2.8  for  the  complete 
equivalent  network).   The  current  electrode  interface  impedances  are 
included  in  Z^  and  Z  .   Z  represents  the  tissue  impedance  to  be 

determined.   The  bridge  balances  are  achieved  by  means  of  Z  and  Z  . 

12 

and  the  detector  may  be  switched  from  position  1  to  position  2  as  needed. 

The  measurement  procedure  begins  with  Z  set  at  zero  and  proceeds 
in  two  successive  steps. 
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Figure  1.2.   Four-Electrode  Bridge 


Step  1:   Adjust  Z  to  null  point  1, 
Step  2:   Adjust  Z  to  null  point  2. 


If  tj_  and  t2  denote  the  times  at  which  step  1  and  step  2  are  performed, 
then  the  balance  conditions  may  be  written 


at  tx:     Zx  ±  £l  +  Za(tx)  =  Zb(tx)  +  Zc(t1) 


at  t. 


Z2  1  e2  +  Zc(t2)  =  Zb(t2)  ♦  Za(t2)  ♦  2jL 
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where  the  e  are  uncertainties  in  the  adjustment  of  Z  and  Z  .  Note 

1      2 

that  each  condition  actually  represents  two  equations,  one  for  the 
resistive  parts  and  one  for  the  reactive  parts  of  the  impedances. 

The  balance  equations  may  be  added  and  the  result  arranged  tc 
give 

Zb(tl}  +  Zb(t2)  =  z2  +  <AZc  "  Aza>  *  (£l  ±   £2} 

where  AZQ  and  AZa  are  the  changes  in  ZQ  and  Za,  if  any,  that  occur 

from  time  t±   to  time  t2>   Ideally,  AZc,  AZa,  *1   and  e2   should  be  zero 

so  that  Z2/2  is  equal  to  Zb  (or  its  average  value  if  it  varies).   In 

practice,  these  deviations  are  finite  and  must  fall  within  limits  much 

smaller  than  Z^   if  the  method  is  to  be  successful. 

At  least  two  factors  are  of  major  concern  in  evaluating  the 

suitability  of  the  bridge  method  for  a  particular  biological  study: 

bridge  resolution  and  constancy  of  Z&,   Zb  and  ZQ  with  time.   At 

frequencies  below  100  Hz  (sometimes  higher),  Zfi  and  ZQ  may  be  of  the 

same  order  as  or  very  much  larger  than  Zb ,  depending  primarily  on  the 

electrode  design.   For  example,  with  large  platinized  platinum  electrodes 

in  vitro,  these  impedances  will  usually  be  on  the  same  order  as  or 

smaller  than  Zfc>.   They  will  also  tend  to  remain  constant  in  the  t  t 

interval.   However,  with  small  stainless  steel  electrodes,  in  situ, 

Za  and  Zc  may  well  be  1000  or  more  times  as  large  as  Zb  and  are  likely 

to  vary  during  the  measurement  interval.   The  first  example  is  easily 

within  bridge  capabilities,  but  the  electrode ' system  requires  careful 

preparation.   In  the  second  case,  very  high  bridge  resolutions  (greater 

than  1  part  in  105)  would  be  required  to  make  the  e  much  smaller  than  Z. 

b 
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Even  then,  the  slightest  variations  in  Z  and  Z  with  time  could 

3.         C 

produce  errors  greater  than  100  percent.   The  problem  is  compounded 

by  the  fact  that  as  the  frequency  is  decreased,  it  usually  takes  longer 

to  achieve  a  bridge  balance. 

Clearly,  there  is  a  whole  range  of  experiments  with  relative 
impedance  values  lying  between  the  extremes  discussed  above.   For  many 
of  these,  the  bridge  method  is  not  suitable  on  the  basis  of  resolution 
and  constancy  arguments  alone.   On-line  bridge  methods  are  also  limiting 
in  that  only  one  frequency  at  a  time  may  be  measured  and  each  measure- 
ment takes  several  or  more  minutes  to  perform.   Hence,  impedance  cannot 
be  monitored  continuously,  nor  can  the  instantaneous  impedance  spectrum 
be  determined. 

An  advantage  of  the  bridge  method  is  that  at  balance,  the  potential 
electrode  in  use  is  approximately  at  ground  potential.   Hence,  the  input 
impedance  of  the  null  amplifier  need  not  be  especially  large.   This 
feature  is  shared  with  the  method  presented  in  Chapters  4  and  5  of  the 
present  work. 

Continuous  methods 

The  use  of  a  four-electrode  method  in  biological  impedance 
measurements  was  first  reported  by  Burger  and  van  Milaan  (1943).   They 
measured  various  segments  'of  the  human  body  and  used  the  results  to 
estimate  the  specific  resistance  of  certain  body  tissues.   The  current 
was  supplied  in  manually  switched  d.c.  pulses  from  a  battery  in  series 
with  a  large  fixed  resistor.   To  find  the  transfer  resistance  of  the 
electrode-tissue  system,  the  voltage  between  potential  electrodes  was 
compared  to  the  voltage  across  a  known  resistor  connected  in  series 
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with  the  current  electrodes.   The  voltages  were  amplified  by  means  of 
a  single-ended  amplifier  and  read  on  a  galvanometer.   Since  the  current 
source  and  the  electrode-tissue  system  were  not  independently  grounded, 
the  amplifier  could  be  switched  at  will  to  the  terminal  pairs  of 
interest. 

Rush  et.  al.  (1963)  used  a  method  similar  to  that  of  Burger 
and  van  Milaan  to  generate  the  current  waveform.   However,  the  voltage 
response  was  sensed  by  means  of  a  differential  amplifier.   The 
amplifier  common  (ground)  was  connected  to  a  remote  point  on  the  animal. 
The  common  mode  signal  of  this  configuration  was  of  the  same  order  of 
magnitude  as  the  differential  signal.   Hence,  the  requirements  on  the 
input  impedance  and  common  mode  rejection  ratio  of  the  amplifier  were 
not  difficult  to  satisfy. 

For  reasons  of  convenience,  reliability  and  accuracy,  it  is 
usually  desirable  to  use  ground  as  a  common  reference  for  the  current 
electrode  and  potential  electrode  circuits.   This  is  especially  true 
when  sophisticated  measurement  techniques  are  employed  and  battery 
operation  becomes  impractical.   The  methods  discussed  in  the  remainder 
of  this  section  use  electronic  current  generators  and  common  ground 
connections. 

Among  the  earliest  applications  of  four-electrode  methods  were 
the  now  famous  voltage  clamp  experiments  of  Hodgkin,  Huxley,  and  Cole. 
In  the  original  experiments  (Cole,  1949;  Marmont ,  1949),  the  membrane 
voltage  was  sensed  with  the  same  pair  of  electrodes  that  were  used  to 
provide  the  membrane  current  (a  two-electrode  method).   Later,  Hodgkin, 
Huxley,  and  Katz  (1952)  eliminated  the  error  introduced  by  the  current 


27 


electrode  polarization  voltages  and  some  of  the  intervening  medium 
by  adding  a  second  pair  of  electrodes  to  measure  the  membrane  voltage 
(a  four-electrode  method). 

The  basic  voltage  clamp  scheme  has  been  carried  out  in  various 
ways.  One  of  the  later  arrangements  used  by  Moore  and  Cole  (1963)  is 
depicted  in  simplified  form  in  Figure  1.3.  The  potential  amplifier 


Figure  1.3.   Voltage  Clamp  Scheme 


senses  the  membrane  voltage  differentially  by  means  of  two  micropipette 
electrodes.   The  output,  V  ,  is  compared  to  the  control  signal,  V  , 
in  the  control  amplifier  which  adjusts  the  current  of  the  internal 
current  electrode  to  maintain  V  =  V  .   The  external  current  electrode 

o     *»• 
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consists  of  three  separate  cylindrical  pieces.   Current  is  measured 
only  in  the  center  piece  where  the  field  is  uniform  and  calculable. 
The  input  of  the  current  sensing  circuit  is  essentially  at  ground 
potential.   The  external  current  electrodes  in  this  system  have  a 
relatively  large  area  and  are  usually  platinized  platinum  (sometimes 
Ag-AgCl). 

The  common  mode  signal  to  the  potential  amplifier  is  the 
voltage  occurring  between  the  external  potential  electrode  and  ground. 
The  differential  signal  is  essentially  the  membrane  voltage.   The 
highest  ratio  of  common  mode  to  differential  signal  occurs  at  the  peak 
of  the  membrane  current.   A  typical  ratio  is  20  mv/60  mv  or  1/3  (Moore 
and  Cole,  1963,  Figure  10).   Hence,  for  a  1  percent  common  mode  error, 
the  rejection  ratio  of  the  potential  amplifier  need  only  be  about  33, 
which  is  not  difficult  to  achieve.   The  low  common  mode  signal  found 
here  is  in  contrast  to  the  relatively  high  ratios  that  occur  in  a 
large  number  of  four-electrode  applications,  especially  those  concerned 
with  in  situ  tissues.   Higher  common  mode  to  differential  ratios  not 
only  require  higher  common  mode  rejection  ratios,  but  also  higher  input 
impedances  in  the  amplifiers.   These  topics  are  discussed  in  detail 
along  with  other  problems  in  chapters  to  follow. 

In  1955,  Freygang  and  Landau  reported  a  clever  four-electrode 
scheme  to  measure  the  specific  resistance  of  the  cerebral  cortex  in 
anesthetized  animals.   One  current  electrode  was  placed  in  a  well  of 
physiological  fluid  overlying  an  exposed  section  of  cerebral  cortex. 
The  other  electrode  was  placed  in  the  pharynx.   In  this  way,  a  fairly 
uniform  current  density  was  established  in  the  cortex  under  the  well. 
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The  current  electrodes  were  made  of  platinized  platinum  and  had  fairly- 
large  contact  areas  to  minimize  polarization  impedance.   The  voltage 
across  the  cortex  was  sensed  by  means  of  two  micropipette  electrodes 
which  were  connected  to  a  differential  amplifier.   The  applied  current 
was  a  square  wave  derived  from  a  phase  splitter  with  an  output  impedance 
greater  than  20  kilohms.   Pulse  durations  of  0.3  ms  to  0.7  ms  were  used. 

In  a  typical  run,  the  cortical  voltage  was  20  mv  or  less,  while 
the  voltage  between  current  electrodes  was  about  1.8  volts.   The 
common  mode  rejection  ratio  of  the  differential  amplifier  was  given  as 
2,000.   Although  the  current  was  derived  from  a  balanced  source,  the 
common  mode  signal  to  the  amplifier  was  not  necessarily  inconsequential 
because  of  probable  differences  in  the  impedance  of  the  current  paths 
above  and  below  the  region  of  the  potential  electrodes.   To  illustrate 
this  point,  suppose  the  common  mode  signal  was  one  fourth  the  voltage 
between  current  electrodes  or  approximately  0.45  volts.   Then  the 
common  mode  rejection  error  would  be  somewhat  greater  than  one  percent 
of  the  cortical  voltage.   The  fundamental  frequencies  of  the  applied 
square  waves  were  greater  than  700  Hz.   On  the  basis  of  amplifier 
rejection  ratio  alone,  it  is  unlikely  that  reliable  measurements  could 
have  been  made  at  much  lower  or  much  higher  frequencies  with  the  system 
described. 

Most  of  the  four-electrode  studies  published  to  date  have 
occurred  since  1960.   In  the  first  of  these  works,  Burger  and  van  Dongen 
(1960-1961)  modified  the  original  four-electrode  method  of  Burger  and 
van  Milaan  (1943)  to  include  sinusoidal  frequencies  from  20  Hz  to  5  kHz. 
The  basic  scheme  is  shown  in  Figure  1.4,  where  the  rectangle  with  four 
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Figure  1.M-.   Four-Electrode  Comparison  Sch 


eme 


terminals  represents  the  electrode-tissue  system.   The  procedure  was 
to  adjust  R  until  the  amplitude  of  Vg  was  equal  to  that  of  Vp .   Then 
the  magnitude  of  the  transfer  impedance  of  the  four-electrode  system 
would  be  approximately  R  divided  by  the  gain  of  the  differential 
amplifier.   The  common  mode  signal  to  the  amplifier  includes  one 
current  interface  voltage  as  well  as  VR  and  some  contribution  from  the 
tissue.   This  signal  is  likely  to  have  been  much  greater  than  the 
differential  signal  in  some  of  the  experiments  reported.   The  authors 
make  almost  no  mention  of  system  performance  or  measurement  difficulties 
encountered.   Also  no  technical  data  on  the  apparatus  were  given.   This 
omission  is  unfortunate  in  view  of  the  significant  differences  noted 
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between  their  results  and  those  of  Schwan  and  Kay  (1956)  which  this 
investigation  was  meant  to  clarify. 

Ranck  (1963a)  used  an  arrangement  similar  to  that  shown  in  Figure 
1.4  to  measure  transfer  impedance  in  rabbit  cerebral  cortex  in  acute 
experiments.   The  voltage  and  current  signals  were  connected  to  the 
vertical  and  horizontal  sweeps  of  an  oscilloscope,  respectively.   The 
magnitude  and  phase  of  the  impedance  were  determined  from  the  parameters 
of  the  resulting  eliptical  pattern  observed  on  the  oscilloscope  screen. 
The  electrode  system  included  three  small  platinum  electrodes  (20  to 
40  micron  diameter)  arranged  in  a  linear  array  on  the  exposed  cortical 
surface.   These  correspond,  in  order,  to  the  upper  three  electrodes 
of  Figure  1.4.   The  fourth  electrode  was  located  remotely  on  the 
animal.   The  cortical  electrodes  were  spaced  at  least  300  microns  apart, 
but  the  array  never  exceeded  1.3  mm  in  length. 

Ranck ' s  highly  assymmetrical  electrode  system  is  reminiscent  of 
the  arrangements  used  for  three-electrcde  measurements  on  individual 
cells  and  possesses  similar  advantages.   The  polarization  impedance  of 
the  remote  electrode  can  easily  be  made  negligible,  and  the  voltage 
gradient  decreases  rapidly  away  from  the  cortical  array.   Hence,  the 
common  mode  signal  to  the  differential  amplifier  is  of  the  same  order 
as  the  differential  signal.   The  common  mode  rejection  ratio  of  the 
amplifier  was  given  as  at- least  1,000,  which  in  this  case  appears  quite 
satisfactory.   Requirements  on  amplifier  input  impedance  were  also  not 
too  difficult  to  meet. 

The  most  informative  discusssion  of  errors  in  a  four-electrode 
impedance  investigation  was  given  by  Ranck  in  a  later  paper  (Ranck, 
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1966).   This  paper  also  was  first  to  report  the  application  of  four- 
electrode  technique  to  chronic  impedance  studies  in  unrestrained  animals. 
The  basic  signal-sensing  scheme  was  as  shown  in  Figure  1.4,  except 
that  the  current  monitoring  resistor  was  replaced  by  an  operational 
amplifier  circuit  with  negligible  input  impedance  (similar  to  that 
shown  in  Figure  1.3).   As  before,  the  voltage  and  current  signals  were 
displayed  as  a  Lissajous  pattern  on  an  oscilloscope.   In  addition,  they 
were  connected  to  a  special  instrument,  the  JB-5  lock-in  amplifier, 
which,  with  auxiliary  circuits,  was  calibrated  to  display  the  in- 
phase  and  quadrature  components  of  the  tissue  voltage  response  with 
respect  to  the  applied  current. 

The  lock-in-amplifier  uses  tuned  circuits  and  special  averaging 
techniques  which  greatly  increase  the  signal-to-noise  ratio  of  the 
displayed  signal.   Hence,  Ranck  was  able  to  measure  voltage  responses 
of  the  order  of  20  microvolts  in  the  presence  of  EEG  activity  of  the 
order  of  300  microvolts.   This  is  not  feasible  with  the  Lissajous 
pattern  method  unless  the  signal  frequency  is  far  removed  from  the 
noise  frequencies  and  a  noise  filter  is  used.   Low  current  densities 
in  the  tissue  (hence  small  voltage  responses)  are  desirable  to  minimize 
the  possible  influence  of  the  test  current  on  the  tissue  impedance. 
This  effect  has  yet  to  be  ascertained. 

The  main  electrode  -assembly  was  chronically  implanted  in 
selected  sites  of  the  rat  brain.   The  submerged  part  consisted  of  four 
parallel  wires  (75  micron)  spaced  uniformly  in  a  linear  array  1  to  1.5 
mm  across.   Each  wire  was  roughly  1  mm  long  and  was  insulated  to  within 
about  100  microns  of  the  tip.   The  wires  were  made  of  platinum  and  the 
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tips  were  platinized  to  reduce  polarization  impedance.   Two  field 
configurations  were  used.   In  the  linear  configuration,  the  brain 
electrodes  were  connected  in  the  order  shown  in  Figure  1.4,   In  the 
radial  configuration,  a  remote  electrode  was  used  as  the  current  return 
in  place  of  an  end  head  electrode.   The  latter  configuration  did  not 
possess  the  advantage  of  a  modest  common  mode  signal  as  found  in  the 
radial  configuration.   Ranck  made  reference  to  this  fact  and  noted 
greater  difficulties  and  poorer  results  with  the  linear  configuration. 

Ranck 's  discussion  of  errors  has  been  a  source  of  inspiration 
for  some  of  the  more  general  treatments  given  in  the  present  work. 
In  particular,  he  called  attention  to  the  errors  caused  by  parasitic 
resistance  (leakage)  and  capacitance  which  provide  undesirable  coupling 
between  the  current  and  potential  electrode  conductors .   In  the  linear 
configuration,  Ranck  found,  for  example,  that  at  frequencies  below  100 
Hz,  the  parasitic  resistance  had  to  be  greater  than  108  ohms  to  realize 
a  voltage  sensing  error  of  less  than  1  percent.   For  frequencies  below 
about  10  Hz,  the  requirement  became  greater  than  10   ohms.   He  noted 
the  difficulty  of  maintaining  resistances  of  this  magnitude  on  the  head 
of  animals  in  chronic  experiments.   These  requirements  also  came 
uncomfortably  close  to  the  nominal  input  resistance  of  the  response 
amplifier,  which  was  given  as  greater  than  10  ohms. 

Ranck  neglected  the  phase  angle  of  the  electrode  polarization 
impedance  in  his  error  analysis.   But,  he  did  make  use  of  this 
phenomenon  in  testing  the  electrode  assemblies  before  implantation. 
Hence,  each  assembly  was  checked  for  leakage  in  0.2  percent  saline  at 
1.5  Hz  and  the  phase  angle  of  the  voltage  response  was  taken  as  the 


34 


most  sensitive  indicator  of  error.   The  phase  angle,  if  any,  was  due 
to  the  effect  of  electrode  polarization  capacitance.   The  reliability 
of  the  results  obtained  with  this  method  of  electrode  testing  is 
questionable  for  reasons  explained  in  the  following  paragraphs . 

As  Ranck  indicated,  the  error  due  to  a  given  amount  of  leakage 
between  a  current  electrode  and  a  potential  electrode  is  proportional 
to  the  product  of  the  electrode  impedances  divided  by  the  transfer 
impedance  to  be  measured.   Mathematically 
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where  Z  ,  Zp  and- Z  designate  the  current  electrode,  potential 
electrode,  and  transfer  impedances,  respectively.   At  the  test  frequency 
of  1.5  Hz,  Z_  and  Zp  were  largely  determined  by  interface  impedance. 
The  interface  impedance  is  not  a  sensitive  function  of  saline  concentra- 
tion, since  it  depends  mainly  on  a  hydrogen  transfer  mechanism  with 
the  water.   Therefore,  |ZQ[  and  |Zp|  would  be  expected  to  have  much 
higher  values  in  the  tissue  than  in  the  saline  test  solution  due  to 
the  masking  effect  of  cellular  components  (Schwan,  1963). 

The  0.2  percent  saline  used  for  the  electrode  tests  produced 
approximately  the  same  |Z„|  as  was  found  in  the  tissue.   However, 
since  |Z  I  and  |Z  |  were  probably  much  lower  in  the  saline  than  in  the 
tissue,  the  simulated  error  would  likewise  be  expected  to  be  much 
smaller,  but  by  the  same  factor  squared.   Suppose  that  isotonic  saline 
(0.9  percent)  were  used  for  the  electrode  tests.   Then,  the  ratio  of  IZ  I 
to  |Z  |  in  the  saline  would  more  closely  simulate  that  found  in  the 
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tissue  (assuming  the  cell  masking  factor  was  roughly  proportional  to 
the  ratio  of  tissue  to  tissue  fluid  impedance).   In  this  case,  the 
simulated  error  would  still  be  too  low,  but  only  by  the  factor  in  |zp| 
to  the  first  power.   The  foregoing  considerations  are  rather  significant 
because  they  imply  that  electrode  test  results  obtained  in  pure  saline 
do  not  provide  a  reliable  indication  of  leakage  errors. 

The  choice  of  test  solution  is  also  important  with  respect  to 
errors  caused  by  imperfect  common  mode  rejection  in  the  differential 
amplifier.   The  ratio  of  common  mode  signal  to  differential  signal 
is  approximately  given  by  the  ratio  of  |Z  |  to  |Z_j.   As  previously 
indicated,  the  ratio  found  in  tissues  is  more  closely  simulated  in 
isotonic  saline  than  in  0.2  percent  saline.   In  the  latter  case,  the 
ratio  would  be  too  low  to  give  a  reliable  indication  of  the  common 
mode  rejection  error. 

The  phase  angle  of  the  voltage  response  was  also  the  most 
sensitive  indicator  of  error  at  higher  frequencies  where  polarization 
impedance  was  negligible  and  electrode  impedance  was  essentially 
resistive.   Here,  the  phase  angle,  if  any,  was  due  entirely  to  capacitive 
coupling  between  the  current  and  potential  electrode  conductors.   The 
input  capacitance  of  the  response  amplifier  was  given  as  30  picofarads 
which  is  equivalent  to  about  5  megohms  at  1  kHz.   This  reactance 
significantly  loaded  the  potential  electrodes  and,  in  the  radial  configu- 
ration, caused  a  phase  error  in  the  apparent  response  of  approximately 
9  degrees^ while  the  magnitude  error  was  less  than  2  percent. 

The  effects  of  parasitic  resistance  and  capacitance  were  greatly 
magnified  in  the  linear  configuration  by  the  high  ratio  of  common  mode 
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signal  to  differential  signal  which  varied  from  approximately  40  at  1 
kHz  to  600  at  1.5  Hz.   This  high  ratio  was  also  troublesome  with 
respect  to  the  common  mode  rejection  capability  of  the  amplifier.   To 
obtain  common  mode  error  signals  of  less  than  1  percent,  the  common 
mode  rejection  ratio  would  have  had  to  exceed  limits  ranging  from  4,000 
at  1  kHz  to  60,000  at  1.5  Hz.   The  common  mode  rejection  ratio  of  the 
system  used  was  given  as  only  better  than  1,000.   On  the  basis  of  the 
foregoing  discussion  and  in  concurrence  with  some  of  Ranck's  remarks, 
it  may  be  concluded  that  the  results  obtained  with  the  linear  configura- 
tion were  unreliable  at  nearly  all  frequencies.   Much  of  this  problem 
lay  in  the  extreme  demands  that  this  configuration  made  on  the  instru- 
mentation components.   Better  performance  could  have  been  achieved  with 
the  instrumentation  technique  described  in  Chapters  4  and  5  of  the 
present  work. 

Summary 

In  a  broad  sense,  electrode  polarization  leads  to  a  similar 
problem  in  every  type  of  electrode  system,  which  becomes  worse  as  the 
electrode  dimensions  are  reduced  and  the  frequency  is  lowered.   Briefly 
stated,  the  problem  is  to  find  an  unknown  quantity,  the  tissue  impedance, 
within  arbitrary  accuracy  limits,  by  taking  the  difference  between  two 
quantities  which  may  become  nearly  equal  at  low  frequencies.   The 
differencing  process  may  be  done  by  manual  and  computational  means 
(as  in  two-electrode  and  three-electrode  systems)  or  continuously  and 
automatically  by  electronic  means  (as  in  four-electrode  systems).   The 
experimental  evidence  indicates  that  the  former  methods  are  not 
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reliable  unless  the  common  quantity  is  significantly  smaller  than  the 
difference  quantity.   In  contrast,  the  four-electrode  electronic 
methods  are  reliable  and  capable  of  high  accuracy  even  when  the  common 
quantity  is  much  larger  than  the  difference  quantity.   This  subject  is 
considered  in  detail  in  Chapters  3  and  4  under  the  title  of  common 
mode  rejection. 

The  advantages  of  the  four-electrode  method  are  obtained  at  the 
cost  of  increased  complexity  and  additional  sources  of  error.   However, 
a  properly  designed  and  well-understood  four-electrode  experiment  can 
be  no  more  difficult  to  perform  than  a  two-electrode  experiment,  all 
factors  considered.   The  basis  for  this  statement  may  be  found  in  the 
remaining  chapters  of  this  report,  which  are  summarized  in  the 
following  paragraphs. 

Chapter  2. — Chapter  2  is  concerned  with  basic  analytical 
techniques  for  describing,  analyzing  and  designing  the  electrode- 
tissue  system.   Special  matrix  parameters  are  defined  for  the  system, 
and  the  concept  of  interface  impedance  is  clarified.   Certain  boundary 
effects  caused  by  interface  impedance  are  explained  and  techniques  for 
correcting  them  are  suggested.   Further  insight  is  provided  by  means 
of  a  convenient  equivalent  network,  which  is  derived  in  terms  of  the 
matrix  parameters  of  the  system.   To  conclude  the  chapter,  the  matrix 
parameters  are  derived  from  field  theory  for  two  different  electrode 
arrangements  and  illustrated  numerically. 

Chapter  3. — Chapter  3  is  concerned  with  the  most  common  method 
for  amplifying  the  signal  of  the  potential  electrodes.   Four  major 
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types  of  measurement  error  are  defined  and  analyzed.   Limit  expressions 
are  derived  for  the  parameters  of  the  instrumentation  system  and 
associated  paraphernalia.   The  results  are  illustrated  numerically, using 
the  examples  of  Chapter  2.   Besides  being  of  value  in  itself,  Chapter 
3  provides  important  definitions  and  background  for  Chapter  4. 

Chapter  4. --In  Chapter  4,  a  new  instrumentation  technique--the 
virtual-ground  system—for  extracting  the  signal  of  the  potential 
electrodes  is  introduced.   The  four  major  types  of  measurement  error 
defined  in  Chapter  3  are  analyzed  for  the  new  system.   Limit  expressions 
are  derived  for  the  system  parameters  and  numerical  results  are 
illustrated  using  the  examples  of  Chapter  2.   These  results  are 
compared  with  those  of  the  common  method  of  Chapter  3.   Two  additional 
topics  are  treated:   the  stabilization  of  the  virtual-ground  system  for 
broad  band  measurements  and  the  operation  of  the  negative  feedback 
filter  loop  for  noise  reduction. 

Chapter  5. --In  this  chapter,  an  experimental  system  using  the 
virtual-ground  configuration  is  described.   The  system  includes  a 
current  control  scheme  which  was  taken  from  the  literature  and  adapted 
for  broadband  impedance  measurements.   Several  experiments  are  reported 
to  verify  system  performance  and  indicate  the  quality  of  results  already 
obtained  in  investigations  on  live  animals. 


CHAPTER  2 

MATRIX  AND  NETWORK  DESCRIPTION  OF  FOUR-ELECTRODE  SYSTEMS 

Voltage  and  Current  Conventions 

It  is  convenient  to  assign  the  symbols  Ql,  PI,  P2,  and  Q2  to 
the  terminals  of  a  four-electrode  system,  where  Q  designates  the  current 
or  excitation  electrodes  and  P  designates  the  potential  or  response 
electrodes.   Let  the  complex  terminal  voltages  and  currents  describing 
the  state  of  the  system  be  as  defined  in  Figure  2.1. 
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Figure  2.1.   Voltage  and  Current  Conventions  for  the 
Electrode-Tissue  System 
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According  to  the  laws  of  network  theory 


V2  +  VP  +  Vl  "  VQ  =  ° 


(2.1) 


and 


XQ1  +  V  +  XP1  +  XP2  =  ° 


(2.2) 


which  shows  that  only  three  voltages  and  three  currents  may  be  defined 
independently.   Note  that  the  V  and  I  quantities  usually  possess  in- 
phase  and  quadrature  components  which  vary  with  frequency.   Equations 

(2.1)  and  (2.2)  each  imply  two  equations,  one  for  the  in-phase  and  one 
for  the  quadrature  components,  with  all  quantities  being  evaluated  at 
the  same  frequency. 

The  Impedance  Parameter  Matrix 

In  a  linear  system  with  sinusoidal  excitation,  the  complex 
voltages  and  currents  are  related  by  a  set  of  linear  simultaneous  equa- 
tions with  complex  coefficients.   It  follows  from  Equations  (2.1)  and 

(2.2)  that  three  independent  equations  are  sufficient  to  describe  the 
four-electrode  system.   Using  vector  and  matrix  notation  with  Vp(f ) 
vn(f)5  lQ1(f),  Ipi(f),  and  Ip2(f)  as  the  working  variables,  the  system 
equations  may  be  written 
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where  the  m  coefficients  are,  by  definition,  parameters  of  the  electrode- 
tissue  system  and  have  the  dimensions  of  complex  impedance.   This 
particular  matrix  form  was  chosen  to  facilitate  the  inclusion  of  inter- 
face effects  and  to  simplify  the  analysis  of  errors  given  in  later 
chapters.   The  minus  signs  are  included  so  that  the  real  part  of  each 
parameter  will  be  positive  in  a  passive  system.   For  example,  in  a 
purely  resistive  system,  all  the  parameters  would  be  real  and  positive. 

The  m  parameters  of  a  given  system  may  be  determined  experi- 
mentally or  estimated,  a  priori,  from  known  properties  of  the  tissue 
and  electrodes.   In  either  case,  one  refers  to  the  basic  definitions 
listed  in  Table  2.1,  which  may  be  obtained  by  inspection  from  Equation 
(2.3). 


TABLE  2.1 
DEFINITIONS  OF  THE  m  PARAMETERS 


Constraint   Ipi  =  Ip2  =  0   IQ1  =  I?2  =  0   IQ1  =  Ipi  =  0 


Excitation 

Q1-Q2 

VP 
mll  = 

V 

Pl- 

Q2 

P2- 

Q2 

m12  E 
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Ip2 

The  immediate  objective  of  most  impedance  measurements  is  to 
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identify  the  transfer  impedance,  m^.   This  parameter  is  estimated  by 
determining  the  ratio  of  an  arbitrary  measure  of  Vp  (the  nominal  voltage 
response)  to  an  arbitrary  measure  of  IQ1  (the  nominal  exciting  current) 
under  the  assumption  that  IPl  and  Ip2  are  zero.   As  shown  in  later 
chapters,  the  other  m  parameters  may  play  a  significant  role  in 
determining  the  accuracy  of  the  result. 

In  a  system  composed  of  bilateral  elements,  the  Reciprocity 
Theorem  from  network  theory  shows  that  only  six  of  the  nine  m  parameters 
are  independent.   Three  of  them  may  be  eliminated  by  means  of  the 
following  relations  (Appendix  1). 

m22  ~   m23  +  mn  (2  4) 

m32  =  ml2  +  m33  "  mll  (2.5) 

m2i  =  m31  +  m23  +  m±1  (2_6) 

The  m  parameters  are  determined  by  the  impedance  characteristics  of  the 
tissue  and  the  electrode-tissue  interfaces  and  by  the  system  geometry. 
These  factors  are  considered  in  the  remainder  of  the  chapter. 

Analysis  of  the  m  Parameters 

The  significance  of  the  m  parameters  may  be  determined  from  a 
study  of  the  general  field  relationships  that  exist  within  the 
electrode-tissue  system.   Consider  a  metal  electrode  in  an  electrolytic 
medium  in  the  presence  of  an  arbitrary  current  density  field  as  shown 
in  Figure  2.2.   Let  s  denote  an  imaginary  surface  which  surrounds,  and 
is  arbitrarily  close  to,  the  electrode  such  that  it  includes  essen- 
tially only  that  part  of  the  tissue  which  is  involved  in  interfacial 
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Figure  2.2.   Electrode-Tissue  Interface  in  an 
Arbitrary  Current  Field 

phenomena.   Then,  the  potential  of  the  electrode  with  respect  to  an 
arbitrary  point  of  reference  in  the  tissue  mass  may  be  written 

V  =  U  -  E  (2.7) 

where  E  is  the  average  potential  on  s  with  respect  to  the  reference 
point  and  U  is  the  average  potential  difference  across  the  interface. 
It  is  apparent  that  E  is  a  function  of  the  current  density  field  and 
the  tissue  impedance,  while  U  is  a  function  of  the  normal  current 
density  on  s  and  the  interface  impedance. 

At  any  position  on  s,  the  potential  difference  across  the  inter- 
face may,  in  principle,  be  written  as  the  product  of  the  complex 
impedance  per  unit  area,  x(s)>  and  the  normal  component  of  the  current 
density,  J"n(s),  at  that  point.   Then  the  average  potential  difference 
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is  given  by 

"sl|sXJnds  (2.8) 

where  the  integral  is  taken  over  all  of  s,  A  is  the  total  surface  area 
of  s  and  Jn  is  taken  as  positive  out  of  s. 

A  more  convenient  form  of  Equation  (2.8)  is  obtained  if  the 
following  substitutions  are  made 

X  =  X  +  n(s)  (2.9) 

Jn  =  Jn  +  v(s)  (2.10) 

where  x  and  Jn  are  the  average  values  of  interface  impedance  and  normal 
current  density  on  s,  and  n(s)  and  v(s)  are  deviations  from  the  averages, 
respectively.   Note  that  by  definition 


i 


's 
and 


n(s)ds=0  (2.11) 


rv  \s)   as  =  o  (2.12) 


v  (s)  ds  =  0 


If  I  designates  the  total  electrode  current  in  the  external  circuit, 
then  from  Equations  (2.10)  and  (2.12) 


Jn  =  t  (2.13) 


A 
With  Equations  (2.9)  through  (2.13),  Equation  (2.8)  may  be 


may  be  rewritten 


U  =  JL    I  +  l  I  nv  ds  (2.i4) 

A      A  Vo 
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The  first  term  on  the  right  in  Equation  (2.14)  expresses  the 
average  potential  difference  across  the  interface  due  to  current  flowing 
through  the  electrode  and  the  external  circuit.   The  second  term 
expresses  the  average  potential  difference  across  the  interface  due  to 
non-uniformity  in  the  interface  impedance.   For  example,  if  the  inter- 
face impedance  were  truly  uniform,  the  integral  would  be  zero  because 
n(s)  would  be  zero  over  all  of  s.   In  most  cases,  the  major  determining 
factor  in  n(s)  is  the  inherent  inhomogeneity  of  the  electrode  surface 
(Flasterstein,  1966b). 

The  relative  importance  of  the  terms  in  Equation  (2.14)  depends 
on  how  the  electrode  is  used.   For  the  current  electrodes  of  a  four- 
electrode  system,  the  second  term  is  undoubtedly  negligible  compared 
to  the  first.   However,  for  the  potential  electrodes,  where  I  is 
nominally  zero,  the  second  term  may  be  significant  (Schwan,  1963).   In 
cases  where  there  is  a  choice  of  position,  the  potential  electrodes 
should  be  placed  where  the  current  density  intercepted  by  the  electrodes 
will  be  low  (in  a  region  of  low  potential  gradient)  and  the  voltage 
between  the  electrodes  is  high.   In  this  way,  v  is  kept  to  a  minimum 
and  the  inhomogeneity  term  in  Equation  (2.14)  is  more  likely  to  be 
negligible  in  comparison  to  the  response  voltage.   For  example,  in  the 
parallel  wire  configuration  shown  in  Figure  2.3,  the  potential  electrodes 
are  placed  outside  rather  than  between  the  current  electrodes,  where 
the  same  potential  difference  may  be  recorded  and  the  current  density 
is  smaller. 

The  Effective  Interface  Impedances 

In  a  system  of  four  electrodes,  Equation  (2.7)  holds  for  each 
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Figure  2.3.   Parallel  Wire  Four-Electrode  System 

electrode.   Using  the  appropriate  subscripts  and  summing  the  electrode 
potentials  according  to  Kirchhoff's  Voltage  Law  and  Figure  2.1,  the 
vector  equivalent  of  Equation  (2.7)  is  found  to  be 


(upi  -  uP2) 


(UQ1  "  UQ2)    +     En  (2.15) 


(uQ1  -  uP1) 


where  Ep,  Eq,  and  E1   are  the  average  potential  differences  between  the 
s  surfaces  of  electrode  pairs  P1-P2,  Q1-Q2,  and  Ql-Pl,  respectively. 

The  present  investigation  is  primarily  concerned  with  the  effects 
of  the  first  term  on  the  right  in  Equation  (2.14).   To  simplify  the 
presentation,  the  term  due  to  inhomogeneity  will  be  neglected  from  this 
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point  on.   With  this  omission  and  adding  the  appropriate  subscripts 
in  Equation  (2.14),  the  average  interfacial  potential  differences  of 

» 

the  system  in  Figure  2.1  become 
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where,  for  brevity,  the  ratio  —  is  replaced  by  Z  for  each  interface. 
Henceforth,  Zq^,  Zpi,  Zp2,  and  ZQ2  are  referred  to  as  the  effective 
interface  impedances . 

Using  Equations  (2.16)  through  (2.19),  Equation  (2.15)  may  be 
written 
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(  2 .  20) 


where  I_j  has  been  eliminated  by  means  of  Equation  (2.2).   Through  the 
effective  interface  impedances,  Equation  (2.20)  gives  explicit  repre- 
sentation to  the  interface  potential  differences  and  is  shown  diagra- 
matically  in  Figure  2.4. 
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Figure  2.4.   Electrode-Tissue  System  Showing  the 
Effective  Interface  Impedances 


The  h  Coefficients 


In  Figure  2.4,  Ql,  PI,  P2  and  Q2  may  be  viewed  as  the  terminals 
of  a  hypothetical  electrode-tissue  system  whose  terminal  voltages  are 
equal  to  the  average  surface  voltages  of  the  real  system.   These 
voltages  may  be  related  to  the  electrode  currents  by  means  of  a  matrix 
equation  as  follows : 
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—          — 

hi 

XP1 

Ip2 

hll   h12   -h13 


h21   h22    h23      X^  (2.21) 


h31  "h32   -h33 


where  the  definitions  of  the  h  coefficients  follow  from  those  given  for 
the  m  parameters  in  Table  2.1  with  the  appropriate  substitutions. 

From  Equations  (2.20)  and  (2.21),  it  is  evident  that  the 
impedance  parameter  matrix  in  Equation  (2.3)  is  given  by 


mll   mi2  "m13 


m21    m22    m93 


m31  ~m32  "m33 


ZP1   "ZP2 


zQ1+zQ2   zQ2   zQ2 


JQ1 


Note  that  for  all  frequencies 


-zpi   0 


h21    h22 


h31   "h32 


hll    h12    "hi3 


J23 


-h 


33 


(2.22) 


mll  "  hll 


(2.23) 


and 


m33  "  h33 


(i.ih-) 


At  low  frequencies  (e.g.,  below  ICO  Hz),  the  elements  in  the  interface 
impedance  matrix  are  much  larger  than  those  in  the  h  matrix  in  all 
except  the  11  and  33  positions.   At  high  frequencies  (e.g.,  above 
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100  kHz),  the  elements  in  the  interface  impedance  matrix  tend  to  become 
small  compared  to  those  in  the  h  matrix.   However,  in  many  cases,  the 
interface  impedances  may  not  become  negligible  before  the  errors  of 
measurement  exceed  permissible  limits.   This  possibility  is  embodied 
in  the  results  of  succeeding  chapters. 

It  may  appear  that  interface  effects  are  expressed  entirely  by 
the  interface  impedance  matrix.   This  statement  is  only  approximately 
true.   The  geometry  of  the  current  field  is  determined  mainly  by  the 
electrode  configuration.   However,  it  is  also  influenced  in  varying 
degree  by  the  frequency-dependent  boundary  conditions  imposed  by  the 
impedance  of  the  electrode-tissue  interfaces.   This  factor  is  not  likely 
to  be  important  except  with  respect  to  its  direct  effect  on  h,  -,  .   For 
the  sake  of  completeness  and  to  demonstrate  the  complexity  of  the 
four-electrode  field  problem,  a  brief  discussion  of  possible  boundary 
effects  is  in  order.   The  examples  that  follow  are  of  a  hypothetical 
nature  and  have  not  been  studied  in  experimental  or  analytical  detail. 
The  existence  of  boundary  condition  effects  is  evidenced  by  certain 
anomalous  behavior  of  electrodes  in  saline  which  is  difficult  to  explain 
in  terms  of  other  sources  of  error. 

Examples  of  boundary  condition  effects 

To  indicate  how  interfacial  boundary  conditions  may  affect  h,,, 
consider  again  the  electrode  configuration  of  Figure  2.3.   The  exposed 
tips  of  the  potential  probes  are  located  near  the  surface  of  the  tissue 
mass  where  the  current  density  field  is  fairly  planar.   Assume  that  a 
sinusoidal  current  of  constant  amplitude,  I_.  ,  excites  the  system.   At 
high  frequencies,  where  the  interface  impedance  is  small,  the  current 
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density  distribution  along  Ql  and  Q2  is  determined  primarily  by  the 
configuration  of  the  metalic  boundaries  and  is  indicated  qualitatively 
in  Figure  2.5  (high).   The  rise  toward  the  electrode  tips  is  due  to 
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Figure  2.5.   Electrode  Current  Density  in  the  System  of 
Figure  2.3 


fringing  into  the  tissue  mass.   At  low  frequencies,  the  interface 
impedances  of  Ql  and  Q2  are  large  compared  to  the  tissue  impedance, 
and  they  tend  to  smooth  the  current  density  distribution,  as  shown  in 
the  same  figure  (low).   Since  the  total  current  is  constant,  the  areas 
between  the  low  and  high  frequency  graphs  are  equal,  and  there  is  more 
current  at  low  frequencies  in  the  field  of  PI  and  P2  than  at  high 
frequencies.   The  importance  of  this  effect  should  depend  on  the  rela- 
tive dimensions  of  the  current  electrodes  and  become  smaller  as  the 
electrode  length  is  increased. 
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A  useful  expression  for  h-,-,  is 


(2.25) 


where  T  is  a  geometrical  factor  and  $  is  the  specific  impedance  of  the 
tissue.   Since  r  reflects  the  current  density  distribution,  it  varies 
with  frequency  and  for  the  configuration  of  Figure  2.3  takes  the 
general  form  shown  in  Figure  2.6.   At  low  frequencies,  where  the  effect 
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Figure  2.6.   Frequency  Dispersion  of  the  Geometrical 
Factor  for  the  System  of  Figure  2.3 


of  interface  impedance  is  much  greater  than  that  of  the  tissue  impedance, 
the  graph  is  flat.   When  the  impedances  are  comparable,  r  decreases 
with  frequency  and  assymptotically  approaches  a  value  corresponding  to 
zero  interface  impedance. 

Since  the  electrode-tissue  system  is  assumed  to  be  linear  and 
passive,  it  should  obey  the  phase-magnitude  relationships  of  network 
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theory  (Bode,  1945).   Hence,  the  magnitude  of  hj,  cannot  vary  with 
frequency  without  a  concomitant  change  in  its  phase  angle  as  well. 
This  implies  that  there  can  be  no  frequency  dependence  of  |r|  unless  r 
is  a  complex  number.   Indeed,  careful  consideration  of  the  phenomenon 
illustrated  by  Figure  2.5  indicates  that  the  phase  angle  of  the  current 
density  may  vary  along  the  electrode  length  even  though  the  phase  of 
the  total  current  is  fixed.   The  phase  of  the  current  toward  the  tips 
of  Ql  and  Q2  would  tend  to  be  positive,  while  that  toward  the  tips  of 
PI  and  P2  would  tend  to  be  negative,  with  respect  to  the  phase  of  the 
total  current.   Hence,  in  agreement  with  network  theory,  a  decrease 
in  Ih-Ql  (or  |r|)  with  frequency  would  be  accompanied  by  a  negative 
phase  angle  in  h^   (or  r).   In  view  of  these  considerations,  the 
ordinates  in  Figures  2.5  and  2.6  should  be  interpreted  as  the  magnitude 
of  the  quantities  indicated.   Further,  the  areas  between  the  high  and 
low  curves  in  Figure  2.5  are  not  equal  since  phase  angles  must  be  taken 
into  account. 

Experimentally,  r  is  often  estimated  by  testing  the  electrode 
system  in  saline  of  known  specific  resistance.   The  frequency  dependence 
of  r  may  be  expected  to  vary  with  the  ratio  of  interface  impedance  to 
saline  impedance.   For  example,  as  saline  concentration  is  increased, 
this  ratio  increases,  and  the  dispersion  in  r  may  be  expected  to  occur 
at  higher  frequencies.   If  the  frequency  dependence  of  T   is  significant 
in  a  given  system,  it  becomes  necessary  to  use  a  saline  concentration 
which  produces  a  similar  ratio  of  interface  impedance  to  specific 
impedance  as  occurs  in  the  tissue  to  be  measured.   This  problem  is 
complicated  by  the  masking  effect  of  cells  in  the  tissue  on  the  inter- 
face impedance.   Hence,  rather  than  use  saline  of  approximately  the 
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same  specific  impedance  as  the  tissue,  the  saline  impedance  should  be 
like  that  of  the  extracellular  fluid  which  is  much  lower. 

Another  possible  boundary  condition  effect  on  hlx  is  due  to 
non-uniformity  of  the  interface  impedance  along  the  electrode  surface. 
Non-uniformity  at  the  potential  electrode  interfaces  has  already  been 
discussed  in  connection  with  Equation  (2.14).   When  it  occurs  at  the 
current  electrode  interfaces,  it  may  alter  the  current  distribution  in 
varying  and  indeterminate  ways.   For  example,  the  low  frequency  plot 
in  Figure  2.5  might  actually  have  an  irregular  shape  which  varies 
spontaneously  with  time.   Further,  if  the  phase  angle  of  the  interface 
impedance  varies  over  the  electrode  surface,  r  may  become  a  complex 
number  introducing  an  erroneous  phase  angle  in  h. 

The  guard  electrode  principle.— In  the  system  of  Figure  2.3, 
the  constancy  and  definition  of  the  current  field  in  the  potential 
electrode  region  may  be  improved  by  using  two  electrically  separate 
sections  for  one  of  the  current  electrodes.   The  two  sections  should 
appear  continuous  except  for  a  minute  break  occurring  at  an  arbitrary 
point  below  the  potential  electrode  tips.   The  controlled  exciting 
current  is  caused  to  flow  through  the  upper  section,  while  the  potential 
of  the  lower  section  is  made  to  follow  the  potential  developed  by  the 
upper  section  by  means  of  external  circuitry.   From  the  field  viewpoint, 
the  two  sections  appear  as  a  single  electrode,  but  with  the  current  in 
the  region  of  interest  under  direct  control.   The  guard  electrode 
principle,  as  it  is  commonly  called,  may  be  applied  in  many  different 
configurations  (see,  for  example,  Van  Harreveld,  et  al.,  1963;  Graham, 
1965b). 
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The  Equivalent  Network 

The  significance  of  interface  impedance  and  geometrical  factors 
is  enhanced  if  the  four  electrode  system  is  viewed  in  terms  of  an 
equivalent  network.   Consider  the  hypothetical  system  described  by 
Equation  (2.21).   As  indicated  earlier  for  the  m  parameters,  in  a 
system  obeying  reciprocity,  only  six  of  the  nine  h  coefficients  are 
independent.   Such  a  system  may  be  represented  by  a  network  with  six 
independent  impedance  elements.   A  convenient  form  for  this  network  is 
shown  in  Figure  2.7  where  W  denotes  complex  impedance.   This 
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Figure  2.7.   Equivalent  Network  for  the  Hypothetical 
System 


configuration  was  chosen  because  it  places  in  evidence  certain  implicit 
properties  of  the  electrode-tissue  system.   For  example,  Wpi  and  Wp2 
may  be  interpreted  as  electrode  impedances  which  characterize  the  effect 
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of  drawing  current  through  the  potential  probes.   The  remaining  four 
elements  split  the  exciting  current  into  two  parts:   one  passing 
through  W-j_,  Wp,  and  W2,  which  generates  the  response  voltage  between 
points  1  and  2,  and  one  passing  through  WQ,  which,  in  a  sense,  is 
wasted.   In  many  geometries,  the  effect  of  WQ  is  negligible,  and  the 
network  simplifies  to  a  double  tee  form. 

The  matrix  equation  governing  the  network  of  Figure  2.7  is 


YW- 


Y(W1  +  W2  +  Wp) 


YW, 


Wl 
(YWp  +  p-Wp  +  Wpi) 


Y(W2  +  Wp) 

w    w 

-^2  +  D^P  +  WP1> 


W2 
-(Wp2  +  fV?) 


YW_ 


Wl 
iT*2 


LQl 


LPl 


LP2 
(2.26) 


where 


Y 


wx  +  w2  +  wp  +  WQ 


(2.27) 


and 


D  =   Wx  +  W2  +  Wp  +  WQ 


wx  +  w2  +  wp 

1  -  Y 


(2.28) 


Setting  the  matrices  of  Equations  (2.21)  and  (2.26)  equal,  the  following 


expressions  are  found  for  the  network  elements 


W    -  v, 

WP1  "  n12  "  "11 


W 


h23 


'33 


P2    13    llh 


(2.29) 
(2.30) 


31 


wp  =  7  hn 


(2.31) 
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W1=Ih31  (2.32) 


W2=7h23  (2'33) 


..    1  h23h31  ,_  _,  , 

Wn  =  Z~ T (2.34) 

Q   Y   h33 


and 


7=1  +  f4^  <2-35> 

Y       h23h31 


In  general,  the  network  elements  are  frequency  dependent  through  the 
tissue  impedance  and  the  interfacial  boundary  conditions. 

As  noted  earlier,  the  immediate  objective  of  four-electrode 
measurements  is  to  identify  h-,-,  .   In  terms  of  the  equivalent  network 


h11  =  YWp  (2.36) 


where  Y,  as  defined  in  Equation  (2.27),  is  the  fraction  of  I-,  which 
flows  through  Wp  when  Ipj_  and  Ip«  are  zero.   As  shown  in  the  examples 
at  the  end  of  the  chapter,  the  value  of  Y  depends  on  the  electrode 
configuration,  but  it  may  not  be  greater  than  unity. 

To  complete  the  equivalent  network,  the  effective  interface 
impedances  are  added  as  indicated  in  Figure  2.8.   At  low  frequencies 
(e.g.,  below  100  Hz),  Zpi  and  Zp2  are  much  larger  than  Wp  and  Wp2, 
respectively.   At  high  frequencies,  the  interface  impedances  should 
become  negligible  in  the  limit.   However,  as  indicated  for  the  matrix 
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Figure  2.8.   Complete  Equivalent  Network  of  the  Electrode- 
Tissue  System 

parameters  of  Equation  (2.22),  this  condition  may  not  occur  before  the 
errors  of  measurement  exceed  permissible  limits . 

In  the  topics  treated  later  in  this  work,  the  matrix  parameters 
are  used  almost  exclusively  because  they  are  more  basic,  easier  to 
manipulate,  and  lead  to  more  concise  results.   However,  some  phenomena 
are  more  readily  understood  when  viewed  in  terms  of  the  equivalent 
network. 

Examples  of  Impedance  Parameter  Matrices 


The  major  results  of  this  chapter  are  illustrated  in  the 
examples  that  follow.   To  simplify  the  computation,  the  configurations 
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are  idealized,  and  the  effect  of  interfacial  boundary  conditions  on  the 
h  coefficients  is  neglected.   It  is  assumed  that  all  electrodes  have 
the  same  average  interface  impedance  per  unit  area  ( x ) ,  and  the 
impedance  medium  is  homogeneous  with  respect  to  the  dimensions  of  the 
main  field  region.   The  specific  impedance  of  the  medium  is  denoted  by 
*. 

The  h  coefficients  are  found  from  the  basic  definitions  in  Table 
2.1  with  appropriate  modifications  (compare  Equations  [2.3]  and  [2.21]). 
The  effective  interface  impedances  follow  from  the  definitions  in 
Equations  (2.16)  through  (2.19).   Finally,  the  network  elements  may 
be  calculated  from  Equations  (2.29)  through  (2.35),  if  desired.   For 
convenience,  the  results  are  tabulated  as  matrices  following  the  form 
of  the  right  hand  terms  in  Equation  (2.22). 

In  each  example,  numerical  values  are  illustrated  for  the  matrix 
elements  at  low  and  high  frequencies.   The  computations  make  use  of 
the  approximate  values  of  £  and  $  shown  in  Table  2.2,  which  are  based 
on  observations  by  the  author  using  stainless  steel  electrodes  in  0.9 
per  cent  (isotonic)  saline.   The  factor  e~J6  is  defined  in  Chapter  1. 

TABLE  2.2 

APPROXIMATE  VALUES  OF  £  AND  *  FOR  STAINLESS  STEEL  ELECTRODES 

IN  ISOTONIC  SALINE 

Frequency x  ohm-cm2 A   ohm- cm 


3  Hz  1,000  e-372°  60  ej0° 

100  kHz  0.5  e"372°  6Q  ej0° 


Here,  it  indicates  a  phase  angle  of  minus  72  degrees.   The  complex 
number  equivalent  is 
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e-j72°  _  cos  72o  _  jsin  72o 
0.31   -  jO.95 

The  examples  are  calculated  for  the  system  shown  in  Figure  2.9,  which 
consists  of  four  hemispherical  electrodes  set  into  the  surface  of  a 
semi-infinite  impedance  medium.   Electrode  radii  and  separation  are 
represented  by  r  and  d  with  appropriate  subscripts,  respectively. 
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Figure  2.9.   Four-Electrode  System  with  Hemispherical 
Electrodes 
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Example  1 

The  approximate  element  expressions  in  Table  2.3  are  determined 
for  the  following  geometrical  conditions 


drVd2Ed 


and 


r       =  r       =  r       =r       =  r,<<d 
Ql  PI  P2  Q2    "   r      Q 


The  numerical  results  correspond  to  d=0.2  cm  and  r=0.01 


cm. 


Frequency 


TABLE  2.3 
MATRIX  ELEMENTS  FOR  EXAMPLE  1 
Interface  Impedance 


h  Coefficients 


All 

0 

1/2 

-1/2 

l/2d 

l/2r 

-l/2r 

X/?rr2 

1 

1/2 

1/2 

Ip/TT 

1/r 

l/2r 

l/2r 

1/2 

-1/2 

0 

l/2r 

-l/2r 

-l/6d 

3  Hz 

0 

1.6 

-1.6 

47.7 

955 

-955 

106e^72 

3.2 

1.6 

1.6 

1910 

955 

955 

1.6 

-1.6 

0 

955 

-955 

-15.9 

-J 

100  kHz 

0 

800 

-800 

47.7 

955 

-955 

e-J72 

1600 

800 

800 

1910 

955 

955 

800 

-800 

0 

955 

-955 

-15.9 

In  this  example,  Y  is  approximately  unity  at  all  frequencies, 


Example  2 


The  geometrical  conditions  for  this  example  are 


d-^  =  dp  =  d 


rQl  =  rPl  =  rP2  ~   r<<d 


d2>>d 


rQ2  =  2d 
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The  approximate  element  expressions  are  given  in  Table  2.4,  where  the 
numerical  results  correspond,  as  before,  to  d=0.2  cm  and  r=0.01  cm. 


Frequency- 


All 


X/2TT 


3  Hz 


106e^'72 


TABLE  2.4 
MATRIX  ELEMENTS  FOR  EXAMPLE  2 
Interface  Impedance 


h  Coefficients 


0 
1/r2 


1/r' 


1/U&4 


l/rz   -1/r' 


l/4d: 

0 


ip/2-rr 


0     1.6     -1.6 
1.6    0.001    0.001 
1.6    -1.6       0 


In  this  example,  Y  is  2/3  at  all  frequencies 


l/2d  1/r  -1/r 
1/r  3/2d  1/d 
1/r    -1/r   -l/2d 


23.8    955    -955 
955     71.4   47.7 
955    -955    -23.8 


100  kHz 

0 

800 

-800 

23.8 

955 

-955 

e-j72 

800 

0.5 

0.5 

955 

71.4 

47.7 

800 

-800 

0 

955 

-955 

-23.8 
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Concluding  Remarks 

The  analysis  in  this  and  succeeding  chapters  uses  methods  of 
linear  mathematics  which  assume  that  the  electrode-tissue  system  is 
linear  throughout.   In  practice,  while  the  exciting  and  observed 
variables  may  be  linearly  related,  the  voltages  and  currents  at  the 
current  electrode  interfaces  may  not  necessarily  be  so,  depending  on 
the  current  densities  and  frequencies  encountered.   However,  this 
discrepancy  detracts  little  from  the  value  of  the  linear  approach  in 
understanding  and  estimating  system  performance.   A  linear  analysis 
also  serves  as  a  basis  for  the  inclusion  of  nonlinear  effects. 


CHAPTER  3 


ERRORS  IN  THE  IDENTIFICATION  OF  TRANSFER  IMPEDANCE  BY 
THE  CLASSICAL  METHOD 


Introduction 

The  classical  mode  of  operation  in  four-electrode  measurements 
is  to  ground  one  current  electrode  and  record  differentially  from  the 
potential  electrodes.   The  purpose  of  this  chapter  is  to  define  and 
analyze  the  errors  of  this  method  and  provide  a  basis  of  comparison 
for  the  more  advanced  technique  in  the  next  chapter. 

Errors  in  the  identification  of  m^   include  four  major  types. 
Briefly,  they  are 

1.  loading  error 

2.  common  mode  rejection  error 

3 .  current  control  error 

4.  noise 

Each  of  these  errors  is  discussed  and  the  first  three  are  analyzed  in 
detail  in  this  chapter. 

A  basic  block  diagram  for  impedance  measurements  is  shown  in 
Figure  3.1.   The  signal-processing  and  display  systems  are  irrelevent 
to  the  present  study  and  have  been  omitted.   The  connecting  system 
comprises  all  the  conductors,  connectors,  and  supporting  materials  that 
link  the  electrode-tissue  system  to  the  instrumentation.   The  electrode- 
tissue  system  includes  the  tissue  and  only  those  parts  of  the  metal 
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Figure  3.1.   Basic  Block  Diagram  for  Impedance  Measurements 

electrodes  which  are  in  contact  with  the  tissue.   The  exciting  system 
usually  includes  a  signal  generator  and  means  for  controlling  or  sensing 
the  applied  current.   The  amplifying  system  amplifies  the  voltage  response, 

Loading  Error 

The  term  loading  admittance  will  be  used  to  refer  to  the 
extraneous  admittances  in  an  impedance  experiment  that  permit 
undesirable  currents  to  flow  in  the  electrode-tissue  system.   For 
clarity,  it  may  be  viewed  as  consisting  of  four  arbitrarily  defined 
classes  of  components  : 
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1.  transverse  admittance 

2.  cable  admittance 

3.  amplifier  input  admittance 
•».  ground  admittance 

Transverse  admittance . —Transverse  admittance  occurs  in  the 
connecting  system  and  denotes  parasitic  current  paths  between  any  or 
all  of  the  signal  conductors.   It  arises  from  fixed  and  indeterminate 
properties  of  the  insulators  which  support  the  electrodes  and  various 
connectors.   Indeterminate  properties  are  primarily  due  to  bulk  and 
surface  contamination  of  the  insulators  by  such  media  as  blood, 
saline,  water,  and  solder  flux.   Since  it  has  a  somewhat  amorphous 
origin,  transverse  admittance  is  appropriately  described  by  means  of 
a  3  by  3  matrix,  which  may  be  converted  to  an  equivalent  network. 
If  the  transverse  admittance  obeys  reciprocity,  the  matrix  and  the 
network  each  contain  six  independent  elements. 

Cable  and  amplifier  input  admittance  .—These  classes  of 
admittance  are  commonly  specified  in  terms  of  finite  elements  between 
ground  and  the  signal  conductors.   Elements  may  also  occur  between 
the  conductors  directly. 

Ground  admittance  .—Ground  admittance  encompasses  all  stray 
current  paths  occurring  between  the  mass  of  the  tissue  and  ground. 
In  a  practical  sense,  it  may  arise  unintentionally  by  contact  with 
grounded  objects  such  as  equipment,  people,  and  wet  surfaces.  For 
present  purposes,  the  manner  in  which  the  ground  admittance  is 
distributed  need  not  be  known.   It  is  hereafter  represented  as  a 
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finite  element,  connected  between  ground  and  an  arbitrary  point  in 
the  tissue. 

The  first  three  classes  of  loading  admittance  may  be  viewed  as 
a  four-terminal  network  in  parallel  with  the  electrode-tissue  system. 
This  network  is  indicated  by  Ql' ,  PI',  P2',  and  Q2'  in  Figure  3.2.   In 
the  same  figure,  ground  admittance  is  represented  by  Y  . 
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Figure  3.2.   Schematic  Diagram  of  the  Electrode-Tissue  and 
Loading  Admittance  Systems 


The  desired  parameter  of  the  electrode-tissue  system,  m ..,  is 
estimated  as  the  ratio  of  the  voltage  response,  V  ,  to  the  source 
current,  I_.   According  to  the  definition  of  m   in  Table  2.1,  the 
result  is  based  on  the  assumption  that  Ip.  and  Ip9  are  zero,  and  that 


I„,  Iq,  and  -Iqo  ars  equal.   This  assumption  is  valid  only  if  the 


68 


loading  currents  I   ,  Ipi,  Ip2  and  I  are  zero.   Each  non-zero  loading 
current  causes  an  error  in  the  estimated  value  of  m,  -,  ,  which,  by- 
definition,  is  the  loading  error  of  that  current. 

The  errors  of  admittance  classes  1,  2,  and  3  are  analyzed 
together  in  the  next  section.   For  clarity,  the  error  due  to  ground 
admittance  is  treated  separately. 

Errors  Due  to  Admittance  Classes  1,  2,  and  3 

It  is  convenient  to  represent  the  first  three  classes  of 
loading  admittance  by  the  equivalent  network  defined  in  Figure  3.3 
which  contains  six  independent  elements.   The  current  vector  of  this 
network  is 


XQ1 


lk 


lk 


■3  YQ    (Y1+V 


Y     y    -(Y  +Y  ) 

p    4   v  r  V 


(Yp+Y2+Y3)    Y2   -(Y2+Y3) 


V 

VQ 

vl 

L     —J 

(3.1) 


With  reference  to  Figure  3.2,  the  current  vector  of  the 
electrode-tissue  system  is  given  by 


Ql 


lk 


Pi 


r 
pi 


(3.2) 


P2 


JP2 
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Figure  3.3.   Equivalent  Network  for  Loading  Admittance 
Classes  1,  2,  and  3 


and  from  Equation  (2.3),  the  voltage  vector  it 


mll    m12    "m13 


m21    m22     m23 


m    -  m      -m 
31     32      33 


V 

XP1 

I 

P2 

(3.3) 


It  follows  from  Equations  (3.2)  and  (3.3)  that 
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Vp   =  mnIQ  -  ra1±VQ1  -  m^I^  +  m^I^ 


(3.4) 


where  the  terms  involving  I*    Ipi  and  Ip2  are,  by  definition,  loading 
errors  in  Vp  which  ideally  should  be  zero.   Dividing  these  terms  by 
Vp  and  substituting  for  the  currents  from  Equation  (3.1),  the  relative 
loading  errors  in  Vp  may  be  written 


■mllIQl 


Ql 


'11 


VQ         Vl     Vl 

v7  YQ  +  (1  +  v^)Y3  +  TT  Yi 
P         vp      P 


(3.5) 


■m12IPl 


PI 


12 


V     V 

Vp    Vp;  4 


Vl 

— -  Y  +  Y 

Vp  ll  P 


(3.6) 


and 


m13Ip2 


P2 


=  m 


13 


(^ 


1)Y   -  (1  +  rr-)Y      -   Y^ 
2        vp   3    P 


(3.7) 


The  operating  condition  of  primary  interest  in  four-electrode 
measurements  occurs  when  a   ,  api  and  ap2  are  quite  small.   Under  this 
condition,  it  is  possible  to  approximate  the  voltage  ratios  in  the 
bracketed  terms  of  Equations  (3.5),  (3.6),  and  (3.7)  by  expressions 
involving  only  the  m  parameters.   The  procedure  is  justified  by  the 
fact  that  when  the  a   are  small,  the  errors  in  the  approximating 
expressions  will  also  be  small  as  explained  in  the  next  paragraph. 

According  to  Equation  (3.3),  V  ,  V   and  W±   each  include  a 
contribution  from  I   and  I   .   If  one  considers  the  geometrical  and 
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interfacial  properties  of  four-electrode  systems,  it  is  not  difficult 
to  see  that  V±   and  especially  Vq  are,  in  most  cases,  less  sensitive  than 
Vp  with  respect  to  Ipi  and  Ipr   In  the  analysis  that  follows,  the 
error  in  Vp  due  to  these  currents  is  required  to  be  less  than  0.05V 
in  magnitude.   It  is  reasonable  to  assume,  therefore,  that,  under  this 
condition,  the  corresponding  errors  in  V  and  V1  will  similarly  be 
less  than  0.05  V  and  0.05  V.^  respectively. 

To  determine  the  voltage  ratios,  the  errors  in  Vp,  V„  and  V, 
due  to  Ipi  and  Ip2  will  be  neglected  in  accord  with  the  foregoing 
discussion.   With  this  understanding,  Equation  (3.3)  yields 


V„  a  mnnI 


P  "  VQ1  0.8) 


VQ  *  'VqI  =  (m31  +  m23  +  "ll^Qi  (3.9) 


and 


Vl  =  ffl31IQl  (3.10) 

where  the  dependent  parameter,  m^,  has  been  replaced  by  means  of 
Equation  (2.6).   It  follows  that  the  voltage  ratios  are 


Vp  ~  m  (3.11) 


and 
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VQ   m31  +  m23  +  mn 

_  = j£ (3.12) 

P         mll 


Using  Equations  (3.11)  and  (3.12),  the  expressions  for  the 
relative  errors  become 


V  ~~   "[(mil  +  m31  +  m23)YQ  +  (mil  +  m31)Y3  +  m31Yl]  (3"13) 


~m12 
"PI  "  -^   [(mH  +  m23)Y4  "  m31Yl  +  milV  (3'1^) 


m13r 
aP2  *   ^m23Y2  "  (mil  +  m31)Y3  '  milV  (3'15> 


Magnitude  and  phase  errors 

Before  continuing  with  the  analysis  of  loading  error,  it  is 
appropriate  to  determine  how  the  error  is  manifested  in  the  experimental 
results.   The  experimental  estimation  of  the  desired  impedance 
parameter,  m   ,  is  defined  by 


VP 
mll  '=   IJ  (3-16) 


If  the  total  relative  loading  error  in  Vp  due  to  admittance  classes 
1,  2,  and  3  (i.e.,  the  sum  of  aQ1,  api  and  ap2)  is  designated  by  a 


73 


then,  from  Equation  (3.4) 


A.  mi  1 

m..  =  i±-  (3.17) 

H   1  -  a 
L 


It  follows  from  Equation  (3.17)  that  the  relative  magnitude 


.   A 

error  in  m, ,  is 


mi  n   _  m-i  t        t 
H=— LJil  =  -r—^ r-1  (3.18) 

lmnl  li  -  aLl 


and  the  phase  error  is 


/A 


$   E  /Sll  "      All   =   -       A1  -  aL)  (3.19) 


To  study  the  dependence  of  the  magnitude  and  phase  errors  on  ar  , 

L 

consider  the  phasor  diagram  shown  in  Figure  3.4.   If  a,  is  allowed  to 

Li 

rotate  in  the  complex  plane,  the  algebriac  signs  of  M  and  $  vary  as 
indicated  in  Table  3.1.   For  a  given  | ot  |,  the  relative  magnitude 

i_i 

error  becomes  maximum  when  the  angle  of  a  is  0  or  180  degrees ,  and 

the  phase  error  becomes  maximum  when  the  angle  of  aT  is  ±(90  -  sin~^la  I) 

degrees.   For  small  |oT | ,  the  maximum  errors  may  be  written 


M  *   ±|a J  (3.20) 


with  zero  phase  error,  and 


Figure  3.4.   Phasor  Relations  of  a   in  the  Complex  Plane 

Li 


TABLE  3.1 


SIGN  OF  MAGNITUDE  AND  PHASE  ERRORS  IN  £ 


11 


Quadrant  of  ctr  ~ 

u 


II 


III 


IV 


"More  precisely,  the  sign  of  M  changes  when  the  angle  of  a.  passes 
through  " 

±  [90  -   -  cos_1(l  -  |  |aL|2)]  degrees, 
which  for  |aL|  <0.05  is  within  2.1  degrees  of  the  quadrature  axis. 
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II 


1  -  a 


III 


IV 
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±57.3    |aj    degrees  (3.21) 


laLl 


j2 

with  approximately  relative  magnitude  error. 

2 

In  the  computations  that  follow,  |a, |  will  be  limited  to  less 
than  0.05,  corresponding  to  magnitude  and  phase  errors  within  ±5  percent 
and  ±2.9  degrees  of  zero,  respectively.   In  practice  (e.g.,  when 
measuring  very  low  frequency  dispersion  phenomena)  a  smaller  limit  for 
|ojj  may  be  necessary.   In  that  case,  it  is  a  simple  matter  to  modify 
all  of  the  results  by  the  appropriate  constant  factor. 

Low  frequency  errors 

If  reasonable  care  is  taken  in  the  construction  and  maintenance 
of  the  connecting  system,  low  frequency  errors  do  not  become  significant 
until  the  low  frequency  approximations  of  the  m  parameters  indicated 
in  connection  with  Equation  (2.22)  become  valid.   At  low  frequencies, 
the  Y  coefficients  will  normally  be  pure  conductances.   Making  the 
appropriate  substitutions  in  Equations  (3.13),  (3.14),  and  (3.15),  and 
denoting  Y  by  G  for  conductance,  the  low  frequency  relative  loading 
errors  become 


V  =  "  C(ZQ1  +  ZQ2)GQ  +  ZQ1G3  +  ZQ1G1]  (3"22) 

aPl  =  1T~  EZq2G4  "  ZQ1G1  +  hllGpl  0.23) 

ZP2 

aP2  =  n^[ZQ2G2  "  ZQ1G3  "  hllGP^  <3-2*> 
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Since  the  effective  interface  impedances  are  generally  much  larger 
than  hi;L,  a  q±   will  ordinarily  be  negligible  compared  to  ap,  and  ap2. 
For  the  same  reason,  the  terms  including  Gp  are  likely  to  be  negligible 
also. 

When  the  angle  of  h^  is  small,  the  phase  angle  of  a.  is  given 

approximately  by  the  sum  of  the  angles  of  two  interface  impedances 

plus  0  or  180  degrees,  depending  on  the  final  sign  of  the  result.   For 

example,  for  the  phase  angles  given  in  Table  2.2,  the  angle  of  a   may 

Li 

fall  in  the  first  or  third  quadrant.   In  this  case,  the  magnitude  and 
phase  errors  are  both  positive  or  both  negative  as  shown  in  Table  3.1, 
respectively. 

If  electrolytic  residue  contaminates  the  connecting  system, 
metal-electrolyte  interfaces  may  form  at  the  signal  conductor  surfaces 
and  act  in  series  with  the  residue  mass.   In  this  case,  the  Y 
coefficients  may  be  primarily  conductive,  or  they  may  include  large 
susceptive  components.   The  outcome  depends  on  the  relative  importance 
of  the  effective  impedances  of  the  residue  mass  and  the  interfaces 
formed.   Clearly,  unless  these  impedances  are  quite  high,  the  analysis 
of  small  errors  developed  here  is  not  valid. 

If  the  expressions  for  a    api  and  ap2  are  combined,  it  is 
evident  that  the  total  relative  error,  ct-,  comprises  six  terms,  each 
multiplied  by  a  different  G  coefficient.   Since  the  G  coefficients  are 
independent  finite  parameters  of  the  physical  system  (see  Figure  3.3), 
it  is  desirable  to  determine  what  constraints  they  should  satisfy  to 
assure  a  certain  level  of  experimental  accuracy.   The  simplest  approach 
is  to  place  a  suitable  upper  limit  on  the  magnitude  of  each  of  the  six 
terms  in  a^.   Since  these  terms  possess  different  algebraic  signs  and 
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phase  angles,  the  resulting  upper  limit  on  |a  |  can  only  be  estimated 
in  the  general  case.   The  terms  including  G  ,  Gp,  G   ,   and  G^  carry 
negative  signs,  while  those  of  G  and  G  carry  positive  signs. 
Further,  the  G3  and  G^  terms  tend  to  cancel  the  Gx  and  G2  terms;  and 
the  terms  including  GQ  and  Gp  are  usually  much  smaller  than  the  others. 
Accordingly,  the  upper  limit  for  j ctL [  may  be  determined  by  setting 
the  terms  of  G3  and  G^  (or  G±   and  G2)  at  their  upper  limits  and 
neglecting  the  remaining  terms.   If,  for  example,  the  magnitude  of 
each  of  the  six  terms  in  a]j  is  limited  to  less  than  0.025,  then  |a  I 
may  be  assumed  to  be  less  than  0.05. 

Combining  Equations  (3.22),  (3.23)  and  (3.24)  and  applying  the 
0.025  limit  as  described  yields  the  following  requirements  for  the  G 
coefficients : 


GQ1>4°lZQl+ZQ2l 


(3.25) 


G  1  >  40 
4 


ZP1ZQ2 


(3.26) 


G-1  >  40 


(3.27) 


G"1  >  40 
3 


11 


(3.28) 
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G  -1  >  40 


zpiV 


11 


(3.29) 


G~  >  4o!z   +  z 
P       '  PI    P2 


(3.30) 


where  the  reciprocals  are  used  in  order  to  express  the  results  in 
familiar  ohmic  dimensions.   Table  3.2  illustrates  the  requirements  on 
the  G  coefficients  at  a  frequency  of  3  Hz  for  the  numerical  examples 
at  the  end  of  Chapter  2. 

TABLE  3.2 

LIMITS  OF  CLASS  1,  2,  AND  3  CONDUCTANCE  FOR  THE  EXAMPLES  OF 

CHAPTER  2 


,-1 

JQ 

,-i 
% 

,-i 


Example  1      Example  2 


130  M 


64  M 


>    2.2  '  1012     2.6  •  109 


>    2.2  '  1012     2.6  •  109 


>    2.2  •  1012     4.4  •  1012 


ohms  at  3  Hz 


>    2.2  •  1012     4.4  •  1012 


130  M 


High  frequency  errors 


130  M 


Although,  in  the  limit,  the  interface  impedances  become  negli- 
gible with  frequency,  this  condition  may  not  occur  before  the  error  due 
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to  parasitic  capacitance  in  the  Y  coefficients  exceeds  permissible 
limits.   Thus,  in  the  analysis  that  follows,  both  components  of  the 
m  parameters  are  included  and  the  Y  coefficients  are  represented  by 
el90S,  where  S  is  the  susceptance  (magnitude)  and  e^90  is  the  90- 
degree  phase  factor  indicating  pure  capacitance.   The  high  frequency 
relative  loading  errors  are 


aQl  =  _eJ9°  C(ZQ1  +  h31  +  hll  +  ZQ2  +  h23)SQ 


+  (ZQ1  +  h31  +  hn)S3  +  (ZQ1  +  h31)S1]         (3.31) 


api  =  "ej9°  Zpih^hl2  t(hn  +  zQ2  +  h23)s. 


-    (ZQ1  +  h31)S1  +  h11Sp]  (3.32) 


ap2   =   e^90  — _ [(zQ2  +  h23)S2   -    (ZQ1  +  h31  +  hu)S3 


husp]  (3.33) 


In  addition  to  the  e^   factor,  the  phase  angle  of  the  total  error  is 
clearly  a  function  of  the  .phase  angles  in  the  interface  impedance 
modified  by  the  h  parameters.   One  is  led  to  the  conclusion  that  ctT 
may  fall  in  any  quadrant  of  the  complex  plane,  depending  on  the 
strengths  of  the  various  factors . 

The  combination  of  Equations  (3.31),  (3.32)  and  (3.33)  to  form 
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c^  produces  six  terms,  each  containing  a  different  S  coefficient. 
These  terms  are  completely  analogous  to  those  of  the  G  coefficients 
studied  earlier  and  may  be  treated  in  the  same  manner.   At  high 
frequencies,  the  terms  of  S  and  Sp  are  more  significant  than  were 
their  low  frequency  counterparts,  but  they  may  still  be  assumed  to 
be  smaller  than  the  other  high  frequency  terms.   Hence  the  0.025  limit 
is  applied  again  to  each  term  and  U  |  is  reasoned  to  be  less  than 
0.05.   Accordingly,  the  requirements  for  the  S  coefficients  are 


SQ  >    4°|Zm  +  hQ1  +  hin  +  ZAO  +  h, 


Ql    31 


'11 


JQ2  '  "231 


(3.34) 


S~x   >    40 


ZP1  +  h12 


1  + 


ZQ2  +  h23 


11 


(3.35) 


S~x  >  40 


Z^  +  h 

P2          13 

ZQ2  +  h23 


11 


(3.36) 


S"1  >  40 


ZP2  +  h13  +  hll 


1  + 


ZQ1  +  h31 


'11 


(3.37) 


S  *  >   40 


ZP1  +  h12  "  hll 


ZQ1  +  h31 


11 


(3.38) 


5P   >  "°   ZP1  +  h12  +  ZP2  +  h 


13 


(3.39) 
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The  results  of  these  relations  for  the  examples  in  Chapter  2  are 
summarized  in  Table  3.3.   To  show  the  effect  of  interface  impedance 
and  the  importance  of  the  h  parameters  at  high  frequencies,  the  limits 
were  also  computed  with  the  interface  impedances  neglected.   The 
limits  are  expressed  in  ohms  at  100  kHz  in  the  upper  part  of  the 
table  and  in  picofarads  of  capacitance  in  the  lower  part,  where 


C  =  (2tt  •  105)_1S 


and  2tt  ■  10  is  the  radian  frequency  at  100  kHz. 


(3.40) 


Transverse  admittance  errors .--Connectors  and  implanted  electrodes 
are  usually  supported  by  means  of  a  potting  material  such  as  acrylic 
or  epoxy  resin.   This  support  introduces  interelectrode  capacitances 
of  between  0.5  and  10.0  picofarads.   The  choice  and  handling  of  support 
materials  requires  special  consideration,  especially  when  limits  on 
the  reciprocal  interelectrode  conductances  are  above  109  ohms. 

With  respect  to  the  results  in  Tables  3.2  and  3.3,  the  limits  in 
Example  1  are  precarious  for  G"1,  G^1,  G"1,  G^1 ,  C4,  C2,  C3  and  Cr 
In  Example  2,  the  limits  are  precarious  for  G^1,  G"1,  C3  and  Ci.   These 
results  demonstrate  that  for  more  reliable  performance  at  both  low  and 
high  frequencies,  it  is  advantageous  to  minimize  the  current  density 
at  the  current  electrodes.   If  a  larger  electrode  were  used  for  Ql  in 
Example  2,  all  of  the  G  and  C  limits  would  have  been  within  convenient 
range . 


Cable  and  amplifier  input  admittance  errors. --The  input  circuit 
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TABLE  3.3 

LIMITS  OF  CLASS  1,  2,  AND  3  SUSCEPTANCE  AND  CAPACITANCE 
FOR  THE  EXAMPLES  OF  CHAPTER  2* 


,-1 


Example  1 Example  2 


120k 


(78k) 


60k 


(40k) 


1.8M 


(800k) 


170k 


(110k) 


1.7M 


(760k) 


110k 


(76k)      ohms  at  100  kHz 


1.8M 


(840k) 


3.4M 


(1.6M) 


1 
,-1 


1.7M 


(720k) 


110k 


(76k) 


3.4M 


(1.6M) 


110k 


(76k) 


13 


(20) 


26 


(39) 


0.90 


(1.9) 


9.0 


(14) 


0.90 


(2.0) 


0.85 


(1.8) 


14 


(20) 


0.46 


(1.0) 


picofarads 


0.95 


(2.2) 


0.46 


(1.0) 


14 


(20) 


14 


(20) 


:Values  in  parentheses  are  obtained  by  neglecting  interface  impedance. 
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of  a  typical  differential  amplifier  may  be  represented  by  three  admittance 
elements,  as  indicated  in  Figure  3.5  (Graham,  1965a).   Two  of  the  elements, 


Figure  3.5.   Input  Circuit  of  a  Differential  Amplifier 

Ya  and  Y5'  are  generally  equal.   In  most  cases,  a  fcot  or  more  of 
shielded  cable  is  required  to  connect  the  amplifier  to  the  P  terminals 
of  the  electrode-tissue  system.   If  the  two  cable  conductors  are 
enclosed  in  a  single  shield,  a  component  of  the  cable  admittance  adds 

to  Ya'  Yb'  and  Yd  to  form  the  maior  Part  of  Y4>  Y2>  and  Yp  (see  Figure 
3.3),  respectively.   If  separate  shields  are  used,  the  cables  only 
contribute  to  Y  and  Jf.. 

The  results  of  Tables  3.2  and  3.3  for  the  elements  of  L,  Y  , 
and  Yp  illustrate  requirements  on  cable  and  amplifier  input  admittance 
parameters.   With  respect  to  the  amplifier,  the  limits  in  Example  2 
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are  not  difficult  to  meet,  especially  with  field  effect  transistor 
input  stages.   In  Example  1,  the  limits  on  G^1  and  G^1  are  borderline, 
while  those  on  C^  and  C2  may  be  very  difficult  to  satisfy  in  practice. 
With  respect  to  the  cables,  the  requirements  in  Example  1  on  C  and  C. 
are  too  stringent  to  meet  unless  special  provision  is  made  to  compensate 
the  cable  capacitance.   One  such  technique  is  described  in  Chapter  5. 

Another  source  of  cable  admittance  is  the  cable  connecting 
terminal  Ql  to  the  exciting  system.   It  usually  constitutes  the  major 
part  of  YQ.   With  reference  to  Table  3.2,  the  requirements  on  G-1  are 
easy  to  meet.   In  Table  3.3,  however,  special  care  may  be  necessary  to 
satisfy  the  requirement  on  CQ  in  Example  1.   The  Ql  cable  is  frequently 
not  shielded. 

Errors  Due  to  Ground  Admittance 

Since  the  error  due  to  ground  admittance  depends  on  how  the 
stray  ground  current  is  distributed  with  respect  to  the  main  current 
field,  a  detailed  analysis  in  the  general  case  is  not  feasible. 
However,  it  is  possible  to  characterize,  in  a  general  way,  the 
susceptibility  of  a  system  to  stray  ground  current.   For  that  purpose, 
assume  that  all  of  the  ground  admittance,  denoted  as  Y  in  Figure 
3.2,  is  concentrated  at  a  point  of  maximum  voltage  magnitude  in  the 
tissue  mass.   For  this  extreme  condition,  the  voltage  and  current  are 
indicated  by  ^(max   and  1^^,  respectively,  and  may  be  expressed  in 
the  following  manner 


y-'.w 
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where 


F 
g 


V 
gmax 


JQ1 


(3.42) 


If  the  loading  errors  are  small,  F  is  primarily  a  function  of 
electrode-tissue  system  parameters  (as  shown  in  the  next  paragraph) 

and  may  be  viewed  as  a  figure  of  merit  for  the  system.   The  larger  F 

g 
is,  the  greater  is  the  ratio  of  stray  ground  current  to  exciting 

current  for  a  given  Y   and  the  greater  is  the  chance  for  error. 

The  maximum  voltage  in  the  tissue  mass  corresponds  to  points 

adjacent  to  the  Ql  electrode  and  is  given  by 

Vgmax  =  VQ  "  ZQ1JQ1  (3-43) 

If  Vq  is  eliminated  by  means  of  Equation  (3.9),  the  expression  for 
Vgmax  becomes 

Vgmax  '  (m31  +  m23  +  mll  "  Zqi^QI  (•.**) 

which  is  valid  when  the  loading  error  in  V  is  small.  Substituting 
Equation  (3.44)  into  Equation  (3.42)  and  replacing  the  m  parameters 
by  their  constituents  from  Equation  (2.22)  yields 

Fg  =  lh31  +  hll  +  ZQ2  +  h23 I  (3-45) 

Let  the  relative  loading  error  in  V  due  to  ground  admittance 
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be  designated  by  Q   (analogous  to  a   ,  a   ,  and  a   ).   Then,  a  may 
be  included  in  ctL  and  Equations  (3.17)  through  (3.21)  may  be  applied 
to  determine  the  corresponding  magnitude  and  phase  errors  in  m, ,  .   It 
is  not  difficult  to  see  that  [a  |  will  be  maximum  when  all  of  the 

to 

ground  admittance  is  concentrated  at  the  point  where  V     occurs . 

gmax 

For  this  condition,  |0  |  is  approximately  equal  to  the  current  ratio 

given  in  Equation  (3.41).   To  illustrate  requirements  for  Y  ,  assume 

g 

|0g|  is  not  to  exceed  0.01.   Then  from  Equation  (3.41),  the  condition 
on  Y  becomes 

© 


|Y-1|  >  100F 
1  g  '       I 


(3.46) 


Table  3.4  summarizes  the  results  of  Equation  (3.46)  at  3Hz  and  100  kHz 
for  the  examples  of  Chapter  2.   As  before,  the  symbols  G,  S,  and  C 
are  used  to  indicate  conductance,  susceptance,  and  capacitance, 
respectively. 


TABLE  3.4 
LIMITS  OF  GROUND  ADMITTANCE  FOR  THE  EXAMPLES  OF  CHAPTER  2- 
Example  1 Example  2 


g 


3-1 
g 


*g 


160M  160k 

(200k)         (100k) 


230k 

(200k) 


100k 


(100k) 


5.8 


16 


(7.8) 


(16) 


ohms  at  3  Hz 


ohms  at  100  kHz 


picofarads 


*Values  in  parentheses  are  obtained  by  neglecting  interface 
impedance . 
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Since  the  extreme  conditions  assumed  in  deriving  the  tabulated 
results  are  not  likely  to  occur  in  practice,  the  values  shown  are 
probably  more  severe  than  necessary.   However,  even  allowing  for  one 
or  two  orders  of  magnitude,  the  requirements  in  Example  1  may  be 
inconvenient  to  satisfy  under  certain  experimental  conditions  (e.g., 
during  surgical  procedures ) . 

Common  Mode  Rejection  Error 

If  the  differential  amplifier  of  Figure  3.5  is  connected  to 
the  test  system  of  Figure  3.2,  the  common  mode  signal  to  the  amplifier 
is  Vj.   The  effective  differential  input  signal  is  given  by 

VPp  =  VP  (1  +  V  (3.47) 


where 


_V2 
"p  =^v7  (3.48) 


and  p  is  the  common  mode  rejection  ratio  of  the  amplifier. 


The  second  experimental  estimation  of  nu,  is  defined  as  m, 
and  includes  both  loading  and  common  mode  rejection  error.  It  is 
given  by 


A   _  VPp 

BU  =  J-^  (3.49) 

Q 
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which,  using  Equations  (3.47)  and  (3.16),  becomes 


A 


mil  =  mil  (1  +  ao)  <3-5°) 


Magnitude  and  phase  errors. --The  relative  error  in  the  magnitude 


of  BL^  as  an  estimate  of  m   is 


«  =  I1  +  apl  "  1  (3.51) 


and  the  phase  error  is 


$  =  /         +  ao)  (3.52) 


For  a  given  |ap|  the  relative  magnitude  error  becomes  maximum  when 
the  angle  of  ap   is  0  or  180  degrees,  and  the  phase  error  becomes 
maximum  when  the  angle  of  o  is  ±(90  +  sin_1|a  |)  degrees.   For  small 


| a  |,  the  maximum  errors  may  be  written 


M  *   *|opl  (3.53) 


with  zero  phase  error,  and 


$  *  ±57. 3  I ap  I  degrees  (3.54) 


-la  I2 


with  approximately  £ —  relative  magnitude  error.   These  expressions 

are  analogous  to  those  found  for  a..      Since  the  common  mode  rejection 
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ratio  is  somewhat  more  controllable  than  the  loading  admittance 

parameters,  it  is  reasonable  to  place  a  more  severe  limit  on  la  I. 

P  ' 

For  the  computations  to  follow,  |o  |  will  be  limited  to  less  than  0.01, 
corresponding  to  magnitude  and  phase  errors  within  ±1  percent  and  ±0.6 
degrees  of  zero,  respectively. 
From  Equation  (2.1) 


!i_^Q   Vl  (3.55) 

Vp  "  Vp   Vp 


which  with  the  small  error  approximation  of  Equations  (3.11)  and. 
(3.12)  becomes 


VJ_   a  ^23  (3.56) 

vp   mu 


Substituting  Equation  (3.56)  into  Equation  (3.48)  and  replacing  m„„ 

^_  o 

by  its  constituents  from  Equation  (2.22)  yields 


M  11 


If  jo  |  is  required  to  be  less  than  0.01,  then  p  must  satisfy  the 
following  inequality 


p   >  100 


ZQ2  +  h23 
hll 


(3.58) 
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Table  3.5  summarizes  the  results  of  Equation  (3.58)  for  the  examples 
in  Chapter  2. 

TABLE  3.5 

MINIMUM  VALUES  OF  COMMON  MODE  REJECTION  RATIO  FOR  THE 
EXAMPLES  OF  CHAPTER  2* 


Frequency 

Example  1 

Example  2 

3  Hz 

3.4M 

(2.0k) 

4.2k 

(200) 

100  kHz 

3.0k 

(2.0k) 

200 

(200) 

*Values  in  parentheses  are  obtained  by  neglecting 
interface  impedance . 

The  requirements  for  Example  1  in  Table  3.5  would  be  extremely 
difficult  to  satisfy  in  practice.   Low  noise  broad  band  differential 
amplifiers  with  high  common  mode  rejection  ratios  and  high  input 
impedances  are  rather  difficult  to  obtain.   In  addition,  special 
precautions  are  necessary  to  realize  the  potential  capabilities  of 
such  amplifiers.   Some  units  are  supplied  with  high  impedance  probes 
which  usually  lead  to  common  mode  rejection  ratios  well  below  that  of 
the  amplifier  alone. 

Current  Control  Error 

Ideally,  one  would  like  to  establish  the  applied  signal  current 
independently  of  the  properties  of  the  electrode-tissue  system  and 
have  available  a  convenient  measure  of  the  current  as  a  reference  in 
measuring  the  voltage  response.   The  simplest  approach  is  to  use  a. 
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voltage  generator  in  series  with  a  large  resistor  as  indicated  in 
Figure  3.6.   If  I  were  determined  entirely  by  E   (within  an  arbitrary 


wwv — — 


to  Ql 


Figure  3.6.   Approximation  of  a  Current  Source 


real  constant),  then  E  could  be  used  as  the  reference  signal.   Further- 
more, this  condition  would  permit  controlled  nonsinusoidal  currents 
to  be  applied  when  information  is  desired  at  several  frequencies 
simultaneously.   The  errors  of  this  method  are  analyzed  in  the 
following  paragraphs.   The  results  are  used  again  in  Chapters  4  and 
5  to  evaluate  the  new  techniques  presented  there. 

In  a  physical  system,  parasitic  capacitance  will  generally 
reduce  the  effective  value  of  R  and  produce  phase  shift  at  the  high 
frequencies.   To  account  approximately  for  this  effect  in  the  present 
analysis,  let  the  resistor  value  be  multiplied  by  the  factor 
(1  +  jf/f  )   •   This  form  results  if  the  parasitic  capacitance  is 
assumed  to  act  entirely  in  parallel  with  R  .   Assume  that  E  serves 
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as  the  reference  signal  for  measuring  the  voltage  response,  V  ,  of 
the  electrode-tissue  system.   Then  the  experimental  estimation  of  I 
is  defined  by 


Vr*  (3.59) 


This  expression  may  also  be  written 


A  I 

Xn  =  Xr,    (-tTTFTp-  +  O  (3.60) 


Q     Q    l+jf/fQ     Q' 


where 


■v.  -°  ^ 


The  third  experimental  estimation  of  m   ,  which  includes  loading, 
common  mode  rejection,  and  current  control  error,  is  defined  by 


&     V 


£    =  _pP 

A 

I 


Bn  =  jr^-  (3.62) 


Q 


where  V   is  as  defined  in  Equation  (3.47).   Using  Equations  (3.60) 
and  (3.49),  Equation  (3.62)  may  be  written 


m   =  . ± (3.63) 


11 


1  +  3f/f0 


+  Of 
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It  is  clear  that  the  current  control  error  arises  from  the  two  terms 
in  the  denominator  of  Equation  (3.63).   As  shown  in  the  following 
sections,  low  frequency  errors  are  determined  by  a  ,  while  high 
frequency  errors  are  determined  by  the  parasitic  capacitance  term. 

Low  frequency  errors 

For  small  loading  errors , 


IQ  -  IQ1  (3.64) 


and  from  Equation  (3.9) 


V0 

-2—  -   m   +  m   +  m  (3.65) 

IQ1    31    23    11 


Substituting  Equations  (3.64)  and  (3.65)  into  Equation  (3.61)  and 
replacing  the  m  parameters  by  their  constituents  from  Equation  (2.22) 
yields 


ZQ1  +  h31  +  hll  +  ZQ2  +  h23  /0  „, 

aQ (3.66) 


At  low  frequencies,  the  effective  interface  impedances  predominate 
and  drive  the  angle  of  a_  well  into  the  fourth  quadrant  of  the  complex 
plane.   As  the  frequency  increases,  |a  |  decreases  and  the  angle  of 
a,-,  moves  toward  zero. 

Since  the  ratio  f/f  in  Equation  (3.63)  is  generally  negligible 


at  frequencies  below  10  kHz,  the  low  frequency  relative  error  in  th 


A 


magnitude  of  m   as  an  estimate  of  in, ,  is 
J-j-  11 


M  = 


|l  +  a 


and  the  phase  error  is 


*  =  -       Al  +  o  ) 
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error 

in  the 

(3.67) 

(3.68) 

The  vector,  1  +  a  ,  is  indicated  in  Figure  3.7, 


IV 


Figure  3.7.   Phasor  Relations  of  a     in  the  Complex  Plane 

It  is  evident  that  M  is  negative  and  $  is  positive. 

If  |a  |  is  limited  to  less  than  0.01,  then  from  Equation  (3.66) 
R  must  satisfy  the  following  inequality 
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Rq  >100|ZQ1  +  h31  +  hu  +  Zq2  ♦  h23|  (3.69) 

The  limits  for  the  examples  in  Chapter  2  are  summarized  in  Table  3.6, 
where  the  100  kHz  calculation  is  included  for  comparison. 

TABLE  3.6 

MINIMUM  VALUES  OF  CURRENT  SOURCE  RESISTANCE* 
FOR  THE  EXAMPLES  OF  CHAPTER  2 

Frequency Example  1      Example  2 


3  Hz  320.M  160M 

(200k)         (100k) 

ohms 

100  kHz       290k  150k 

(200k)         (100k) 


•'Values  in  parentheses  are  obtained  by  neglecting 
interface  impedance . 


High  frequency  errors 

At  high  frequencies,  aQ  generally  becomes  negligible,  so  that 
Equation  (3.63)  may  be  written 


m 


m   (1  +  jf/f  )  (3.70) 


11    11 

It  is  evident  that  parasitic  capacitance  across  the  source  resistor 
causes  positive  magnitude  and  phase  errors  at  high  frequencies.   The 
3  d.b.  frequency,  f  ,  may  be  expected  to  vary  more  or  less  inversely 
with  the  resistance  value  so  that  the  larger  R  mUst  be  to  satisfy  the 
limit  on  cxq,  the  worse  the  frequency  response  becomes.   For  example,  a 
1  megohm  carbon  resistor  might  have  an  f  at  about  200  kHz.   To 
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avoid  a  phase  error  greater  than  6  degrees,  the  frequencies  for 

impedance  measurements  would  have  to  be  less  than  0.1  f   or  20  kHz 

o  * 

in  this  example.   For  the  values  in  Table  3.6,  the  useful  frequency 
range  is  proportionately  smaller. 

To  conclude  the  analysis  of  loading,  common  mode  rejection,  and 
current  control  error,  it  is  instructive  to  consider  all  three  error 
factors  explicitly  in  one  expression.   Thus,  combining  Equations  (3.17), 
(3.50)  and  (3.63),  the  third  estimation  of  m, ,  may  be  written 

mu  -   mn(l  ♦  jf/f0)  u  -  ah){l  I   ^1  ;  jf/fQ])         (3-?1) 

Equation  (3.71)  may  be  simplified  if  the  denominator  factors  are 
replaced  by  equivalent  series  expansions  in  the  numerator.   Then,  if 
each  a  is  assumed  to  be  small  (e.g.,  less  than  0.05),  the  factors  in 
the  numerator  may  be  multiplied  together,  and  the  resulting  second 
order  terms  may  be  neglected.   This  procedure  leads  to  the  following 
approximate  expression  for  the  third  estimation  of  m,  : 


11  S  mll(l  +  ^f/fo)(l  +  ao  +  aL  "  aQ[1  +  if/fo^  (3-72) 


The  most  uncertain  factor  in  Equation  (3.72)  will  generally  be 
This  follows  from  the  fact  that  transverse  and  ground  admittance 
components  are  easily  influenced  by  indeterminate  properties  of  the 
experimental  environment.   In  contrast,  common  mode  rejection  and 
current  control  capabilities  are  relatively  definite  characteristics 


aL 


m 
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of  the  instrumentation  system  which  can  be  evaluated  for  a  given 
electrode-tissue  system.   The  effect  of  the  frequency  factor  multiplying 

, ,  can  usually  be  eliminated  by  introducing  suitable  compensation 
at  the  input  signal  generator,  or,  more  conveniently,  in  the  gain  of 
the  response  amplifier.   The  frequency  factor  multiplying  a  cannot  be 
eliminated  and  its  importance  depends  on  a_. 

For  convenience,  assume  the  frequency  factor  multiplying  m-,  ,  has 
been  compensated  and  the  frequency  factor  multiplying  a 0  is  included 
with  aQ.   The  total  relative  error  is  then  given  approximately  by  the 
algebraic  sum  of  the  three  a  terms ,  which  may  be  reduced  to  magnitude 
and  phase  error  in  a  manner  similar  to  that  for  a   in  Equations  (3.50) 
and  (3.54).   The- relative  importance  of  the  a  terms  depends  on  the 
properties  of  the  instrumentation  system  and  the  electrode-tissue 
system  used.   Each  case  may  be  evaluated  in  the  light  of  the  analyses 
presented  in  this  and  succeeding  chapters. 

Noise  in  the  Classical  System 

A  detailed  analysis  of  the  noise  problem  is  beyond  the  scope 
of  the  present  work.   This  section  gives  a  brief  account  of  the  major 
types  of  noise  found  in  biological  impedance  measurements.   It  also 
serves  as  a  basis  for  the  noise  analysis  presented  in  Chapter  4. 

In  linear  studies  of  tissue  impedance,  the  response  voltage 
is  normally  much  less  than  one  millivolt  and  frequently  less  than  one 
hundred  microvolts,  in  amplitude.   Noise,  which  may  be  significant  at 
these  levels,  originates  in  the  instrumentation  system,  in  the 
external  environment,  and  in  the  electrode-tissue  system  itself.   These 
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classes  of  noise  are  considered  in  the  following  sections  with  the  aid 
of  Figure  3.8. 


Figure  3.8.   Noise  Sources  in  the  Classical  Method 
Instrumentation  Noise 

Instrumentation  noise  consists  primarily  of  shot  and/or  thermal 
noise  which  originates  in  the  active  and  passive  components  of  the 
initial  amplifier  stages.   It  may  be  distinguished  from  other  noise 
by  its  broad  bandwidth  (limited  primarily  by  the  method  of  observa- 
tion) and  by  its  more  limited  and  consistant  range  of  amplitude.   Part 
of  the  noise  may  be  described  by  means  of  noise  current  sources 
connected  between  the  amplifier  input  terminals  and  from  each  input 
to  ground  (not  shown  in  Figure  3.8).   The  rest  may  be  described  by 
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means  of  a  noise  voltage  source  in  series  with  the  input  terminals 
(N  in  Figure  3.8).   Since  most  of  the  noise  energy  occurs  at  frequencies 
above  10  kHz,  where  the  terminal  impedances  of  the  electrode-tissue 
system  are  not  especially  large,  the  noise  current  sources  are  usually 
less  significant  than  the  noise  voltage  source.   For  bandwidths  up  to 
1  MHz,  the  equivalent  input  noise  of  most  amplifiers  is  less  than  30 
microvolts  peak-to-peak  (as  observed  on  an  oscilloscope  trace). 

Environmental  Noise 


em 


Noise  originating  in  the  external  environment  enters  the  system 
by  means  of  electrostatic  and/or  electromagnetic  coupling.   The  syst 
is  most  sensitive  to  this  noise  at  the  cables  of  PI  and  P2,  where  the 
impedance  level  is  high,  the  signal  level  is  low,  and  the  exposure  is 
relatively  great.   Environmental  noise  may  also  enter  throughout  the 
volume  of  the  impedance  subject  and  by  means  of  the  Ql  cable.   The 
discussion  that  follows  is  concerned  only  with  the  PI  and  P2  cables. 

In  addition  to  shielding  the  cables,  electrostatic  noise  may  be 
reduced  by  arranging  them  symmetrically  with  respect  to  possible  noise 
sources,  which  takes  advantage  of  the  differential  nature  of  the 
amplifier.   To  reduce  electromagnetic  coupling,  the  loop  area  between 
cables  should  be  kept  to  a  minimum.   Both  objectives  may  be  achieved 
by  twisting  the  cables  together.   Electrostatic  noise  also  depends  on 
the  impedance  levels  in  the  electrode-tissue  system  and  since  it 
usually  is  a  more  prevalent  problem,  an  analysis  follows. 

Electrostatically  coupled  noise  may  be  described  by  means  of 
noise  currents  feeding  the  PI  and  P2  cables.   For  simplicity,  assume 
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the  currents  are  equal  as  Indicated  by  I   in  Figure  3.8.   If  V 

x  Px 

and  V,>x  designate  the  differential  and  common  mode  voltages  created  by 
I  ,  then  the  effective  differential  noise  voltage  is 


'2x 


vpxP  =  vpx  +  -f-  (3-73) 


where  p  is  the  common  mode  rejection  ratio  of  the  differential 

amplifier. 

Since  the  external  circuit  impedances  to  ground  at  PI  and  P2 

are  quite  high,  the  noise  components  of  I-,  and  IDO  are  essentially  I  . 

riri  x 

Using  this  valve  in  Equation  (2.3)  yields  for  the  noise  component  of 
VP 


VPx  =  Ix(m12  "  m13>  0. 74) 


The  noise  component  of  V-  may  be  written  according  to  Equation  (2.1)  as 


V2x  =  VQx  -  Vlx  "  VPx  (3.75) 


If  V   ,  Vlx,  and  Vpx  are  replaced  according  to  Equation  (2.3),  and 
the  dependent  parameters,  nun  and  n^j  are  eliminated  by  means  of 
Equations  (2.4)  and  (2.5),  the  expression  for  V„  becomes 

V2x  =  V"l3  +  2m23  +  2m33>  (3-76) 

Substituting  Equations  (3. 74)  and  (3.76)  into  Equation  (3.73) 
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and  replacing  the  m  parameters  by  their  constituents  from  Equation 
(2.22),  the  effective  differential  noise  voltage  may  be  written 


Pxp 


Ql 


PI 


zp2  +  h12  -  h13 


11 


Zp2  +  2Zq2  +  h13  +  2h23  +  2h33 


Ph 


11 


(3.77) 


where  the  expression  is  divided  by  Vp  (as  given  in  Equation  (3.8))  to 
obtain  the  noise  in  relative  terms.   Note  that  I_,  is  the  nominal 
exciting  current  and  Vp  is  the  nominal  voltage  response.   As  a  rule, 
the  cancellation  occurring  in  the  numerator  of  the  left  term  in  the 
brackets  is  neither  complete  nor  easily  controlled  (even  in  symmetrical 
systems)  and  cannot  be  used  as  a  reliable  index  of  noise  reduction. 
This  statement  is  especially  true  at  noise  frequencies  where  the 
interface  impedances  are  large.   The  right  term  in  the  brackets  is  not 
likely  to  be  important  unless  p  is  abnormally  small. 

It  is  instructive  to  determine  the  requirement  on  the  signal- 
to-noise  current  ratio  to  meet  a  given  limit  on  the  signal-to-noise 
voltage  ratio  in  the  response.   Assume  that  | Vp |  is  required  to  be 


greater  than  10 |V 


Pxp 


Then  according  to  Equation  (3.77),  the  currents 


must  satisfy  the  following  inequality; 


Ql 


>  10 


ZP1  ~  ZP2  +  h12  "  h13 
hll 


Zp2  +  2ZQ2  +  h13  +  2h23  +  2h33 


Ph 


11 


(3.78) 
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Since  an  upper  limit  is  placed  on  I_,  by  the  requirements  of  linearity 

and  passivity  in  the  tissue  to  be  measured,  the  inequality  can  not  be 

satisfied  by  increasing  I   indefinitely.   Thus,  special  precautions 

are  often  necessary  to  minimize  I  .   The  greatest  problem  is  usually 

presented  by  sixty-cycle  power  lines  and  equipment.   In  many  situations, 

Equation  (3.78)  cannot  be  satisfied  for  this  type  of  noise,  unless  the 

potential  cables  are  shielded  to  reduce  I  . 

x 

Electrode-Tissue  Noise 

Noise  of  the  electrode-tissue  system  originates  independently 
in  the  tissue  mass  and  at  each  of  the  electrode-tissue  interfaces. 
These  two  classes  of  noise  are  explained  in  the  following  paragraphs. 

Tissue  noise. — Noise  potentials  reflect  the  spontaneous  and 
induced  electrochemical  activity  of  living  cells .   It  is  well  known  that 
potentials  of  myogenic  and/or  nervous  origin  occur  in  most  body  tissues 
and  differ  in  form  and  amplitude  from  place  to  place.   In  the  brain, 
for  example,  these  potentials,  called  the  electrocortigram,  may  make 
impedance  measurements  especially  difficult.   The  electrocortigram 
appears  as  a  random,  rhythmic,  relatively  continuous  activity,  which 
changes  in  character  and  amplitude  with  the  physiological  state  of  the 
individual.   Maximum  peak  amplitudes  generally  range  from  100  to  150 
microvolts.   Most  of  the  energy  in  the  electrocortigram  lies  in  the 
low  frequency  band  between  1  and  30  Hz . 

Interface  noise. — Noise  potentials  also  arise  from  the 
spontaneous  electrochemical  activity  of  the  metal-electrolyte  interface 
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and  include  d.c.  and  a.c.  components  (Flasterstein,  1966a,  b).   In 
theory,  thermodynamics  predicts  the  occurrence  of  a  d.c.  equilibrium 
potential  between  the  metal  and  the  electrolyte.   In  practice,  a 
number  of  factors  relating  to  the  inherent  inhomogeneity  of  the 
interface  prevent  the  attainment  of  a  true  equilibrium.   Consequently, 
the  d.c.  level  is  unstable  and  varies  randomly  in  time.   The  d.c. 
potential  of  two  apparently  identical  electrodes  in  electrolyte  may 
differ  by  more  than  several  tenths  of  a  volt  with  variations  (positive 
and  negative)  of  the  same  order.   Superimposed  on  this  slow  noise, 
there  often  occur  faster  a.c.  and  pulse  potentials  with  frequency 
components  in  the  electrocortigram  spectrum.   The  peak  amplitude  of 
the  a.c.  and  pulse  potentials  may  reach  several  hundred  microvolts. 

It  is  evident  from  Figure  3.8  that  each  input  of  the  differential 
amplifier  receives  noise  from  the  Q2  interface  and  the  corresponding  P 
interface  in  addition  to  noise  generated  in  the  tissue  mass.   The 
differential  signal  includes  the  noise  of  both  P  interfaces.   To 
prevent  the  instrumentation  system  (amplifiers  and/or  other  equipment) 
from  being  overdriven,  d.c.  imbalances  and  slow  fluctuations  should  be 
rejected  in  the  early  stages  of  amplification.   For  example,  a  system 
with  any  member  set  to  near-full  range  for  signals  on  the  order  of 
100  microvolts,  would  have  to  reject  99  percent  or  more  of  the  noise 
to  assure  reliable  operation.   A  versatile  system  should  be  able  to 
handle  signals  down  to  0.1  Hz  and  have  the  option  of  d.c.  coupling. 
To  achieve  these  aims,  without  degrading  amplifier  performance,  is 
likely  to  prove  difficult  because  of  the  large  filter  components 
required.   One  may  also  object  to  the  inconvenience  of  making  modifi- 
cations inside  the  amplifier. 
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When  the  impedance  information  of  interest  occurs  at  frequencies 
well  above  the  spectrum  of  electrode-tissue  noise,  a  low  frequency 
filter  usually  suffices  to  establish  an  adequate  signal-to-noise  ratio 
before  application  to  the  measuring  instruments .   If  signal  frequencies 
overlap  the  noise  spectrum,  then  special  averaging  techniques  may  be 
necessary  to  improve  the  signal-to-noise  ratio.   These  techniques  make 
use  of  the  random  character  of  the  noise  and  the  deterministic  nature 
of  the  signal,  and  usually  require  extended  observation  times. 


CHAPTER  4 


THE  VIRTUAL-GROUND  CONFIGURATION 


Introduction 


The  analysis  of  the  classical  method  in  Chapter  3  shows  that 
to  achieve  a  given  experimental  accuracy,  the  common  mode  rejection 
ratio,  the  loading  admittance  and  the  current  source  resistance 
must  satisfy  requirements  that  depend  on  the  parameters  of  the 
electrode-tissue  system.   For  small  electrode  systems  (for  example 
in  animal  experiments  in  vivo)  these  requirements  are  often  too 
severe  to  satisfy  in  practice.   This  statement  is  illustrated  by  the 
numerical  results  in  Chapter  3  for  the  examples  of  Chapter  2.   Another 
problem  which  lacks  a  convenient  solution  in  the  classical  method  is 
presented  by  electrode-tissue  noise. 

In  the  virtual-ground  configuration,  described  in  this  chapter, 
the  parameters  of  the  electrode-tissue  system  continue  to  influence 
the  experimental  accuracy.   However,  the  requirements  on  instrumenta- 
tion parameters  are  much  easier  to  meet  than  in  the  classical  method. 
The  proposed  system  also  provides  a  convenient  means  of  controlling 
the  effect  of  electrode-tissue  noise. 

The  basic  connections  of  the  virtual-ground  configuration  are 
illustrated  in  Figure  4.1.  A  high  gain  inverting  amplifier  is  shown 
connected  between  electrodes  P2  and  Q2  of  the  electrode-tissue  system. 
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Figure  4.1.   Basic  Connections  of  the  Virtual-Ground  Configurati 


on 


The  signal  current  is  injected  at  Ql  and  flows  through  the  tissue  to 
Q2,  where  it  returns  to  ground  via  the  output  stage  of  the  amplifier. 
Since  the  gain  is  high,  Vp2  is  much  smaller  than  V    so  that  the 
latter  tracks  V2  of  the  e.lectrode-tissue  system  (with  inverted  sign). 
If  the  gain  is  sufficiently  high,  Vp2  will  also  be  negligible  compared  to 
Vp  and  the  voltage  response  may  be  measured  between  PI  and  ground  with 
a  single-ended  amplifier.   Otherwise,  the  response  must  be  measured 
differentially  between  PI  and  P2.   Both  cases  are  analyzed  in  this 
chapter. 
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The  most  suitable  amplifiers  for  the  virtual-ground  connection 
belong  to  a  class  of  devices  known  as  operational  amplifiers.   Nominal 
open-loop  gain  and  phase  shift  characteristics  for  the  type  of 
operational  amplifier  of  interest  here  are  sketched  in  Figure  4.2 
(Philbrick,  1966).   The  figure  also  shows  common  ranges  for  the  gain 
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Figure  4.2.   Nominal  Gain  and  Phase  Shift  Characteristics 
of  Operational  Amplifiers 


magnitude,  the  break  frequency  and  the  unity  gain  frequency.   For  an 
inverting  amplifier,  an  additional  180  degrees  must  be  added  to  the 
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phase  plot  to  give  the  total  phase  shift  from  input  to  output.   It 
is  evident  that  the  total  internal  phase  lag  of  the  amplifier  does 
not  reach  360  degrees  until  well  beyond  the  unity  gain  frequency. 
This  characteristic  of  operational  amplifiers  makes  it  possible  to 
construct  stable  wideband  circuits  having  tight  negative  feedback 
loops.   The  matter  of  stability  is  considered  in  a  later  section. 

The  essential  elements  for  the  analysis  of  the  virtual- ground 
configuration  are  indicated  in  Figure  4-.  3.   The  noise  sources  (Ix,  Nv, 
and  Nr)  are  treated  in  a  separate  section  and  should  be  neglected 
unless  otherwise  noted. 

Common  Mode  Rejection  Error 

Single-Ended  Operation 

With  reference  to  Figure  4.3,  the  input  signal  to  the  single- 
ended  amplifier  at  PI  is 


VP1  =  VP(1  +  M  C».l) 


where 


vP2 


Sv   =  V^  (^-2) 


Note  that  gy  plays  the  same  role  here  as  a  does  in  Equation  (3.47) 
By  inspection,  from  Figure  4.3 


VP2  =  V2-KvVp2+Z0Ig2  (4.3) 

where  -Kv  is  the  complex  open-loop  unloaded  gain  and  ZQ  is  the  output 
impedance  of  the  virtual-ground  amplifier.   For  small  loading  error 
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IQ2  *  IQ1  (4.4) 

where  IQ1  is  as  defined  in  Figure  3.2.   Combining  Equations  (4.3)  and 
(4.4),  solving  for  Vp2  and  substituting  the  result  into  Equation  (4.2) 
yields 


V9    Z 
2  +   ° 


vP   vP/IQl 

K   =   —  (4.5) 

1  +  Kv 

Using  the  small  loading  error  approximations  of  Equations  (3.8) 
and  (3.56),  and  replacing  m23  by  its  constituents  from  Equation  (2.22), 
Equation  (4.5)  may  be  written 

m23  +  Zo   Z0?+h23+Zn 
Pv  "  (l+Kv)m1:L  "   (1+Kv)hn  (4.6) 

Note  that  the  factor  (l+K,)  plays  the  same  role  here  as  the  rejection 
ratio  of  the  differential  amplifier  does  in  Equation  (3.57).   As  a 
rule,  the  unity  term  is  negligible  in  comparison  to  K^. 

If  |ey|  is  restricted  to  less  than  0.01,  then  by  Equation 
(4.6),  the  amplifier  must  satisfy  the  following  inequality 


|Ky|  >  100 


ZQ2+h23+Zo 


hll 


(4.7) 


Except  for  the  term  Zq,  this  relation  places  the  same  requirement  on 
|Kv|  as  was  found  for  |p|  In  Equation  (3.58).   It  is  evident  that  the 
gain  requirement  is  reduced  if  ZQ  is  negligible  compared  to  ZQ2+h23. 
In  addition,  the  smaller  Zq  is  at  high  frequencies,  the  more  parasitic 
capacitance  the  amplifier  can  tolerate  at  its  output  terminal  without 
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becoming  unstable.   The  results  of  Table  3.5  are  directly  applicable 
to  Equation  (4.7)  if  |p|  is  replaced  by  | Kv |  and  ZQ  is  neglected. 

The  unity  gain  frequency  of  wideband  operational  amplifiers  is 
commonly  between  30  MHz  and  100  MHz.   Since  the  gain  rolls  off  at 
about  6  db  per  octave,  it  typically  falls  in  the  range  of  300  to  lpOO 
at  100  kHz.   Amplifiers  with  unity  gain  frequencies  near  or  exceeding 
100  MHz  are  likely  to  incur  formidable  stability  problems  in  the 
virtual-ground  connection  and  may  not  be  practical  for  general  appli- 
cations (see  section  on  stability).   If  lpOO  is  taken  as  a  practical 
upper  limit  on  |kJ  at  100  kHz,  Equation  (4.7)  may  be  rearranged  to 
give 


<  10  (4.8) 


This  relation  restricts  the  type  of  electrode-tissue  systems  which 
can  be  measured  by  the  single-ended  method  with  better  than  1  percent 
accuracy  at  100  kHz.   With  reference  to  the  examples  of  Chapter  2, 
only  Example  1  exceeds  the  limit  with  a  parameter  ratio  of  30. 

Differential  Operation 

When  the  requirements  on  | Ky |  cannot  be  met  with  single-ended 
operation,  a  differential- amplifier  may  be  employed  at  PI  and  P2. 
With  reference  to  Figure  4.3,  the  common  mode  signal  to  the  differential 
amplifier  would  be  Vp2.   The  effective  differential  input  signal  would 
be 

VPP  =  VP  (1  +  V  (4'9) 


ZQ2 

+ 

h23 

+ 

Zo 

hll 
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where 

V 


P2 
PA  -  —  (4.10) 


P 


pv 


and  p  is  the  common  mode  rejection  ratio  of  the  differential  amplifier. 
Note  that  6p  plays  the  same  role  in  Equation  (4.9)  as  6  does  in 
Equation  (4.1)  and  ap  does  in  Equation  (3.47).   It  follows  from 
Equations  (4.10)  and  (4.2)  that 


?p  =  f-  (4.11) 


which,  using  Equation  (4.6),  becomes 


z  +h23+z 

5p  =  -t ; (4.12) 

M   p(l+K  )h 

v  11 


Equation  (4.12)  is  to  be  compared  with  Equations  (3.57)  and  (4.6). 
Placing  the  0.01  limit  on  |fi  |  in  Equation  (4.12)  yields  the 


following  requirement  on  p  and  K 


|pK  |  >  100 


ZQ2+h23+Zo 


hll 


(4.13) 


where  the  unity  term  is  omitted  as  justified  earlier.   This  expression 
is  to  be  compared  with  Equations  (3.58)  and  (4.7). 

The  results  in  Table  3.5  apply  directly  to  Equation  (4.13) 
if  |p|  is  replaced  by  |pKj  and  Z  is  neglected.  Recall  that  the 
classical  differential  system  and  the  single-ended  virtual-ground 
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system  are  not  satisfactory  in  the  case  of  Example  1.   For  the  present 
system,  the  requirements  for  Example  1  are  easily  satisfied  with 
relatively  modest  amplifiers.   As  an  illustration,  consider  a  virtual- 
ground  amplifier  with  | Ky [  equal  to  2  •   10 5  at  3  Hz  and  300  at  100 
kHz.   Then,  |p|  need  only  exceed  17  and  10  at  the  two  frequencies, 
respectively. 

Loading  Error 

The  four  classes  of  loading  admittance  defined  in  Chapter  3 
also  occur  in  the  virtual-ground  system.   With  respect  to  the 
transverse  admittance  (class  1)  and  the  component  of  cable  and 
amplifier  input  admittance  (classes  2  and  3)  that  may  occur  between 
PI  and  P2,  the  loading  error  analysis  of  Chapter  3  and  the  results 
of  Tables  3.2  and  3.3  apply  directly.   These  admittances  may  be 
assumed  to  be  included  in  the  four  terminal  block  of  Figure  4.3.   The 
remainder  of  the  loading  admittances  (class  4  and  most  of  classes  2 
and  3)  is  referenced  to  ground  and  requires  special  consideration. 

Cable  and  Amplifier  Input  Admittance 

Cable  and  amplifier  input  admittance  is  indicated  by  YQ1 ,  Ypi, 
and  Yp2  in  Figure  4.3.   From  the  figure,  the  corresponding  loading 
currents  are 

IQ1  =  VQ1YQ1  (4.14) 

^l  =  VP1YP1  (4.15) 
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Ip2   =   Vp2Yp2  (4.16) 


and 


m31 


VQ1   =   VVVP2   a    (~7     +   1   +    6V)VP  (4.17) 


rall 


OR 


where  V±   has  been  replaced  by  the  small  loading  error  approximati 
of  Equation  (3.11)  and  Vp2  has  been  eliminated  by  means  of  Equation 
(4.2). 

Using  Equations  (4.1),  (4.2),  and  (4.17),  the  expressions  for 
the  loading  currents  may  be  written 


it     m3l 
^1*  (i^+1  +  Wqi  (4.18) 


Ip!  =  (1  +  6y)VpYpi  (4.19) 


JP2  =  Wp2  (4.20) 


If  """Ol'  IP1  '  and  lV0   are  considered  as  additions  to  I '  ,  I '   and  T ' 
v     x  cz  Ql   Pi'     "P2 

in  Equation  (3.4),  respectively,  the  new  relative  error  terms  to  be 
added  to  aQ1,  opi,  and  ap2  in  Equations  (3.13),  (3.14),  and  (3.15)  are, 

«         -   ""ll^l        - 

V   =  ~Tp S"Cm31  +    (l+ev)mll]YQi  C+-2D 

_mi2IPl 
6P1   =  — "   "m12d  +   3y)Ypi  (4.22) 
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t" 


13  p2 
3P2  ~  —\T-   =  mi3BvYP2  C^-23) 


respectively.   These  terms  may  now  be  analyzed  following  the  method 
for  the  loading  errors  (i.e.,  the  a)  in  Chapter  3. 
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If  the  magnitudes  of  gQi,  Bpi,  and  Bp2  are  limited  to  les 
than  0.025,  then,  from  Equations  (4.21),  (4.22),  and  (4.23),  the 
requirements  on  the  loading  admittances  become, 


IYq1|  >  4°lZQl+h31+hll+Bvhlll  (*•») 


l*p£l  >  ^0|zpi+h12|  | i+e  |  (4.25) 


|y;J|  >  4o|zp2+h13|  |ev|  (4>26) 


where  the  m  parameters  have  been  replaced  by  their  constituents.   Note 
that  for  single-ended  amplification,  | gy  |  is  normally  negligible 
compared  to  unity.   With  differential  amplification,  |g  |  may  not  be 
negligible,  but  it  will  almost  certainly  be  less  than  unity  in  a  well- 
designed  system.   Comparing  Figure  4.3  with  Figure  3.3,  it  is  evident 
that  YQ1,  Ypi,  and  Y^  correspond  to  the  cable  and  amplifier  components 
of  YQ>  Yi+>  and  Y2>  respectively.   Accordingly,  Equations  (4.24),  (4.25) 
and  (4.26)  are  to  be  compared  at  low  frequencies  to  Equations  (3.25), 
(3.26),  and  (3.27);  and  at  high  frequencies  to  Equations  (3.34),  (3.35), 
and  (3.36),  respectively. 

To  illustrate  the  admittance  requirements  of  Equations  (4.24), 
(4.25)  and  (4.26)  for  the  examples  at  the  end  of  Chapter  2,  it  is 
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first  necessary  to  determine  suitable  values  for  g  .   For  this  purpose, 
assume  Zq  is  negligible  and  let  Ky  be  2  •  105e']'°  at  3  Hz  and  300e"]'90 
at  100  kHz.   The  electrode-tissue  system  parameters  are  given  in  Tables 
2.3  and  2.4.   The  values  of  3y  have  been  calculated  from  Equation  (4.6) 
and  are  given  in  Table  4.1.   It  is  evident  that  for  Example  1,  |g  | 

TABLE  4.1 
VALUES  OF  B  FOR  THE  EXAMPLES  OF  CHAPTER  2* 

Frequency Example  1 Example  2 

3  Hz  .  0.17e~^72         0.00021e~3'69 

(O.OOOleJ'0)  (O.OOOOlei0) 

100  kHz  O.lOe^58  0.0067e^90 

(0.067eJ90)  (0.0067ei90) 


*Values  in  parentheses  are  obtained  by  neglecting  interface 
impedance . 

is  too  large  for  single-ended  amplification  to  be  used.   This  case  was 
noted  in  the  section  on  differential  operation  where  it  was  found  that 
the  rejection  ratio  of  the  differential  amplifier  would  only  have  to 
be  greater  than  17  at  3  Hz  and  greater  than  10  at  100  kHz. 

The  limits  of  cable  and  amplifier  input  admittance  for  the 
examples  of  Chapter  2  are  given  in  Table  4.2  which  is  based  on  the 
values  of  Sy  in  Table  4.1.   If  the  limits  in  Table  4.2  are  compared 
with  the  corresponding  limits  in  Tables  3.2  and  3.3,  it  is  evident 
that  the  requirements  on  cable  and  amplifier  input  admittance  have 
been  significantly  relaxed. 
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TABLE  4.2 

LIMITS  OF  CABLE  AND  AMPLIFIER  INPUT  ADMITTANCE  FOR  THE 
EXAMPLES  OF  CHAPTER  2* 


> 

Example  1 

64M 

(40k) 

Example  2 

G-1 

Qi 

64M 

(40k) 

gpi 

> 

64M 

(38k) 

64M 

(38k) 

ohms  at  3  Hz 

GP2 

> 

11M 

(3.8) 

13k 

(.38) 

S"1 
Ql 

>  • 

58k 

(40k) 

58k 

(40k) 

S"1 

bpl 

> 

58k 

(38k) 

58k 

(38k) 

ohms  at  100  kHz 

s-1 

P2 

> 

5.8k 

(2.6k) 

380 

(260) 

CQ1 

< 

27 
(39) 

27 
(39) 

CP1 

< 

27 
(41) 

27 
(41) 

picofarads 

cp2 

< 

270 

(600) 

4,000 

(6,000) 

••Based  on  3y  from  Table  4.1.   Values  in  parentheses  are 
obtained  by  neglecting  interface  impedance . 
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Ground  Admittance 

As  in  Chapter  3,  assume  that  all  of  the  ground  admittance, 
denoted  by  Yg  in  Figure  4.3,  is  concentrated  at  a  point  of  maximum 
voltage  magnitude  in  the  tissue  mass.   In  the  virtual-ground  system, 
the  maximum  voltage  may  occur  at  points  adjacent  to  the  Ql  and/or  the 
Q2  electrode.   For  points  adjacent  to  Ql,  the  figure  of  merit  defined 
by  Equation  (3. 42)  is 


■gl 


Wqi 


Ql 


m31+(l+E3v)rairZQl 


(4.27) 


where  V   has  been  eliminated  by  means  of  Equations  (4.17),  and  (3.8) 
For  points  adjacent  to  Q2,  the  figure  of  merit  is 


g2 


-KvVP2+(WZo^2 


'Ql 


-K  g  m,  +Z  +Z 
v  v  11  Q2  o 


(4.28) 


where  Vp2  and  1^  have  been  eliminated  by  means  of  Equations  (4.2), 
(3.8),  and  (4.4). 

Replacing  the  m  parameters  by  their  constituents  from  Equation 
(2.22),  and  substituting  for  $v  from  Equation  (4.6),  Equations  (4.27) 
and  (4.28)  may  be  written 


Fgl  =  lh31+hll+6vh 


11 


(4.29) 


and 


g2 


"h23  + 


WZo 
Kv 


(4.30) 
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It  is  evident  that  both  F   and  F 


are  smaller  than  F  of  the  classi- 


gl      g2 

Equation  (3.46) 
■gl  ""**  Lg2   in  Place  of  Fg  for  the  examples  of  Chapter  2,  are 
summarized  in  Table  4.3.   As  before,  ZQ  is  neglected  and  ^  is 


cal  method  given  in  Equation  (3.45).   The  results  of 
using  F  ,  and  F 


assumed  to  be  2  •  loV-u  at  3  Hz  and  300  e~J90  at 


100  kHz. 


TABLE  4.3 


LIMITS  OF  GROUND  ADMITTANCE  FOR  THE  EXAMPLES  OF  CHAPTER 


G;    >   100k 


gl 


\l     > 


3g2 


Example  1       Example  2 


(100k) 


96k 


(96k) 


S"J    >   100k 


(100k) 


96k 


(96k) 


98k 


(98k) 


4.8k 


(4.8k) 


98k 


(98k) 


4.8k 


(4.8k) 


ohms  at  3  Hz 


ohms  at  100  kHz 


Cgl    * 
Cg2     < 


(15) 


(16) 


16 


(16) 


320 


(320) 


picofarads 


'^Values  in  parentheses  are  obtained  by  neglecting  interface 
impedance . 


Current  Control  Error 


Consider  again  the  current  control  source  shown  in  Figure  3.6, 
but  with  VQ  replaced  by  VQr   The  error  analysis  of  Chapter  3  may  be 
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applied  directly  to  the  virtual -ground  system  if  a  of  Equation  (3.60) 
is  replaced  by  B  ,  where 

BQ  -    R  (4.31) 

Q 

Using  Equations  (4.17),  (3.64),  and  (3.8),  the  expression  for  g 
becomes 


.      m31  +   (l+gy)mH 

BQ  =       R (».32) 

q 

Replacing  m31  and  mu  by  their  constituents  from  Equation  (2.22), 
Equation  (4.32)  yields 


«  ~  ZQ1  +  h31  +  hH  +  Bvhll 

PQ r — (4.33) 

q 

It  is  evident  that  6Q  is  smaller  than  a  as  given  for  the  classical 
method  in  Equation  (3.66).   If  |6q|  is  limited  to  less  than  0.01,  the 
requirement  on  R  becomes 

q 


I 

Rq   >   10°    ZQ1   +   h31   +   hH   +    Vn 


(4.34) 


which  is  to  be  compared  with  Equation  (3.69).   The  limits  on  R  for 

q 

the  examples  of  Chapter  2  are  summarized  in  Table  4.4  and  are  to  be 
compared  with  those  of  Table  3.60. 
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TABLE  4.4 

MINIMUM  VALUES  OF  CURRENT  SOURCE  RESISTANCE* 
FOR  THE  EXAMPLES  OF  CHAPTER  2 

Frequency Example  1 Example  2 

3  Hz  160M 


100  kHz         147k 


160M 

(look)        (look) 

145k 
(100k)  (98k) 


ohms 


*Based  on  ^   fr0m  Table  4.1.   Values  in  parentheses 
are  obtained  by  neglecting  interface  impedance. 

Note  that  for  Example  1,  the  virtual-ground  method  reduces  the  resistance 
requirement  to  approximately  half  that  of  the  classical  method.   For 
the  assymetrical  electrode-tissue  system  of  Example  2,  the  improve- 
ment is  negligible. 

Stability  Considerations 

If  the  virtual-ground  system  is  to  be  used  in  the  single-ended 

mode  for  broadband  measurements,  the  unity  gain  frequency  of  the 

virtual-ground  amplifier,  in  most  cases,  would  have  to  be  greater  than 

4  MHz  to  meet  common  mode  rejection  requirements  at  signal  frequencies. 

In  such  systems,  special  techniques  must  be  employed  to  assure 

stability  for  various  operating  conditions.   These  techniques  and 

l 
their  rationale  are  the  subject  of  this  section. 

The  loop  gain  of  the  virtual -ground  amplifier  circuit  shown 
in  Figure  4.3  may  be  written 


Kv  "  ~Vf)Vf)  (^.35) 
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where  1^  is  the  complex  voltage  transfer  function  from  Q2  to  P2  of  the 
network  consisting  of  the  electrode-tissue  system  and  all  of  its 
loading  admittances.   As  before,  -Ky  is  the  complex  open-loop  unloaded 
gain  of  the  virtual-ground  amplifier.   The  effect  of  amplifier  output 
inpedance  is  considered  in  a  later  paragraph.   It  is  reasonable  to 
assume  that  the  magnitude  of  K^  decreases  monotonically  with  frequency. 
If  fu  designates  the  frequency  at  which  |Ky|  becomes  unity,  the 
condition  for  stability  is 


fu  K    f*  <^36) 


where  f  is  any  frequency  at  which  the  sum  of  the  phase  lags  in  K 

v 

and  r  becomes  180  degrees. 

Nominal  magnitude  and  phase  characteristics  of  K  are  sketched 
in  Figure  4.2.   Over  most  of  the  spectrum  |  Ky  |  is  a  linear  function 
of  frequency.   The  phase  lag  is  fairly  constant  except  where  the 
slope  of  |Kj  is  changing.   At  higher  frequencies  where  stability 
becomes  a  problem,  the  phase  lag  is  at  least  90  degrees.   The 
frequency  at  which  |Kyj  crosses  unity  is  designated  by  f  .   Near 
fu,  the  phase  lag  begins  to  increase  sharply  but  does  not  pass 

through  180  degrees  until  beyond  f  . 

J  u 

Since  rv  describes  a  passive  resistance-capacitance  network, 

|r  |  is  always  less  than  unity,  so  that  f'  is  always  less  than  f 

u  u* 

It  follows  that  only  three  different  sequences  are  possible  for  the 
critical  frequencies.   With  frequency  increasing  from  left  to  right, 
these  sequences  are 
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fu>  fu>  fu  (4.37) 


fu'  fu'  f*  (4.38) 


fu>  f7r>  fu  (4.39) 


The  first  sequence  is  the  unstable  condition  to  be  avoided.   The 
second  sequence  is  stable,  but  it  places  f  in  a  region  where  the 

IT 

phase  lag  in  Ky  is  approaching  180  degrees  and  leaves  only  a  small 

margin  for  the  phase  lag  of  py.   Therefore,  it  would  not  occur,  in 

practice,  unless  f  were  relatively  low  (e.g.,  below  4  MHz).   A  low  f 

u 

proportionately  reduces  |Ky|  at  signal  frequencies,  which  is  undesirable 
The  third  sequence  represents  the  most  efficient  condition  and  is 
considered  further  in  the  section  on  stabilization  techniques. 

The  phase  lag  in  Ty  is  caused  by  YQ1 ,  Ypi,  Yp2  and  Y  which  are 
capacitive  at  high  frequencies  and  significantly  load  the  electrode- 
tissue  system.   For  a  given  set  of  loading  admittances  (of  all  four 
classes),  the  resulting  phase  lag  in  ry  depends  on  the  h  parameters 
(or  the  W  elements)  of  the  electrode-tissue  system.   Thus,  a  given 
instrumentation  system  may  operate  in  a  stable  mode  with  some 
electrode-tissue  systems,  and  not  with  others,  or  it  may  not  be  stable 
with  any. 

To  illustrate  the  foregoing  discussion,  consider  the  simplified 
equivalent  network  shown  in  Figure  4.4.   The  W  elements  correspond  to 
those  of  Figure  2.7  and  the  capacitors  represent  the  major  components 
Of   Yq-l,  Ypi,  Yp2  and  Yg,  respectively.   Symbols  v,  r  and  e  have  been 
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W, 


yvw\ 


Figure  4.4.   Simplified  Equivalent  Circuit  for  the  Virtual- 
Ground  Amplifier  Feedback  Loop 


introduced  for  convenience  in  the  next  section.   The  transverse 
loading  admittance,  the  effective  interface  impedances  and  W  have 
been  neglected.   It  is  evident  that  there  may  be  a  phase  lag  from 
Q2  to  b  and  an  additional  lag  from  b  to  v.   The  lag  from  b  to  v  will 
usually  be  the  larger  of  the  two,  but  it  may  not  exceed  90  degrees. 
The  total  lag  from  Q2  to  v  becomes  worse  (and  f  decreases)  as  W 
and  W   become  large  and  W  ,  W   and  W  become  small.   It  may  exceed 
90  degrees. 

The  output  impedance  of  the  virtual-ground  amplifier  should  be 
negligible  compared  to  the  capacitive  reactance  at  Q2  to  avoid  an 
additional  phase  lag  in  the  loop.   Output  impedances  below  50  ohms 
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are  desirable.   If  necessary,  this  condition  may  be  achieved  by  adding 
an  emitter  follower  circuit  at  the  amplifier  output  terminal. 

Stabilization  Techniques 


The  techniques  described  in  this  section  are  suggested  as  a 
means  of  achieving  the  condition  in  Equation  (4.39)  with  as  large 
an  f  as  practical.   For  clarity,  the  stabilization  process  may  be 
viewed  in  two  steps : 

Make  f  as  large  as  possible 

Select  f„  and  reduce  f '  below  f„ 
u  u        ?r 

These  steps  are  explained  in  the  following  paragraphs. 


Maximizing  f 


As  a  rule ,  f  will  lie  in  a  region  where  the  phase  lag  in  K 
is  approximately  90  degrees  (or  somewhat  larger).   Consequently,  the 
phase  lag  in  r  at  f  will  also  be  approximately  90  degrees  (or 
somewhat  smaller).   The  task  of  maximizing  f  is  equivalent  to 
minimizing  the  capacitive  effects  that  cause  the  lag  in  T    .   For 
this  purpose,  resistors  are  added  in  series  with  each  of  the  loading 
capacitors  as  indicated  by  Rqj_>  Rpj_»  Re  and  Rr  in  Figure  4.5.   For 
the  present  discussion,  neglect  the  presence  of  Rpo-   If  these 
resistors  are  sufficiently  large,  the  lag  from  Q2  to  b  and  from  b  to 
v  will  each  be  reduced  to  safe  values.   However,  upper  limits  are 
forced  on  Rq-^   and  Rp^  by  the  loading  error  conditions  which  must  be 
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»--v&vwt — vfyw 


Figure  4.5.   Stabilization  Techniques  for  the  Virtual- 
Ground  Amplifier  Feedback  Network 

satisfied  at  signal  frequencies.   Equations  (4.24)  and  (4.25)  are 
applicable  if  Cr  and  Ce  are  used  to  calculate  Ypi  and  YQ1,  and  RQ1 
and  Rpi  are  considered  as  part  of  Zq]_  and  Zpi,  respectively. 

In  the  simplified  circuit  of  Figure  4.5,  the  lag  from  Q2  to  v 
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would  peak  at  some  intermediate  frequency  (or  frequencies)  and  become 
zero  at  higher  frequencies.   Actually,  all  of  the  loading  capacitance 
cannot  be  compensated  with  series  resistance.   The  scheme  shown  in 
Figure  4.5  pertains  primarily  to  cable  shields  and  other  conducting 
objects  (table  tops,  cages,  etc.)  which  are  in  proximity  to  critical 
parts  of  the  system  and  require  a  ground  connection  to  reduce  extraneous 
noise.   These  items  usually  introduce  capacitance  in  the  range  of  20 
to  200  picofarads,  which  if  left  uncompensated  would  reduce  f  to 

TT 

undesirably  low  values  (for  example,  below  4  MHz). 

There  remain  scattered  throughout  the  Ty   network  small  parasitic 
elements  which  cause  the  lag  from  Q2  to  v  to  continue  to  increase  at 
higher  frequencies.   Included  in  this  group  is  the  input  capacitance 
of  the  virtual-ground  amplifier.   These  elements  impose  a  natural  upper 
limit  on  fr      Fortunately,  they  can  usually  be  held  to  small  values 
(typically  0.5  to  10  pf)  so  that,  in  most  cases  of  interest,  ff  would 
be  greater  than  10  MHz. 


Reducing  f1 
a — u 


As  indicated  earlier,  f^  will  always  be  less  than  f  .   Thus 
if  fu  were  not  much  greater  than  f f ,  the  existing  attenuation  of  Tv 
might  be  sufficient  to  place  f<  below  f ^ ,  resulting  in  a  stable 
condition.   This  possibility  unnecessarily  restricts  f  and  is  not 
reliable  since  f%   and  ry  will  vary  with  different  electrode-tissue 
systems.   In  some  cases,  such  a  system  might  actually  be  operating  in 
the  condition  of  Equation  (4.38).   Furthermore,  as  f  is  decreased 
toward  f^,  f^  also  decreases  because  of  the  increasing  phase  lag 
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in  K  .   It  is  evident  that  a  larger  f  may  be  used  if  fr  is  reduced 
v  u  u 

by  additional  means. 

To  reduce  f^,  it  is  convenient  to  use  the  lag  attenuator  network 
consisting  of  Cv,  Ry,  Wp2  and  Rp2  as  indicated  in  Figure  4.5. 
Normally,  Rp2  and  Ry  are  the  controllable  elements,  while  C   is  a 
cable  shield  capacitance  and  Wp2  is  an  electrode-tissue  system 
parameter.   Typical  magnitude  and  phase  characteristics  of  the 
transfer  function  from  b  to  v  of  this  network  are  sketched  in  Figure 
4.6. 

The  break  frequencies  of  the  lag  network  are  given  by 


1   2TT(Wp2  t  Rp2  +  Rv)Cv 


(4.40) 


f„  =  i_ 

2    2*RVCV 


(4.41) 


and  the  attenuation  factor  is 


(4.42) 


The  peak  phase  lag  occurs  at 


f  =  /  f  f 
P     "12 


(4.43) 


and  is  given  by 
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Figure  4.6.   Typical  Magnitude  and  Phase  Characteristics 
of  a  Lag  Network 


8  s  (90  -  2tan_1/T")  degrees 


(4.44) 


which  may  be  rearranged  to  give 


5  =  [tan  (45  -  ^P)]2 


(4.45) 
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The  major  design  parameters  are  ^  and  ?.   The  selection  of  these 
factors  is  not  unique.   It  is,  however,  subject  to  certain  con- 
straints which  are  discussed  in  the  following  paragraphs. 

Restrictions  on  ^.--At  signal  frequencies,  loading  error  in 
the  P2  arm  must  be  negligible.  The  condition  in  .Equation  (4.26)  is 
applicable  and  may  be  rearranged  to  give 

2*|zp2  ?  h13  ;  Rp2|C;  >40fleJ  (*.«) 

where  Rp2  has  been  included  as  part  of  the  electrode-tissue  system, 
Yp2  has  been  replaced  by  the  admittance  of  C  ,  and  R  has  been 
neglected.   The  worst  case  occurs  at  the  highest  signal  frequency  of 
interest,  hereafter  designated  by  fg.   At  ffi  the  left  term  in  Equation 
(4.46)  will  usually  approximate  ^  as  given  in  Equation  (4.40).   Thus 
to  a  first  approximation,  f±   should  meet  the  following  condition 

fl>40fsl6vl  (4.47) 

where  ^   is  evaluated  at  fg.   This  condition  holds  whether  the  system 
is  operated  in  the  single-ended  or  the  differential  mode.   In  the 
single-ended  case,  |PJ  must  be  very  small  so  that  f  could,  in 
principle,  be  less  than  f g .   However,  a  factor  of  safety  may  be 
necessary  because  of  the  variable  nature  of  H     Even  in  the 
differential  mode,  | gy |  would  normally  be  less  than  unity  so  that  the 
limit  in  Equation  (4.47)  would  never  be  greater  than  40f 
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Restrictions  on  g.— The  peaked  phase  lag  characteristic  of 

Figure  4.6  raises  the  possibility  of  a  new  f  in  the  region  between 

7r  to 

f±   and  f2  in  addition  to  the  inherent  f  at  higher  frequencies.   To 

avoid  this  potentially  unstable  situation,  the  phase  lag  from  Q2  to 

v  at  frequencies  between  f  and  f2  must  be  less  than  90  degrees.   This 

lag  is  likely  to  be  greatest  at  approximately  f  .   If  the  phase  lag 

from  Q2  to  b  at  f  is  designated  as  4>   then  0  must  satisfy  the 
P  P        p  J 

following  inequality 


3p  <  (90  -  <j>  )  degrees  (4.48) 


An  equivalent  inequality  on  £  results  if  Equation  (4.48)  is  substituted 
into  Equation  (4.45) 


c  > 


tan  IE- 


(4.49) 


Numerical  values  for  this  inequality  are  given  in  Table  4.5, 


TABLE  4.5 
LIMITS  OF  £  FOR  GIVEN  VALUES  OF 


i     degrees 
P 

£> 

10 

.0075 

20 

.031 

30 

.072 

40 

.13 

50 

.22 

60 

.33 
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Since  f  must  be  on  the  same  order  as  or  greater  than  f  ,  f  will 
x  s   p 

usually  fall  in  a  region  (e.g.,  above  1  MHZ)  where  0  may  not  be 
neglected.   To  be  sure,  at  least  a  10-degree  factor  of  safety  is 
necessary.   In  many  cases,  much  larger  values  may  be  required.   It 
follows  that  5  cannot,  in  general,  be  made  arbitrarily  small. 

Amplifier  considerations  .  —  Since  the  minimum  value  of  £  is 
restricted,  fu  cannot  be  arbitrarily  large  if  f^  is  to  be  kept  below 
the  inherent  limit  on  f  ^ .   In  a  properly  adjusted  attenuator  network, 
fu  wil1  fal1  above  f2  and  wil1  be  somewhat  lower  than  gf  due  to  the 
additional  attenuation  from  Q2  to  b.   The  input  capacitance  of  the 
virtual-ground  amplifier  should  be  as  small  as  possible  (less  than 
10  pf)  to  minimize  the  loading  of  Ry  and  maximize  the  inherent  limit 
of  f^.   To  accommodate  a  wide  range  of  electrode-system  parameters, 
it  is  desirable  to  make  Rv  and  possibly  Rp2  variable. 

In  summary,  stabilization  of  the  virtual-ground  system  involves 
the  selection  of  at  least  four  resistors  (two  for  the  P2  arm,  and 
one  each  for  the  PI  and  Ql  arms).   Some  resistance  can  often  be  put 
in  the  ground  paths  of  Cr  and  Cg  with  similar  stabilizing  effects. 
The  resistor  selection  is  somewhat  arbitrary  and  will  depend  largely  on 
how  a  particular  system  is  implemented.   The  most  satisfactory  set  of 
values  generally  must  be  found  by  an  educated  trial  and  error 
procedure  based  on  the  theory  presented  and  certain  common  considera- 
tions. 

The  stabilizing  resistors  which  occur  in  series  with  the 
terminals  of  the  electrode-tissue  system  are  indicated  in  Figure  4.7. 
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Figure  4.7.   Electrode-Tissue  System  with  Stabilizing 
Resistors 


It  is  convenient  to  view  these  resistors  as  additions  to  the 
corresponding  effective  interface  impedances,  i.e.,  Z   ,  Z   and  Z 
With  this  simple  modification,  the  analytical  results  derived  in 
previous  sections  may  be  applied  directly  to  the  compensated  system. 
The  same  convenience  applies  to  the  noise  analyses  that  follow. 

Noise  in  the  Virtual-Ground  System 
Instrumentation  Noise 


Instrumentation  noise  arises  in  the  virtual-ground  amplifier 
and  the  response  amplifier.  For  reasons  similar  to  those  discussed 
in  Chapter  3,  the  instrumentation  noise  of  greatest  significance  is 
that  which  may  be  described  by  means  of  a  noise  voltage  source  in 
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series  with  the  amplifier  input  terminals.   These  sources  are  represented 
as  N^  and  Np  in  Figure  4.3.   N  produces  a  noise  component  at  the 
input  to  the  response  amplifier  which  together  with  N  comprises  the 
instrumentation  noise  that  is  amplified  along  with  the  desired 
response  signal. 

The  noise  due  to  Nv  at  the  input  of  the  response  amplifier  may 
well  be  greater  than  Ny  itself.   To  verify  this  statement,  consider 
again  the  simplified  network  in  Figure  4.5.   The  gain  for  N  from 
the  input  of  the  virtual-ground  amplifier  to  Q2  is  given  by 


-K 

v 


*n  =  i  ;  k  r  ■  (4-50) 


v  v 


where,  as  before,  r   is  the  transfer  function  from  02  to  v  and  -K 

v  v 

is  the  open-loop  unloaded  gain  of  the  amplifier.   If  r  designates 

r     & 

the  complex  transfer  function  from  Q2  to  r  (the  input  to  the  response 
amplifier),  then  the  power  spectral  density  of  the  noise  at  r  due  to 
N  is  given  by 

Nrv(f)  =  IVf)rp(f)l%(f)  (4-51> 

where  the  symbol  N  is  taken  to  mean  power  spectral  density.   If 
I KNrr I  is  greater  than  unity  in  a  certain  frequency  band,  then  the 

noise  power  at  r  will  be  greater  than  the  noise  power  of  N  in  that 

v 

band. 

As  a  rule,  the  instrumentation  system  is  designed  to  minimize 
loading  at  the  PI  terminal  so  that  |r  |  will  usually  remain  close  to 
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unity  at  all  frequencies  of  concern.   Therefore,  it  follows  from 
Equation  (4.51)  that  Npv  will  be  greater  than  Ny  at  approximately 
those  frequencies  where  |kn|  exceeds  unity.   In  the  signal  band  (below 
200  kHz)  and  perhaps  somewhat  higher,  | ^ |  is  approximately  equal  to 
unity  (see  Equation  [4.50])  because  |kJ  is  large  and  |r  |  is 
approximately  unity.   At  higher  frequencies,  |iyv|  crosses  unity  at 
f£,  while  |Kv|  crosses  unity  at  fu<   As  indicated  in  the  section  on 
stability,  f£  is  usually  well  below  fu<   Hence  |kn|  may  be  greater  than 
unity  for  frequencies  above  the  signal  band  and  up  to  approximately  f 
(e.g.,  at  frequencies  between  300kHz  and  30MHz).   In  at  least  part  of 
this  range,  the  instrumentation  noise  at  PI  due  to  N  will  generally 
be  enhanced.   Experimental  evidence  indicates  that  this  enhancement 
will  usually  correspond  to  less  than  a  factor  of  two  in  the  noise 
voltage  level  at  PI. 

It  is  evident  that  the  virtual-ground  system  is  inherently 
noisier  than  the  classical  system  with  respect  to  a  given  level  of 
amplifier  noise.   Fortunately,  a  significant  portion  of  the  power  in 
this  noise  lies  at  frequencies  well  above  the  usual  impedance  measure- 
ment band  and  it  can  be  filtered  out  after  amplification.   In  addition, 
because  of  its  regular  random  nature,  instrumentation  noise  may  be 
eliminated  by  means  of  standard  averaging  techniques. 

Environmental  Noise 

The  effects  of  environmental  noise  on  the  virtual-ground 
system  are  similar  to  those  found  in  the  classical  system.   Consider, 
for  example,  the  electrostatically  coupled  noise  currents   I 
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indicated  in  Figure  If. 3.   The  noise  voltage  at  PI  due  to  these  currents 


is 


V, 


/    =  v   +   2x 
PI*    Px   1  +  K 


(4.52) 


v 


This  expression  is  the  dual  of  Equation  (3.73)  with  Vpix  replacing 
Vpxpand  1  +  Ky  replacing  p.   The  values  of  Vpx  and  V^  are  essentially 
as  given  by  Equations  (3.74)  and  (3.76),  respectively.   Finally,  the 
virtual-ground  equivalent  of  Equation  (3.77)  is 


'Plx 


Ql 


ZP1  "  ZP2  +  h12 


'13 


11 


Zp2  +  2ZQ2  +  h13  +  2h23  +  2h33" 
(1  +  K^h" 


(4.53) 


As  before,  the  cancellation  occurring  in  the  numerator  of  the  left 

term  in  the  brackets  is  neither  complete  nor  easy  to  control  and  cannot 

be  used  as  a  reliable  index  of  noise  reduction.   Also,  the  term  on 

the  right  in  the  brackets  is  not  likely  to  be  important  because  K 

v 
will  usually  be  quite  large. 

If  the  signal  voltage  is  required  to  be  greater  than  ten  times 
the  noise  voltage,  then  the  dual  of  Equation  (3.78)  is  for  the  virtual- 
ground  system 
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*»'.    ,  10 


Jx 


SP1   ~   ZP2   +  h12  "  h13 
hll 


+     ZP2  +   2ZQ2  +  h13  +   2h23  +   2h33 
(1  +   Kv)hll 


(4.54) 


As  in  the  classical  system,  this  inquality  can  seldom  be  satisfied 
unless  the  potential  cables  are  shielded. 

■ 
Electrode-Tissue  Noise 

The  electrode-tissue  noise  at  the  input  of  the  response  amplifier 

in  Figure  4.3  is  equal  to  the  noise  in  Vp  (including  two  interfaces), 

plus  a  contribution  from  the  noise  in  V.  due  to  the  fact  that  K  is 

^  v 

not  infinite.   This  is  similar  to  what  occurs  in  the  classical  system 
except  that  Kv  takes  the  place  of  the  amplifier  common  mode  rejection 
ratio.   As  a  rule,  electrode-tissue  noise  in  V2  win  not  saturate  the 
virtual-ground  amplifier  because  it  appears  in  the  feedback  path 
where  the  dynamic  range  is  large  (e.g.,  ±10  volts).   On  the  other  hand, 
noise  in  Vp  may  saturate  the  response  amplifier  and/or  the  instruments 
that  follow  it.   This  problem  is  resolved  conveniently  by  means  of  the 
negative  feedback  filter  loop,  described  in  the  next  section. 

The  Negative-Feedback  Filter  Loop 

The  negative- feedback  filter  loop  provides  a  high  performance 
alternative  to  using  a.c.  coupling  between  the  response  amplifier  and 
the  electrode-tissue  system.   It  permits  the  filter  for  d.c.  and 
noise  rejection  to  be  placed  outside  the  main  signal  channel,  where  it 
introduces  no  spurious  effects  on  the  signal  and  adjustments  are 
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particularly  convenient  to  make.   This  arrangement  is  especially 
desirable  in  a  system  which  is  to  operate  over  a  broad  frequency 
spectrum.   The  feedback  filter  loop  is  also  convenient  for  reducing 
the  operating  point  disturbances  that  accompany  changes  from  one 
electrode-tissue  system  to  another. 

Many  operational  amplifiers  are  available  with  a  noninverting 
input  channel  in  addition  to  the  standard  inverting  connection.   In 
the  virtual-ground  system  of  Figure  4.8,  the  output  of  the  response 


Figure  4.8.   Virtual-Ground  System  with  Negative-Feedback 
Filter  Loop 
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amplifier  is  fed  back  through  a  low-pass  filter  to  the  ncninverting 
side  of  the  virtual-ground  amplifier.   Signals  with  frequencies  well 
within  the  filter  pass  band  are  subject  to  heavy  negative  feedback 
while  those  with  frequencies  well  beyond  the  pass  band  are  not  fedback 
at  all.   In  this  way,  the  system  is  made  to  fully  amplify  the  higher 
frequency  signals,  while  the  lower  ones  are  passed  with  little  or  no 
amplification.   The  feedback  filter  loop  is  relatively  insensitive 
to  extraneous  disturbances  and  greatly  facilitates  the  use  of  variable 
and  large  filter  components. 

With  reference  to  Figure  4.8,  the  output  voltage  of  the  response 
amplifier  is 

Vr  =  "Kr(VP  +  VP2}  (4.55) 

where  -Kp  is  the  complex  gain  of  the  response  amplifier.   Further,  it 
is  not  difficult  to  show  that 


vp2 _____ (1+.55) 

v 


where  p  is  the  common  mode  rejection  ratio  of  the  virtual-ground 

amplifier  and  L(f )  is  the  complex  voltage  transfer  function  of  the 

feedback  filter.   As  before,  Kv  is  the  unloaded  open-loop  complex 

gain  of  the  virtual-ground  amplifier  and  the  output  impedance  is 

neglected.   Substituting  Equation  (4.56)  into  Equation  (4.55)  and 

solving  for  V  yields 
r  " 
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vr  *  rrtr  (VP  +  r  }  <*.«-) 


where  the  1/p  in  the  numerator  and  the  unity  in  the  denominator  of  the 
right  side  of  Equation  (4.56)  have  been  neglected.   The  expression 
for  Vr  is  discussed  in  the  following  paragraphs  for  desirable  and 
undesirable  signals. 

- 

Signals  due  to  IQ1.— For  signals  produced  by  IQ1,  the  term  in 
parenthesis  in  Equation  (4.57)  is  equal  to  Vpj_  when  there  is  no  feedback 
filter  and  is  given  by  Equation  (4.1),  where  3y  is  as  given  in  Equation 
(4.6).   Excitation  for  impedance  measurements  is  normally  produced  at 
frequencies  well  above  the  filter  pass  band  where  LiC«l.   For  this 
case,  Equation  (4.57)  becomes 

V 

which  is  as  though  the  filter  were  absent. 

Electrode-tissue  noise.— The  electrode-tissue  noise  voltages  of 
interest  occur  in  Vp  and  V2  and  are  hereafter  designated  by  Vp  and 
V2e,  respectively.   Looking  out  of  the  terminals  of  the  electrode- 
tissue  system,  the  effective  impedance  presented  by  the  external 
circuit  is  quite  high.   Therefore,  Vpe  and  V2e  correspond  to  what 
would  be  observed  if  the  electrode-tissue  system  were  not  connected. 
For  noise  voltages  falling  well  within  the  filter  pass  band,  LK  »1 
and  Equation  (4.57)  yields  for  the  noise  output 


mi 


Vre  ~-   "  T(VPe  +  ^  C*.59) 

v 

In  a  typical  passive  filter,  L  is  unity  from  d.c.  through  an  arbitrary 
band  of  frequencies.   In  this  range,  the  noise  is  passed  without 
amplification.   At  higher  frequencies,  L  begins  to  decrease  and  the 
noise  amplification  increases,  eventually  catching  up  with  that  of  the 
signal  (Equation  [4.58]). 

To  illustrate  the  advantages  of  the  filter  loop,  let  K  be  1000 
and  assume  that  the  output  stage  of  the  amplifier  has  a  maximum  range  of 
±10  volts.  If  this  system  were  not  compensated,  it  would  tolerate,  at 
most,  ±10  millivolts  of  d.c.  and  low  frequency  electrode-tissue  noise. 
With  the  feedback  filter,  the  system  could  handle  up  to  ±10  volts  of  the 
noise.  In  addition,  "the  signal-to-noise  ratio  is  improved  by  a  factor 
of  1000. 

Instrumentation  noise. --Most  of  the  energy  in  instrumentation 
noise  lies  at  frequencies  well  above  the  pass  band  of  the  feedback 
filter.   Accordingly,  the  analysis  of  instrumentation  noise  given  in 
a  previous  section  applies  directly  in  the  presence  of  the  filter. 

Environmental  noise.— For  environmental  noise  due  to  equal 
currents,  Ix,  the  term  in  parenthesis  in  Equation  (4.57),  is  equal 
to  Vpix  as  given  in  Equations  (4.52)  and  (4.53).   For  noise  frequencies 
well  within  the  filter  pass  band, 


Vrx  =  "  TVPlx  (1-60) 
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It  is  apparent  that  the  feedback  filter  loop  may  be  used  for  supple- 
mentary reduction  of  environmental  noise  provided  that  the  noise 
frequencies  are  well  below  the  signal  frequencies  of  interest. 

Examples  of  filter  loop  operation 

To  illustrate  the  operation  of  the  feedback  filter  loop,  let 
the  feedback  filter  be  as  shown  in  Figure  4.9. 


AA/W 


Figure  4.9.   Simple  Low-Pass  Filt 


er 


For  this  filter, 


1  +  jf/f. 


(4.61) 


and  the  overall  transfer  function  indicated  in  Equation  (4.57)  b 


ecomes 


-K 


-K     1  +  jf/f 


L   1  +  LK    1  +  K    I  +  jf/Q  ;   K  )f 
r        r  r   L 


(4.62) 
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where 


f .  =  i_ 

L    2tiR  C 


(4.63) 


L  L 


If,  instead  of  using  the  low-pass  feedback  filter,  the  high -pass  filt 
shown  in  Figure  4.10  were  placed  in  the  input  path  of  the  response 
amplifier,  the  overall  transfer  function  would  be 


er 


where 


jf/ft 


r  1  +  jf/f 


(4.64) 


f..  = 


H   2*RHCH 


(4.65) 


. 


Figure  4.10.   Simple  High-Pass  Filter 
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The  asymptotic  log  modulus  and  phase  characteristics  of  TL  and 

TH  are  sketched  in  Figure  4.11  for  the  values  (1  +  K  )fr  =  f „  =  1Hz 

r   L    H 

and  Kp  =  1000.   The  characteristics  of  T  and  Tt  are  essentially 
identical  except  at  frequencies  near  10~3Hz  and  below.   In  this  region, 
the  transfer  function  of  the  feedback  method  remains  at  unity,  while 
the  forward-path  filter  continues  to  attenuate.   Since  the  gain  is 
rather  high,  the  advantage  of  extra  attenuation  in  the  latter  case  is 
not  especially  significant.   Note  that  to  achieve  the  same  upper  break 
frequency,  fL  must  be  about  one  thousandth  of  f  with  correspondingly 
larger  circuit  elements.   However,  since  the  feedback  filter  occurs 
outside  the  amplifier  and  one  end  of  the  capacitor  is  grounded,  the 
larger  components  present  no  special  problem. 

A  disadvantage  of  simple  filters  is  that  the  desired  signal 
frequencies  must  be  well  above  the  upper  break  frequency  to  avoid 
phase  distortion.   For  example,  with  reference  to  TL  and  TH,  if  phase 
distortion  is  not  to  exceed  2  degrees,  signal  frequencies  must  be  at 
least  29  times  greater  than  the  break  frequency. 

The  "dead  space"  between  significant  noise  attenuation  and 
the  lowest  allowable  signal  frequency  may  be  reduced  by  adding  a  high 
pass  filter  in  the  forward  path  outside  the  feedback  loop  in  addition 
to  the  feedback  filter.   For  example,  with  the  simple  filters  of 
Figures  4.9  and  4.10,  the 'overall  transfer  function  would  be 


T   -     '**  '  *   5f/fL  Jf/fH  ,     , 

LH   1  +  K    1  +  jf/(i  +  K  )f    1  +  jf/f  (U'56) 

L  r   L  H 
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Figure  4.11.   Overall  Transfer  Function  Characteristics  Using 
Several  Methods  of  Noise  Filtering 
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If  the  denominator  break  frequencies  are  made  equal,  i.e.,  (1  +  K  )f  = 

fH,  then  for  less  than  2  degrees' phase  distortion,  signal  frequencies 

must  be  58  fy  or  greater.   It  follows  that  for  the  same  lowest  allowable 

signal  frequency,  the  break  frequencies  of  the  combined  filters  must 

be  half  as  large  as  the  break  frequency  of  either  filter  acting  alone. 

The  phase  and  log  modulus  characteristics  of  TT „  are  indicated  in 

LH 

Figure  4.11  for  the  case  (1  t  Kp)fL  =  fR  =  0.5  Hz  and  K   =  1000.   It 
is  evident  that  the  noise  attenuation  is  significantly  improved  with 
the  combined  filters  for  the  same  lowest  allowable  signal  frequency. 
Further  reduction  cf  the  filter  dead  space  between  noise  and 
useful  signal  frequencies  may  be  achieved  with  higher  order  and 
nonlinear  filters  in  the  feedback  and  forward  paths.   In  the  design 
of  second  order  and  higher  filters  for  the  feedback  loop,  care  must 
be  exercised  to  prevent  the  appearance  of  underdamped  complex  roots 
in  the  closed-loop  transfer  function.   Such  roots  may  cause  objec- 
tionable  ringing  to  occur  in  the  presence  of  common  disturbances. 


CHAPTER  5 


AN  EXPERIMENTAL  SYSTEM  BASED  ON  THE  VIRTUAL-GROUND 

CONFIGURATION 


Introduction 

The  theory  and  analysis  of  the  virtual -ground  system  were  given 
in  Chapter  4.   Much  of  that  work  was  based  on  experience  obtained  with 
the  experimental  system  described  in  the  present  chapter.   The  system 
was  first  assembled  by  the  author  in  March  of  1967,  and  minor  improve- 
ments have  been  made  as  needed.   Since  that  time  it  has  been  used 
regularly  to  measure  impedance  in  certain  animal  tissues  in  situ.   The 
major  application  has  been  concerned  with  the  impedance  properties  of 
the  cerebral  cortex  in  unrestrained  rats.   This  tissue  is  especially 
difficult  to  measure  accurately  because  electrode  dimensions  must  be 
quite  small.   Hence,  it  provides  a  proper  proving  ground  for  the 
virtual -ground  technique.   Some  of  the  results  of  tissue  measurements 
will  be  reported  later  in  the  chapter. 

The  basic  configuration  of  the  experimental  system  is  shown  in 
Figure  4.3  and  may  be  referred  to  for  the  discussions  that  follow.   The 
system  consists  of  three  major  components,  the  virtual-ground  amplifier, 
the  response  amplifier  and  the  low-pass  filter.   In  addition,  a  special 
current  control  circuit  has  been  designed  to  supply  the  signal  current. 
The  various  circuits  are  discussed  in  detail  following -a  description 
of  the  physical  layout. 
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The  Physical  Layout 

A  close-up  photograph  of  the  experimental  virtual-ground  system, 
with  important  components  labeled,  is  shown  in  Figure  5.1.   Four  solid 
state  operational  amplifier  units,  similar  to  the  one  indicated,  are 
plugged  into  sockets  on  the  circuit  board.   The  second  response 
amplifier  costs  approximately  $85.00.   The  first  response  amplifier, 
the  virtual-ground  amplifier,  and  the  current  control  amplifier  are 
identical  types  and  cost  about  $155.00  each.   The  remaining  components 
of  the  system,  less  power  supply,  total  less  than  $15.00.   The  entire 
system  may  be  powered  from  a  single  supply  with  plus  and  minus  15 
volt  d.c.  outputs  and  a  30  ma  current  capacity.   The  power  unit  shown 
in  the  picture  is  much  larger  than  necessary. 

Connection  to  the  Electrode-Tissue  System 

Another  viev;  of  the  experimental  virtual-ground  system, 
including  additional  paraphernalia,  is  shown  in  Figure  5.2.   The 
electrode-tissue  system  (in  this  case,  rat  plus  electrode  assembly) 
is  connected  to  the  circuit  board  by  means  of  four  miniature  coaxial 
cables.   The  voltage  and  current  cables  are  twisted  in  pairs  to 
minimize  extraneous  electrical  signals  and  coupling  between  the  current 
and  voltage  circuits.   The  cables  extend  about  2.5  feet  from  the 
circuit  board  and  end  in  a  small  female  connector.   The  mating  male 
connector  is  attached  to  the  electrode-tissue  system  through  approxi- 
mately seven  inches  of  unshielded  wire.   This  arrangement  is  light 
and  flexible  and  permits  a  great  deal  of  versatility  in  connecting  to 
various  impedance  subjects.   The  female  connector  also  serves  as  a 


149 


e 

Q) 
■P 

W 

>> 
CO 

T> 
G 
3 
O 

O 
l 

>H 
ffl 

+■> 
b 

> 


IT) 
4-> 
C 
<D 
E 

•H 
ft 
CD 

W 

CD 

ja 

■P 

in 

o 

a 


o 

O 

a. 
D 

CD 
CO 

o 

H 
O 


in 


150 


+-> 

B 

§ 

bO 
C 

03 

i 

< 

■0 

■P 
C 
0) 

e 

•h 

h 

S1 

W 

a; 

-M 

4-1 
O 

s 

•H 
> 

4-» 
<D 
H 

I 


bO 
Em 


151 


convenient  location  for  the  stabilization  resistors  Rp5,  Rpi»  ar>d  Rqi 
as  discussed  in  Chapter  4.   In  the  present  system,  2.2  kilohms  has  been 
a  satisfactory  value  for  Rp„  and  Rp, ,  and  no  resistor  has  been  necessary 
for  RQ1. 

The  female  connector  actually  comprises  two  sets  of  receptacles 
which  have  been  soldered  together.   This  feature  facilitates  changes 
from  one  electrode-tissue  system  to  another  (one  of  which  may  be  for 
reference)  and  helps  to  minimize  amplifier  operating  point  disturbances 
during  the  changes.   The  double-ended  connector  also  aids  in  the 
realization  of  an  important  advantage  of  the  virtual-ground  technique , 
i.e.,  the  ability  to  check  whether  the  common-mode  error  signal  at  P2 
is  indeed  negligible.   Hence,  during  an  experiment,  one  simply 
disconnects  the  unshielded  PI  lead  at  the  electrode-tissue  system  and 
plugs  it  into  the  vacant  P2  receptacle  on  the  female  connector,  leaving 
the  other  leads  undisturbed.   The  system  output  then  reflects  the 
voltage  at  P2,  which  should  be  negligible  compared  to  the  response  at 
PI.   Similarly,  this  technique  may  also  be  used  to  observe  the  voltages 
at  Ql  and  Q2  by  plugging  the  PI  lead  into  the  appropriate  vacant 
receptacle.   The  connection  at  the  electrode-tissue  system  itself  is 
discussed  later  in  the  chapter. 

As  explained  in  Chapter  4,  the  shields  of  the  PI  and  P2  cables 
are  especially  important  to  reduce  environmental  noise.   The  shield  of 
the  Ql  cable  is  usually  not  important  and  that  of  the  Q2  cable  may  be 
completely  omitted. 

The  total  shield  to  center  conductor  capacitance  of  each  coaxial 
cable  is  approximately  130  picofarads.   To  avoid  the  effects  of  this 
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capacitance,  each  shield  is  connected  in  a  special  way  on  the  circuit 
board.   The  shields  of  the  Ql  and  PI  cables  are  maintained  at  the  same 
voltage  as  their  center  conductors  to  minimize  the  effective  loading 
admittance  at  signal  frequencies.   These  connections  are  described 
in  the  circuit  section.   The  shield  of  the  P2  cable  is  connected  to 
ground  by  means  of  a  miniature  500  ohm  potentiometer  which  forms  R 
of  the  stabilization  scheme  discussed  in  Chapter  U.   The  shield  of  the 
Q2  cable  is  generally  not  used. 

Ground  Connection 

The  system  is  connected  to  ground  by  means  of  a  braided  ribbon 
conductor  (the  ground  strap)  which  is  approximately  0.5  inch  wide  and 
3  feet  long.   The  ground  strap  serves  as  the  signal  common  and  is 
connected  to  an  arbitrary  ground  point  on  the  peripheral  equipment. 
No  other  ground  path  is  provided  between  the  circuit  board  and  the 
equipment  (house)  ground.   The  power  supply  common  and  the  metal  case 
of  each  amplifier  are  connected  to  ground  on  the  circuit  board. 


Output  Connection 

The  amplified  voltage  response  is  led  out  of  the  circuit  board 
through  approximately  three  feet  of  RG  58A/U  coaxial  cable  (the  output 
cable)  which  terminates  in  a  BNC  connector.   The  shield  is  connected 
to  ground  through  a  270  ohm  resistor  on  the  circuit  board.   This 
resistor  is  a  precaution  to  insure  the  stability  of  the  output  section 
of  the  response  amplifier  under  the  load  of  the  shield  capacitance. 
The  shield  is  necessary  to  minimize  feedback  from  the  output  cable  to 
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low  level  parts  of  the  system,  which  might  otherwise  cause  instability. 

The  center  conductor  (but  not  the  shield)  of  the  output  cable 
is  connected  to  an  oscilloscope  which  is  used  to  monitor  the  voltage 
response.  From  this  point,  another  coaxial  cable  is  run  to  the  input 
of  the  feedback  filter  and  other  measuring  apparatus.  The  shield  of 
this  cable  is  joined  to  that  of  the  output  cable.  The  output  of  the 
feedback  filter  returns  to  the  circuit  board  by  means  of  a  miniature 
coaxial  cable.  The  cable  shield  is  grounded  at  the  circuit  board  and 
serves  as  the  ground  connection  for  the  filter. 

Input  Connection 

The  current  control  circuit  is  excited  by  means  of  an  external 
voltage  generator.   The  cable  from  the  generator  usually  runs  directly 
to  an  oscilloscope  which  monitors  the  input  signal.   At  this  point,  it 
feeds  the  input  cable  of  the  virtual-ground  system  which  consists  of 
2.5  feet  of  RG  58A/U  coaxial  cable  terminated  in  a  BNC  connector.   The 
shields  of  both  cables  are  grounded  to  the  equipment. 

Circuit  Details 
The  Response  Amplifier 

The  response  amplifier  consists  of  two  sections  as  shown  in 
Figure  5.3.   The  first  section  is  a  common  voltage  follower  configura- 
tion with  a  nominal  noninverting  gain  of  3.13.   A  voltage  divider  is 
connected  at  its  output  to  provide  a  unity  gain  point  for  the  shield 
of  the  PI  cable.   The  input  impedance  of  this  section  corresponds  to 
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Figure  5.3.   The  Response  Amplifier  Circuit 

10   ohms  in  parallel  with  5  picofarads  which  is  more  than  adequate 
for  virtual-ground  operation.   These  values  stem  from  a  field  effect 
transistor  input  stage  which  also  requires  negligible  d.c.  offset 
current  (less  than  150  picoamperes)  from  the  electrode-tissue  system. 
The  second  section  of  the  response  amplifier  is  a  common  inverter 
circuit  with  a  nominal  gain  of  100. 

An  analysis  of  the  follower  and  inverter  configurations  is 
presented  in  Appendices  2  and  3.   The  first  and  second  amplifier  units 
have  unity  gain  frequencies  of  at  least  30  MHz  and  100  MHz,  respectively, 
Hence,  according  to  Equation  A2.1,  the  first  section  should  have  less 
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than  1  percent  decrease  in  gain  out  to  100  kHz.   The  gain  of  the  second 
section  tends  to  decrease  somewhat  as  100  kHz  is  approached  due  to  the 
finite  gain  of  the  amplifier  unit  (see  Equation  [A3.1]).   This  feature 
may  be  used  to  advantage  to  offset  a  slight  current  increase  which 
occurs  due  to  parasitic  capacitance  in  the  current  control  circuit  at 
high  frequencies  (see  Current  Control  Error  in  Chapter  3).   Additional 
flattening  of  the  system  response  (from  input  signal  voltage  to  system 
output  for  a  purely  resistive  electrode-tissue  system)  has  been 
accomplished  by  adding  an  adjustable  series  resistor-capacitor  combina- 
tion in  parallel  with  the  47  kilohm  feedback  resistor. 

The  high  gain  of  the  response  amplifier  predisposes  the  system 
to  instability  due  to  spurious  feedback  from  the  system  output  to  the 
PI  cable  (center  conductor  and  shield)  and  less  importantly  to  the  Ql 
and  P2  cables.   Therefore,  a  shielded  cable  is  used  for  the  output 
channel.   The  shield  is  connected  to  ground  through  a  270  ohm  resistor 
to  insure  the  stability  of  the  output  section  against  the  reactive 
load  of  the  shield  capacitance.   Another  precaution  is  to  connect  a 
suitable  capacitor  from  the  shield  of  the  PI  cable  to  ground.   In  the 
present  system,  a  300  picofarad  capacitor  has  been  satisfactory.   Acting 
in  parallel  with  the  470  ohm  divider  resistor,  this  capacitor  does  not 
affect  the  +1  gain  return  to  the  PI  shield  at  signal  frequencies,  but 
becomes  effective  at  higher  frequencies  where  instability  normally 
occurs.   It  also  serves  to  counter  the  increase  in  high  frequency 
amplifier  noise  which  accompanies  the  +1  connection. 

The  Virtual-Ground  Amplifier 
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The  virtual-ground  amplifier  was  chosen  for  its  low  d.c.  offset 
current,  high  gain  and  broad  bandwidth,  and  is  identical  in  design  to 
the  first  response  amplifier.   This  particular  unit  exhibits  a  frequency 
dependent  output  impedance  which  decreases  from  roughly  50  kilohms  at 
very  low  frequencies  to  less  than  50  ohms  at  100  kHz. 

Recall  that  the  common  mode  rejection  of  the  system  is  described 
by  8v  as  defined  in  Equation  (4.2).   If  | Bv|  is  to  be  less  than  0.01, 
then  according  to  Equation  (4.5),  the  voltages  of  the  electrode-tissue 
system  under  study  must  satisfy  the  following  requirement 


(5.1) 


where  the  ZQ  and  the  unity  terms  have  been  neglected.   For  the  amplifier 
used  here,  the  right  side  of  Equation  (5.1)  becomes  approximately  2,000 
at  very  low  frequencies  and  3  at  100  kHz.   That  these  conditions  are 
met  in  a  particular  experiment  may  be  easily  checked  as  described  in  a 
previous  section. 

The  Feedback  Filter 

The  feedback  filter  is  of  the  form  shown  in  Figure  4.9.   The 
capacitor  is  a  large  fixed  element  of  15  uf  and  is  nonelectrolytic . 
Resistor  values  from  100  megohms  to  100  kilohms  haVe  been  provided  for 
selection  into  the  filter  circuit.   Since  the  response  amplifier  gain 
is  313,  these  values  correspond  to  upper  break  frequencies  of  0.031  Hz 
to  31  Hz  in  the  system  transfer  function  (see  Equation  [4.63]).   At  the 
higher  resistance  setting,  the  system  can  handle  signal  frequencies 


157 


down  to  0.9  Hz  with  less  than  2  degrees' phase  distortion.   At  the  lower 
resistance  setting,  the  signal  frequencies  would  have  to  be  greater 
than  900  Hz,  but  all  noise  below  31  Hz  would  be  attenuated.   VJider 
ranges  of  operation  can  easily  be  achieved  with  other  resistor  and 
capacitor  values. 

It  is  especially  significant  to  note  that  the  input  impedance  of 
the  virtual-ground  amplifier  is  of  the  order  of  1011  ohms  so  that  the 
filter  components  are  not  at  all  loaded.   The  feedback  filter  is 
convenient  for  minimizing  the  system  settling  time  following  a  change- 
over from  one  electrode-tissue  system  to  another.   Hence,  during  such 
changes,  one  uses  the  smallest  resistor  value  in  the  filter. 

The  Current  Control  Circuit 

For  many  applications,  the  resistance  of  a  simple  current  source 
(Figure  3.6)  would  have  to  be  impractically  large,  resulting  in  poor 
frequency  response  and  excessive  driving  voltage  requirements.   The 
scheme  which  was  chosen  for  the  experimental  system  is  capable  of 
reducing  these  problems  by  several  orders  of  magnitude.   The  basic 
configuration  may  be  found  in  most  applications  manuals  for  analog 
computing  amplifiers  (e.g.,  Philbrick  [1966];  Burr-Brown  [1963]). 

It  is  instructive  to  present  the  current  control  circuit  in  terms 
of  the  elementary  principles  involved.   Consider  the  circuit  shown  in 
Figure  5.4.   Let  the  q  elements  represent  a  simple  current  source  and 
let  the  f  elements  designate  a  correction  circuit  applied  as  shown  to 
the  electrode-tissue  system.   It  is  desired  to  control  If  such  that 

E 
Iq  =  jf  (5.2) 

q 
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Figure  5.4.   Principles  of  Current  Control 


which  is  independent  of  the  impedance  occurring  at  Ql.   Since 


i  .'i.Sa 

1   R    R 

q   q 


(5.3) 


it  is  evident  that  the  required  value  of  I,  is  V„n/R  .   Noting  that 

r     Ql   q        • 


E    V 
I  -  £    Qi 

f   Rf   Rf 


C5.H) 


the  specification  for  the  f  generator  may  be  written 


Ef        Rf 

Vm       R 
Ql        q 


(5.5) 
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The  right  side  of  Equation  (5.5)  is  equal  to  the  ideal  transfer  function 
of  a  common  voltage  follower  (Figure  A2)  if 


Rl    Rq 


(5.6) 


Hence  the  complete  current  control  circuit  may  be  drawn  as  in  Figure 
5.5. 
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Figure  5.5.   The  Current  Control  Circuit 


The  exact  output  of  the  current  control  circuit  may  be  written 


q 


(5.7) 
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The  term,  5q,  is  the  current  control  error  and  is  analogous  to  a  and 
of  Chapters  3  and  4.   It  is  given  by 


VRq 


-  1 


(5.8) 


where  Kf  is  the  exact  transfer  function  of  the  voltage  follower  given 

by  Equation  A2-1.   g„  is  the  current  control  error  that  would  result  if 

E_  and  R  were  used  without  the  correction  circuit  and  is  defined  bv 
q      q  -> 

Equation  (4.31).   When  the  amplifier  gain  and  common  mode  rejection 
ratio  are  sufficiently  high  (see  Equations  [A2-4]  and  [A2-5]),  Equation 
(5.8)  reduces  to 


R2/R1 
V*q~ 


(5.9) 


which  goes  to  zero  if  Equation  (5.6)  is  exactly  satisfied. 

As  a  practical  illustration,  assume  that  the  difference  term  in 
Equation  (5.9)  is  adjusted  to  0.01.   Then  if,  as  in  previous  examples, 
the  magnitude  of  the  current  control  error  (here  | 5  | )  is  limited  to 
less  than  0.01,  the  generator  resistance  requirement  becomes 


R  >   Z   th.,  th„  t(h„ 
q   . '  Ql    31    11    v  ll1 


(5.10) 


where  use  has  been  made  of  Equation  (4.33).   This  condition  is  signifi- 
cantly less  stringent  than  that  for  the  simple  current  source  in 
Equation  (4.34).   In  many  electrode-tissue  systems,  8n  is  significant 
only  at  very  low  frequencies  where  the  interface  impedance  becomes  large. 
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In  such  cases,  the  difference  terra  in  Equation  (5.8)  becomes  unimportant 
at  higher  frequencies,  and  a  low  frequency  amplifier  may  suffice  in 
the  current  control  circuit. 

In  the  experimental  system,  R  and  R_  are  10  kilohms,  and  R  and 
Rf  are  1  megohm  each.   It  is  desirable  not  to  make  the  ratio  of  R2  "to  R, 
large  to  avoid  reducing  the  dynamic  range  for  current  control.   Since 
the  amplifier  used  has  an  output  voltage  range  of  approximately  ±10 
volts,  the  present  circuit  is  capable  of  controlling  up  to  ±5  micro- 
amperes into  a  1  megohm  load  at  Ql.   The  current  drawn  by  the  input 
admittances  of  the  amplifier  must  have  negligible  effect  on  circuit 
performance.   This  condition  is  virtually  assured  here  since  the 
impedance  at  either  input  is  determined  by  a  field  effect  transistor 
and  corresponds  to  10   ohms  in  parallel  with  5  picofarads . 

In  the  actual  system,  R  is  on  the  circuit  board  and  is 
connected  to  the  voltage  generator  through  several  feet  of  coaxial 
cable.   The  shield  of  this  cable  is  grounded  at  the  generator.   The 
current  output  of  the  control  circuit  is  fed  to  the  electrode-tissue 
system  connector  by  one  of  the  four  miniature  coaxial  cables  discussed 
earlier.   The  shield  of  this  cable  is  tied  to  a  voltage  divider  at  the 
output  of  the  current  control  amplifier  in  similar  fashion  as  was  done 
for  the  shield  of  the  PI  cable  in  Figure  5.3. 

Experimental  Studies 

The  experimental  virtual-ground  system  provides  current  control 
for  the  applied  signal  and  amplification  for  the  voltage  response. 
Additional  instruments  are  necessary  to  process  and  display  the  signals 


162 


for  impedance  determinations.   In  the  experiments  to  be  described,  the 
input  (E  in  Figure  5.4)  and  output  (Figure  5.3)  voltages  of  the  systen 
were  continuously  monitored  on  a  dual-trace  oscilloscope  (see  Figure 
5.2).   This  display  was  convenient  for  observing  approximate  amplitude 
and  phase  relationships  between  the  two  signals  and  for  ascertaining 
that  the  system  was  functioning  properly. 

Precise  measurements  were  made  with  a  model  JB-6  Lock-in 
Amplifier  from  Princeton  Applied  Research  Inc.  (see  Appendix  4  for 
principles  of  operation).   When  suitably  calibrated,  this  instrument 
continuously  displays,  on  separate  meters,  the  relative  in-phase  and 
quadrature  components  of  the  system  output  with  respect  to  the  system 
input,  or  any  other  reference  signal.   The  system  input  signal  was 
derived  from  the  internal  generator  of  the  JB-6.   Aside  from  the 
convenience  of  a  direct  meter  display,  the  lock-in  amplifier  permits 
measurements  to  be  made  in  the  presence  of  signal-to-noise  ratios  much 
below  unity  (Moore,  1962). 

The  experimental  results  are  reported  in  two  separate  sections 
designated:   Performance  Tests  and  Four-Electrode  Impedance  Measure- 
ments.  The  measurement  techniques  used  are  incorporated  in  the  block 
diagram  of  Figure  5.6.   The  components  enclosed  in  the  dashed  block 
constitute  the  experimental  virtual-ground  system.   Switch  S2  symboli- 
cally represents  the  double  female  connector  and  the  male  connectors 
described  previously.   Additional  information  about  measurement 
techniques  is  given  as  needed  in  the  discussions  that  follow. 

Performance  Tests 
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The  success  of  an  impedance  experiment  depends  on  the  overall 
performance  of  the  combined  electrode-tissue  and  instrumentation  systems. 
The  three  experiments  to  be  described  are  aimed  at  outlining  the  basic 
capabilities  of  the  present  virtual-ground  system.   Simple  circuit 
arrangements  are  used  to  permit  an  independent  evaluation  of  the 
frequency  response,  current  control  and  common  mode  rejection 
characteristics  of  the  system. 

Frequency  response  test 

The  frequency  response  of  the  virtual-ground  system  is  primarily 
of  interest  in  cases  where  nonsinusoidal  input  waveforms  are  to  be 
used.   Hence,  if  a  rectangular  wave  were  applied  to  the  current-control 
circuit  and  the  electrode-tissue  transfer  impedance  (m,  )  wera  a  pure 
resistance,  the  output  waveform  of  the  virtual-ground  system  should 
also  be  a  rectangular  wave,  differing  only  in  amplitude  from  the  input 
signal.   If  the  electrode-tissue  transfer  impedance  includes  reactive 
elements ,  then  the  output  waveform  should  accurately  represent  the 
effects  of  these  elements  so  that  a  spectral  analysis  can  be  made. 

To  measure  the  frequency  response,  a  1/4  watt  carbon  resistor 
of  10  or  100  ohms  was  used  as  the  test  subject  in  place  of  the  electrode- 
tissue  system.   Since  the  resistor  is  a  two  terminal  device,  terminals 
Ql  and  PI,  and  P2  and  Q2,  were  coincident.   This  arrangement  eliminated 
the  possibility  of  current  control  and  common  mode  rejection  errors, 
which  could  have  masked  the  true  frequency  response.   The  reference 
voltage  for  this  test  was  derived  from  a  simple  resistive  attenuator 
as  shown  in  Figure  5.7.   This  circuit  provides  an  adjustable  frequency- 
independent  fraction  of  the  JB-6  output  voltage. 
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Figure  5.7.   Attenuator  Circuit  for  Frequency  Response  Test 

The  first  step  in  the  procedure  was  to  adjust  the  attenuation 
ratio  of  the  attenuator  to  equal  the  overall  transfer  function  of  the 
virtual-ground  system  and  test  resistor.   This  was  done  at  an  inter- 
mediate frequency  (e.g.,  100  Hz)  where  the  transfer  function  was  known 
to  be  real  and  constant.   Then,  the  attenuator  was  maintained  at  this 
setting  to  provide  the  reference  voltage  at  all  other  test  frequencies. 
The  procedure  at  each  frequency  was  to  calibrate  the  JB-6  using  the 
attenuator  output  and  then'  to  switch  over  to  the  virtual-ground  system. 
Magnitude  and  phase  errors  were  detected  as  changes  in  the  JB-6  meter 
readings  from  the  calibration  values. 

At  frequencies  of  2  Hz ,  10  Hz,  100  Hz,  1  kHz,  and  10  kHz,  no 
phase  or  magnitude  errors  could  be  detected.   The  results  at  higher 
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frequencies  are  indicated  in  Table  5.1.   As  a  practical  estimate,  the 
maximum  uncertainty  in  the  tabulated  results  is  approximately  ±1  percent 
for  the  magnitude  errors  and  ±1  degree  for  the  phase  errors.   At 
frequencies  of  10  kHz  and  below,  the  uncertainties  in  the  zero  errors 
observed  are  smaller. 

TABLE  5.1 

FREQUENCY  RESPONSE  ERRORS 


Frequency 

Magnitude  Error 

Phase  Error 

Hz 

percent 

degrees 

40  k 

+1.4 

0 

70  k 

0 

• 

-1.7 

85  k 

+1.3 

0 

100  k 

+1.6 

-1.8 

150  k 

+  8.6 

-3.3 

The  errors  shown  in  Table  5.1  are  probably  within  tolerable  limits 
for  most  nonsinusoidal  measurements  in  tissue  at  the  present  state  of 
the  art.   However,  it  is  reasonably  certain  that  the  frequency  response 
could  be  improved  using  simple  compensation  techniques,  should  the 
occasion  require  it.   One  such  technique  was  mentioned  previously  in 
the  description  of  the  response  amplifier  circuit. 

Common  mode  rejection  test 

Common  mode  rejection  capability  is  of  interest  because  it 
establishes  certain  limits  for  the  type  of  electrode  configuration  which 
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can  be  accurately  measured  with  a  given  virtual-ground  system.   More 
specifically,  it  determines  how  large  the  impedance  ratio  of  Equation 
(4.6)  may  become  before  |3V|  exceeds  an  arbitrary  maximum  value.   It  is 
evident  that  common  mode  rejection  capability  is  directly  proportional 
to  the  gain  of  the  virtual-ground  amplifier. 

To  test  the  common  mode  rejection  capability  of  the  experimental 
system,  the  circuit  shown  in  Figure  5.8  was  used  as  the  test  subject 
in  place  of  the  electrode-tissue  system.   The  calibrate  resistance  and 
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Figure  5.8.   Circuit  for  Common  Mode  Rejection  Test 

the  attenuator  of  Figure  5.6  were  not  used  during  this  test.   At  each 
frequency  of  measurement ,  the  procedure  was  to  set  Z  equal  to  zero  and 
calibrate  the  JB-5  against  R  .   Then  Z  was  increased  in  arbitrary  steps 
until  the  magnitude  or  phase  error  due  to  excessive  common  mode  signal 
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became  large.   The  errors  were  detected  as  changes  in  the  JB-6  meter 
readings  from  their  calibration  values.   At  very  low  frequencies,  Z 
was  derived  from  a  decade  capacitance  box.   At  other  frequencies,  it 
was  derived  from  a  decade  resistance  box.   This  technique  permitted 
phase  shift  to  be  used  as  the  most  sensitive  indicator  of  error  for 
all  measurements. 

The  results  of  the  common  mode  rejection  test  are  summarized 
in  Table  5.2.   Notice  that  with  each  decade  increase  in  frequency  above 
100  Hz,  the  value  of  |z|/Rp,  at  which  a  similar  phase  error  occurs, 
decreases  by  roughly  one  decade.   At  frequencies  of  100  Hz  and  lower, 
the  ratio  corresponding  to  similar  phase  errors  is  much  less  dependent 
on  frequency.   Furthermore,  it  is  evident  that  the  phase  angle  was  the 
predominant  indicator  of  error  at  all  frequencies ,  and  that  this  was 
accomplished  by  making  Z  reactive  at  frequencies  below  100  Hz  and 
resistive  at  100  Hz  and  above.   The  foregoing  observations  are  in 
agreement  with  the  nominal  gain  and  phase  shift  characteristics  given 
for  the  virtual-ground  type  of  amplifier  in  Figure  4.2. 

It  is  apparent  that,  at  100  kHz,  the  present  system  does  not 
possess  the  capability  expressed  by  Equation  (4.8)  if  |g  |  is  required 
to  be  less  than  0.01.   The  equation  holds  in  the  present  case,  however, 
if  the  limit  for  | 3  |  is  taken  as  0.04.   In  other  words,  the  gain  of 
the  virtual-ground  amplifier  in  the  experimental  system  is  approximately 
250  at  100  kHz  rather  than  the  conjectured  practical  upper  limit  of  1000. 

The  frequency  dependence  of  the  common  mode  rejection  capability 
of  the  virtual-ground  system  is  especially  appropriate  for  measurements 
in  electrolytic  (biological)  media.   Consider  for  example,  the  impedance 
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TABLE  5.2 


COMMON  MODE  REJECTION  ERRORS* 


Frequency      | Z | /Rp      Magnitude  Error      Phase  Error 
Hz  percent degrees 


2 

2,700 

8,000 

16,000 

0 
<0.2 
<2 

0 
-2.0 

-4.7 

10 

1,600 
16,000 
32,000 

0 
<0.5 
<1 

0 
-3.1 
-7.4 

100 

1,600 
10,000 

0 
2.7 

0 
4.2 

Ik 

100 
1,000 

0 
<0.5 

0 

3.1 

10k 

10 
100 

0 
<0.5 

0 
2.4 

70k 

1 
10 

0 
<0.5 

0 
1.4 

100k 

1 

5 
10 
20 

0 

<0.1 
<0.5 
<0.5 

0 

0.9 
2.8 
4.8 

>Z  is  capacitive  reactance  at  2  Hz  and  10  Hz,  and  resistance 
at  all  other  frequencies .   Results  were  obtained  with  Rp 
equal  to  100  ohms.   Similar  results  were  obtained  with  Rp 
equal  to  10  ohms . 
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ratio  of  Equation  (4.7).   At  high  frequencies,  the  interface  impedance, 
Zq2»  becomes  small  and  the  ratio  approaches  a  constant  value  determined 
mainly  by  the  geometry  of  the  electrode  system.   As  the  frequency  is 
decreased,  the  ratio  increases  and  Z   eventually  becomes  the  dominant 
factor  in  the  numerator.   Since  Zq2  varies  according  to  a  fractional 
power  of  frequency  (as  noted  in  Chapter  1),  it  increases  more  slowly 
than  K  ,  which  depends  on  the  first  power  of  frequency  (see  Figure  4.2). 
Therefore,  if  Equation  (4.7)  is  satisfied  at  a  given  test  frequency,  it 
will  almost  certainly  be  satisfied  (to  an  even  greater  extent)  at 
frequencies  below  the  test  frequency.   More  precisely,  since  |k.| 
becomes  constant  at  very  low  frequencies ,  the  inequality  will  eventually 
become  invalid  at  some  frequency  below  the  break  in  |k  |. 

When  measuring  an  unknown  electrode-tissue  system,  the  foregoing 
considerations  imply  that  it  is  sufficient  to  ascertain  that  Vp9  is 
negligible  compared  to  Vp  (i.e.,  1 8V  |  falls  below  an  arbitrary  limit) 
at  the  highest  and  lowest  frequency  in  the  measurement  band.   Then  the 
common  mode  rejection  error  will  also  be  negligible  at  all  other 
frequencies  in  the  band.   In  most  cases,  only  the  high  frequency  check 
is  necessary. 

Current  control  test 

Current  control  capability  determines  how  large  the  electrode- 
tissue  impedance  between  the  Ql  terminal  and  ground  (given  by  the  right 
side  of  Equation  (5.10)  may  become  before  the  error  in  the  applied 
current  exceeds  an  arbitrary  maximum  value.   It  is  primarily  important 
at  low  frequencies  where  the  interface  impedance  of  the  Ql  electrode 
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may  become  very  large.   Current  control  capability,  per  se,  is  fairly 
independent  of  frequency  and  is  rarely  a  limiting  factor  at  higher 
frequencies  due  to  the  predominance  of  loading  and  common  mode  rejection 
errors . 

To  test  the  current  control  capability  of  the  experimental  system, 
the  circuit  shown  in  Figure  5.9  was  used  as  the  test  subject  in  place 
of  the  electrode-tissue  system.   The  calibrate  resistance  and  the 


Ql  ♦ 


02 


♦  P2 


Figure  5.9.   Circuit  for  Current  Control  Test 


attenuator  of  Figure  5.6  were  not  used  during  this  test.   At  each 
frequency  of  measurement,  the  procedure  was  to  set  Z  equal  to  zero  and 
calibrate  the  JB-6  against  Rp.   Then  Z  was  increased  in  arbitrary  steps 
until  the  magnitude  or  phase  error  became  large  as  evidenced  by  changes 
in  the  JB-6  meter  readings. 
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The  results  of  the  current  control  test  are  summarized  in  Table 
5.3.   Since  R  ,  in  the  present  system,  is  1  megohm,  the  results  reported 


TABLE  5.3 


CURRENT  CONTROL  TEST  RESULTS- 


Frequency 

1*1 

Magnitude 

Error 

Phase  Error 

Hz 

ohms 

percent 

degrees 

10 

16k 

0 

0 

160k 

0 

0 

1.6M 

-1.8 

2.1 

3 .  2M 

-3.3 

4.3 

100k 

500 

0 

0 

2.7k 

<0.1 

-1.9 

10k 

1.6 

-3.8 

*Z  is  capacitive  reactance  at  10  Hz  and  resistance  at  100 
kHz.   Results  were  obtained  with  Rp  equal  to  100  ohms. 

for  10  Hz  are  in  approximate  agreement  with  the  requirement  expressed 
by  Equation  (5.10)  for  a  current  control  error  of  less  than  1  percent. 
The  errors  obtained  at  100  kHz  are  not  current  control  errors.   They 
represent  combined  loading  errors  and  have  been  included  only  for 
comparison. 

Noise  levels 


Noise  accompanying  the  signal  output  of  the  experimental  virtual- 
ground  system  has  been  observed  on  the  monitoring  oscilloscope.   Using 
a  scope  amplifier  with  a  3  db  bandwidth  of  500  kHz,  the  broadband 
peak-to-peak  instrumentation  noise  has  been  between  20  and  40  microvolts 
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referred  to  the  PI  terminal.   The  peak-to-peak  noise  due  to  sixty- 
cycle  pick-up  has  consistently  fallen  within  the  range  of  the  instru- 
mentation noise.   This  low  sixty-cycle  level  is  obtained  without  the 
use  of  shielded  enclosures,  and  in  the  presence  of  fluorescent  lighting 
under  ordinary  laboratory  conditions. 

Four-Electrode  Impedance  Measurements 

The  measurements  reported  in  this  section  were  all  made  using 
a  substitution  technique  suggested  by  Dr.  Jack  R.  Smith.   With  reference 
to  Figure  5.6,  switch  S,  was  maintained  on  the  system  output.   The 
initial  step  at  each  frequency  of  interest  was  to  calibrate  the  JB-6 
for  phase  measurements  with  S„  set  on  the  calibrate  resistor  (the 
decade  resistance  box  in  Figure  5.2).   This  resulted  in  a  zero  reading 
on  the  quadrature  meter  and  an  arbitrary  reading  on  the  in-phase  meter 
which  was  proportional  to  the  resistance  used.   The  resistance  was 
generally  chosen  to  approximate  the  anticipated  value  of  the  transfer 
impedance  to  be  measured. 

The  next  step  was  to  set  S2  on  the  electrode-tissue  (or  other 
medium)  system  and  record  the  quadrature  and  in-phase  readings  of  the 
JB-6.   The  ratio  of  the  quadrature  to  in-phase  reading  gave  the  tangent 
of  the  transfer  impedance  phase  angle. 

The  final  step  was  to  switch  S„  back  to  the  calibrate  resistor 
and  adjust  the  resistance  until  the  in-phase  reading  was  equal  to  that 
observed  for  the  electrode-tissue  system.   This  resistance  was  equal  to 
the  resistive  part  of  the  transfer  impedance.   Hence,  knowing  the  phase 
angle  and  resistive  part,  the  reactive  part  and  the  magnitude  of  the 
transfer  impedance  could  be  determined. 
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The  substitution  technique  has  the  advantage  that  amplifier  gain 
and  frequency  response  characteristics  need  not  be  known  precisely. 
In  addition,  it  minimizes  the  chances  for  calibration  errors  since 
these  will  tend  to  be  similar  for  both  the  calibration  resistor  and 
the  unknown  electrode-tissue  system.   The  substitution  technique, 
however,  does  not  eliminate  the  effects  of  common  mode  rejection  and 
current  control  error  since  the  calibrate  resistor  is  a  two-terminal 
device  and  the  electrode-tissue  system  is  a  four-terminal  device. 
Therefore,  these  properties  should  be  checked  during  experiment  as 
previously  described. 

Electrode  development 

An  important  part  of  the  current  program  has  been  to  develop 
more  reliable  electrode  systems,  especially,  for  use  in  situ  at  very 
low  frequencies.   The  problems  stem  mainly  from  transverse  admittance 
in  the  electrode  support  materials  (Chapter  3)  and  boundary  condition 
effects  at  the  electrode-electrolyte  interfaces  (Chapter  2),  which 
occur  independently  of  the  instrumentation  technique  used.   The  virtual- 
ground  system  has  made  it  possible  to  evaluate  the  performance  of 
various  electrode  systems  (in  known  test  media)  independently  of  errors 
which  were  previously  introduced  by  less  reliable  instrumentation 
techniques. 

One  of  the  more  recent  electrode  assemblies  developed  during  the 
current  program  is  illustrated  in  Figure  5.10.   The  square  arrangement 
of  the  electrodes  was  designed  by  Dr.  Arnold  H.  Nevis  and  is  the 
outgrowth  of  a  suggestion  by  the  author  to  place  the  potential  electrodes 
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epoxy  pedestal 


8  mm 


2.5  mm 


2.5  mm 


5  mm 


1.  5  mm 


bottom  view 


side  view 


Figure  5.10.   Electrode  Assembly  for  In  Situ  Impedance  Measurements 


in  a  region  of  low  potential  gradient  (see  discussion  of  Equation  [2. 14] 
and  Figure  [2.3]).   Stainless  steel  and  platinum  have  been  used  for  the 
electrode  materials.   In  the  latter  case,  the  surfaces  were  platinized. 
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Experiments  are  in  progress  to  evaluate  the  effects  of  insulating 
various  portions  of  the  electrode  surfaces. 

Each  electrode  is  soldered  to  a  flexible  piece  of  insulated 
wire  approx5_mately  one  inch  in  length.   The  electrodes  and  wires  are 
cast,  as  indicated,  in  a  cylindrical  epoxy  pedestal.   Connections  to 
the  wire  ends  are  made  by  means  of  small  covered  alligator  clips  which 
are  attached  to  the  leads  of  the  male  connector  of  the  virtual-ground 
system  (see  rat's  head  in  Figures  5.1  and  5.2).    By  bringing  the 
exposed  connections  away  from  the  electrode  assembly,  problems  due  to 
transverse  admittance  at  very  low  frequencies  have  virtually  been 
eliminated.   Furthermore,  the  alligator  clips  are  especially  easy  to 
apply  to  unrestrained  animals. 

A  typical  set  of  test  results  for  an  electrode  assembly  of  the 
type  indicated  in  Figure  5.10  using  0.18  mm  diameter  platinized  platinum 
wires  are  given  in  Table  5.4.   The  transfer  impedance  (itvq)  was  found 
to  be  purely  resistive  and  independent  of  frequency,  except  at  100  Hz 
where  an  anomalous  2.1  percent  increase  was  observed.   The  source  of 
this  increase  has  not  been  determined.   Additional  and  more  carefully 
controlled  electrode  studies  will  be  necessary  before  valid  explanations 
can  be  given  for  the  anomalous  results  which  are  sometimes  observed. 
Sources  of  error  may  lie  in  human  factors ,  in  malfunctions  of  the  JB-6 
and  in  boundary  condition  effects  in  the  electrode-saline  system. 

The  fifth  column  in  Table  5.4  indicates  the  ratio  of  common  mode 
to  differential  signal  of  the  electrode-saline  system.   The  sixth 


The  epoxy  pedestal  and  wire  connection  technique  were  developed 
by  William  K.  Converse. 
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TABLE  5.4 


TEST  RESULTS  Or  AN  ELECTRODE  ASSEMBLY  IN  ISOTONIC  SALINE 


Freq 

mll 

Magnitude 

Phase  Error 

Hz 

ohms 

Error 

degrees 

percent 

2 

14.2 

0 

0 

10 

14.2 

0 

0 

100 

It.  5 

+  2.1 

0 

Ik 

14.2 

0 

0 

50k 

14.2 

0 

0 

100k 

14.2 

0 

0 

Vr 


1*1 

ohms 


230 
60 


5.9 


3,200 

850 


84 


column  indicates  the  effective  impedance  at  the  Ql  terminal  which  the 
current  source  must  drive.   It  is  apparent  from  the  performance  data 
in  Tables  5.2  and  5.3  that  the  parameters  of  the  electrode-saline  system 
are  well  within  the  capabilities  of  the  experimental  virtual-ground 
system  except  at  100  kHz.   At  this  frequency,  the  ratio  of  common  mode 
to  differential  signal  is  in  the  range  (greater  than  5)  where  the  common 
mode  rejection  errors  start  to  become  significant.   This  condition 
could  be  improved  with  more  efficient  placement  of  the  potential 
electrodes  (as  indicated  in  Figure  2.3). 

If  the  electrode  assembly  of  Table  5.4  were  used  in  tissue,  the 
values  in  the  |z|  column  would  be  greater  than,  while  the  values  in  the 
[ V„  1  / | Vp  j  column  would  be  similar  to  those  in  isotonic  saline.   Since 
the  |z|  column  is  well  within  system  capabilities,  the  measurement  of 
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tissue  impedance  would  not  be  any  more  challenging  to  the  virtual- 
ground  system  than  the  measurement  of  isotonic  saline  impedance. 

Tissue  impedance  measurements 

The  experimental  virtual-ground  system  has  been  used  successfully 
to  measure  impedance  in  the  cerebral  cortex  of  rats  and  in  the  femur 
bone  of  rabbits.   The  cortical  studies  are  part  of  a  research  program 
conducted  by  Doctors  Arnold  H.  Nevis,  Jack  R.  Smith,  and  George  H. 
Collins.   The  bone  measurements  were  part  of  a  research  project 
conducted  by  medical  student  Jim  Mandel.   These  applications  are 
especially  noteworthy  because  they  represent  two  widely  different 
impedance  levels,  the  bone  impedance  being  roughly  five  times  that  of 
the  cortex.   Several  of  the  cortical  experiments  are  discussed  in  the 
following  paragraphs  to  demonstrate  the  quality  of  results  already 
obtained  with  the  virtual-ground  technique. 

Electrode  assemblies  of  the  type  shown  in  Figure  5.10  have  been 
chronically  implanted  in  the  cerebral  cortex  of  rats.   The  measurement 
techniques  were  described  earlier  in  this  section.   The  experimental 
arrangement  is  shown  in  Figures  5.1  and  5.2.   The  current  density  used 

at  the  current  electrodes  (area  approximately  l.U  mm  )  was  usually  less 

2 

than  1  microampere  peak-to-peak/mm  ,  and  the  response  voltages  were 

usually  less  than  100  microvolts  peak-to-peak.   In  the  experiments  to 
be  described,  the  electrodes  were  made  of  platinum  wire  (0.18  mm 
diameter)  and  were  platinized. 

The  results  plotted  in  Figure  5.11  characterize  seven  separate 
frequency  runs  measured  in  the  same  rat  on  seven  different  days  during 
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a  15-day  interval.   Only  one  of  the  frequency  runs  included  a  measure- 
ment at  2,  50,  500,  5k,  and  50k  Hz  (shown  by  the  small  circles).   The 
vertical  lines  at  10,  100,  Ik,  10k,  and  100k  Hz  indicate  the  spread  of 
the  majority  of  the  seven  experimental  values  obtained  at  each  frequency. 
Individual  values  not  within  the  lines  are  indicated  by  a  symbol 
designating  a  particular  frequency  run  (i.e.,  the  cross,  circle, 
triangle,  and  square). 

The  frequency  run  indicated  by  the  crosses  is  a  special  case 
which  was  obtained  after  the  other  six  runs  and  is  discussed  in  a 
later  paragraph.   The  six  remaining  spectra  were  measured  within  a 
thirteen-day  interval.   In  the  magnitude  plot,  the  magnitude  spread  is 
4  ohms  or  less ,  except  at  10  Hz  where  it  is  5  ohms .   The  largest 
percentage  spread  occurs  at  100  kHz  and  is  approximately  7  percent  of 
the  midspread  impedance  (which  is  about  57  ohms).   In  the  original 
data,  the  individual  magnitude  changes  at  each  frequency  tended  to  be 
in  the  same  direction  so  that  the  shape  of  the  spectrum  did  not  vary 
significantly.   The  total  magnitude  drop  from  10  Hz  to  100  kHz  is 
approximately  20  ohms  or  25  percent  of  the  value  at  10  Hz.   The 
magnitude  decrease  at  2  Hz  is  contrary  to  previous  experience  and  may 
be  due  to  the  increased  uncertainty  and  difficulty  of  operating  the 
JB-6  at  very  low  frequencies  in  the  presence  of  the  electrocortigram. 

In  the  phase  plot, "the  phase  spread  is  less  than  1.5  degrees 
and  all  angles  are  zero  or  negative  (capacitive) ,  except  at  100  kHz. 
The  largest  measured  phase  angle  occurs  at  100  Hz  and  lies  between  -3.6 
and  -4.8  degrees.   The  phase  angles  corresponding  to  the  magnitude  plot 
marked  by  the  crosses  consistently  fell  in  the  majority  group  for  all 
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runs.   At  10  Hz,  five  of  the  seven  runs  gave  zero  phase  angle  (and 
zero  spread).   The  two  anomalous  points  could  conceivably  be  due  to 
greater  variability  of  the  tissue  impedance  at  10  Hz.   However,  it  is 
more  probable  that  these  points  are  due  to  errors  in  the  procedure  of 
operating  the  J3-6. 

The  positive  phase  angles  at  100  kHz  may  be  explained  by  the 
limited  common  mode  rejection  capability  of  the  present  virtual-ground 
system  (see  Table  5.2).   The  electrode  system  used  had  a  common  mode 
to  differential  signal  ratio  greater  than  7  at  100  kHz.   An  experimental 
check  of  the  voltage  at  P2  indicated  a  positive  reactive  component 
which  could  account  for  the  anomalous  phase  angles.   Errors  of  this 
type  are  influenced  by  the  gain  characteristics  of  the  virtual-ground 
amplifier  and  certain  other  parameters  which  may  vary  from  day  to  day. 

It  is  apparent  that  more  experimental  points  are  necessary  to 
fully  evaluate  the  frequency  spectrum,  especially  at  the  very  low 
frequencies.   Unfortunately,  the  very  low  frequency  measurements  are 
the  most  difficult  to  make  due  to  the  presence  of  strong  electrocorti- 
graphic  activity.   They  also  involve  lengthy  calibration  procedures 
because  of  the  extended  averaging  time  and  large  filter  time  constants 
which  the  JB-6  requires  at  these  frequencies. 

On  each  of  the  last  nine  days  in  the  15-day  observation  interval, 
the  rat  received  an  intraperitoneal  injection  of  triethyl  tin  sulfate. 
This  compound  is  known  to  be  toxic  and  after  a  sufficient  number  of 
injections  causes  brain  edema.   The  clinical  symptoms  include  paralysis 
which  begins  in  the  hind  limbs  and  continues  to  spread  forward  to  the 
forelimbs  and  respiratory  muscles  until  the  animal  dies. 
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The  frequency  run  designated  by  the  crosses  was  obtained  two 
days  after  the  last  of  the  other  six  spectra  were  measured.   The 
impedance  rise  was  approximately  correlated  with  the  appearance  of  hind 
limb  paralysis.   As  in  the  case  of  the  day  to  day  variations,  the 
magnitude  increase  due  to  the  injections  was  fairly  uniform  at  all 
frequencies  (approximately  11  ohms)  and  the  shape  of  the  spectrum  was 
preserved.   The  phase  angle  values,  as  previously  noted,  were  not 
affected  by  the  injections. 

It  is  interesting  to  note  the  magnitude  of  the  specific  impedance 
of  the  cortical  tissue  represented  by  the  results  in  Figure  5.11.   When 
tested  in  isotonic  saline  (approximately  60  ohm-cm.  specific  impedance), 
the  electrode  system  indicated  a  transfer  impedance  of  12.5  ohms 
independent  of  frequency.   The  high  and  low  frequency  calculations  for 
the  specific  impedance  of  the  cortex  are  indicated  in  Table  5.5.   The 

TABLE  5.5 

SPECIFIC  IMPEDANCE  OF  RAT  CEREBRAL  CORTEX 


Frequency 
Hz 


Normal 


Paralyzed 


10 


100k 


370 


280 


160 


330 


ohm- cm 


normal  values  are  not  too  different  from  those  found  by  Ranck  (1963a) 
in  acute  experiments  on  rabbit  cortex. 
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The  present  data  are  insufficient  to  determine  the  exact  location 
and  the  value  of  the  maximum  phase  shift.   However,  it  appears  to  be 
no  greater  than  5  degrees  and  to  lie  in  the  50  to  100  Hz  frequency 
range.   Ranck  (1963a)  reported  a  maximum  phase  shift  of  7  degrees  which 
also  occurred  in  the  50  to  100  Hz  range  in  acute  experiments  on  rabbit 
cortex.   He  also  found  that  the  phase  angle  went  to  zero  at  5  kHz  and 
increased  again  at  frequencies  above  20  kHz  (a  possible  3  dispersion). 
In  the  present  case,  the  phase  angle  does  not  approach  zero  until 
frequencies  beyond  50  kHz.   Except  for  this  order  of  magnitude  difference 
in  frequency  of  closure ,  the  results  found  here  in  rat  cortex  in  chronic 
experiments  are  in  general  agreement  with  the  results  found  by  Ranck  in 
rabbit  cortex  in  acute  experiments. 

In  the  complex  impedance  plane,  the  results  of  Figure  5.11  have 
the  general  form  shown  in  Figure  5.12.   The  plot  suggests  that  unlike 
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Figure  5.12.   Polar  Plot  of  the  Impedance  of  Rat  Cortex 
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that  of  muscle  (Fatt,  1963),  the  impedance  of  cerebral  cortex  does 
not  display  distinctive  circular  arc  behavior  at  frequencies  below  100 
kHz.   Furthermore,  it  will  be  necessary  to  obtain  measurements  above 
100  kHz  in  order  to  determine  whether  the  6  dispersion  is  distinct  from 
the  a  dispersion. 

The  consistency  and  reproducibility  of  the  results  presented 
in  Figure  5.11  demonstrate  a  clear  advance  in  the  field  of  in  situ 
impedance  measurements.   The  day  to  day  magnitude  spread  (excluding 
the  rise  due  to  the  injections)  is  well  within  bounds  of  normal 
physiological  variations.   For  example,  the  author  has  monitored  the 
impedance  (at  1  kHz)  of  several  normal  rats  who  accidently  were  not 
fed  for  one  or  two  days.   In  a  typical  case,  the  impedance  dropped  from 
94  to  85  ohms  within  the  ten-minute  interval  after  food  was  put  in  the 
cage.   These  results  were  correlated  with  visible  symptoms  of  anxiety 
and  oscillatory  bursts  in  the  electrocortigram. 

The  similarity  of  the  seven  impedance  spectra  implies  that 
significant  physiological  variations  did  not  occur  during  the  frequency 
runs,  even  though  it  usually  took  from  one-half  to  one  hour  to  complete 
a  run.   Unfortunately,  time  invariance  to  this  degree  may  not  be 
possible  in  many  important  experiments  to  come.   Hence,  if  spectral 
measurements  are  to  be  meaningful  in  such  cases,  it  will  be  necessary 
to  develop  data-handling  techniques  with  sufficient  accuracy  and  speed 
to  obtain  "instantaneous"  spectral  information  either  on-line  or  from 
stored  data.   The  virtual-ground  system  is  one  major  link  in  the  required 
instrumentation  chain. 


RECOMMENDATIONS  FOR  FUTURE  STUDY 

The  foregoing  work  includes  a  body  of  general  analytical  results 
which  express  various  aspects  of  the  problem  of  four-electrode  impedance 
measurements.   It  is  hoped  that  many  of  these  results  will  find  appli- 
cation in  the  design  and  analysis  of  future  electrode  and  instrumenta- 
tion systems. 

Potential  developments  in  the  field  of  biological  impedance 
measurements  appear  to  have  practically  unlimited  scope.   Topics 
suggested  for  immediate  study  are  summarized  in  four  groups  below. 

Electrode  Systems 

a.  Investigation  of  interfacial  boundary  effects  in  various 
electrode  configurations  using  different  electrode  materials. 

b.  Development  of  more  definable  current  fields  using  guard 
electrode  principles. 

c.  Development  of  more  appropriate  media  than  saline  for 
testing  electrodes. 

The  Virtual-Ground  System 

a.  More  precise  compensation  adjustments  to  flatten  and  extend 
frequency  response. 

b.  Connection  for  differential  operation  and/or  the  use  of  a 

wider  band  higher  gain  virtual-ground  amplifier  to  extend  common  mode 

rejection  capability. 
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c.  More  precise  adjustment  of  current  control  circuit  to  increase 
low  frequency  capability. 

d.  More  precise  adjustment  of  +1  gain  return  to  PI  cable  to 
extend  frequency  response. 

e.  Detailed  evaluation  and  minimization  of  inherent  loading 
admittances  of  the  virtual-ground  system  and  associated  paraphernalia. 

f.  Development  of  optimum  filters  for  the  feedback  filter  loop. 

g.  Packaging  of  the  virtual-ground  system  and  feedback  filter. 

Information  Retrieval 

a.  Plotting  of  more  points  in  the  frequency  spectrum,  especially 
at  very  low  and  very  high  frequencies . 

b.  Development  of  accurate  techniques  to  store  and/or  process 
instantaneous  spectral  information. 

c.  Development  of  optimum  signal  waveforms  and  sequences  to 
maximize  information  retrieval  in  a  given  observation  period. 

Biological  Studies 

a.  Application  of  the  principles  presented  in  the  foregoing 
work  to  studies  of  tissues  in  addition  to  brain  and  bone  using  suitably 
designed  electrodes. 

b.  Investigation  of  relationships  between  various  physiological 
variables  and  tissue  impedance  properties. 
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APPENDIX  1 

RECIPROCITY  RELATIONS  OF  THE  IMPEDANCE  PARAMETER  MATRIX 

The  reciprocity  relations  of  Equations  (2.4),  (2.5)  and  (2.6) 
may  be  derived  with  the  aid  of  Figure  Al.   The  currents  at  the  terminals 
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Figure  Al.   Voltage  and  Current  Conventions  for  Reciprocity  Relations 

of  the  electrode-tissue  system  are  related  to  the  loop  currents  as 
follows : 


LQ1 


PI 


LP2 


(0   +  IQ  +  I±) 


(ip    +o    -   h) 


(-1   +0   +  0) 

F 


(Al.l) 
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If  Equation  (Al.l)  is  substituted  into  Equation  (2.3),  the  result  may 
be  expanded  and  arranged  in  the  following  form: 


m12  +  m13 


m22  -  m23 


-m32  +  m33 


m 


11 


mll  "  m12 


m21      m21  ~  m22 


m31      m31  +  m32 


XP 

JQ 

h 



(A1.2) 


This  type  of  matrix,  which  relates  terminal  voltages  to  loop  currents, 
is  a  common  form  in  network  theory  and  is  known  to  be  symmetrical  in  a 
system  obeying  reciprocity. 

Hence,  setting  the  corresponding  elements  on  either  side  of  the 
diagonal  equal  yields 


m22  =  m23  +  mll 


m32  =  m12  +  m33  -  mu 


m21  =  m31  +  m22 


(A1.3) 
(ALU) 
(A1.5) 


These  expressions  correspond  in  order  to  Equations  (2.4),  (2.5) 
and  (2.6),  if  m99  is  eliminated  in  Equation  (A1.5)  by  means  of  Equation 
(A1.3). 
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THE  VOLTAGE  FOLLOWER  AMPLIFIER 


The  configuration  of  a  common  voltage  follower  amplifier  is 
shown  in  Figure  A2.   The  exact  transfer  function  of  this  circuit  may 
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.  Figure  A2.   Voltage  Follower  Amplifier 


be  written 


^   (1  +  1/P)(1  +  Rj/R^ 


*_Bi 


1  + 


1  +  R2/Ri 
K' 


(A2.1) 


where  P  is  the  common  mode  rejection  ratio  and  K'  is  the  loaded  open- 
loop  gain  of  the  operational  amplifier.   If  Z'  designates  the  total 
effective  load  at  the  output,  then 
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where  K  is  the  open-loop  unloaded  gain  and  Z  is  the  output  impedance 
of  the  operational  amplifier. 

Ideally,  P  and  K'  are  infinite  and  the  transfer  function  is  given 
by 


K  =  1  +  R2/R  (A2.3) 


which  is  independent  of  amplifier  characteristics.  To  approach  this 
ideal  with  an  error  of  less  than  1  percent ,  the  f ollowing  conditions 
are  sufficient 

|P|  >  100  (A2.4) 

|K|  >  100(1  +  R  /R  )|l  +  Z0/Z'|  (A2.5) 


APPENDIX  3 


THE  VOLTAGE  INVERTER  AMPLIFIER 

The  configuration  of  a  common  voltage  inverter  amplifier  is 
shown  in  Figure  A3.   The  exact  transfer  function  of  this  circuit  may  be 
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written 


Figure  A3.   Voltage  Inverter  Amplifier 


Eo      "Ro/Ri 

K.  =  r=r-  =  

El       I  +  Rj/Rjl 


1  + 


K' 


(A3.1) 


where  K'  is  as  defined  in  Equation  (A2.2). 
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APPENDIX  k 


PRINCIPLES  OF  THE  LOCK-IN  AMPLIFIER 


The  main  components  of  the  lock-in  amplifier  are  indicated  in 
Figure  A 4.   The  signal  of  the  tuned  oscillator  is  applied  to  the 
variable  phase  shifter  and  via  the  reference  channel  output  to  the 
test  system  of  interest.   The  output  of  the  test  system  is  fed  via 
the  signal  channel  input  to  the  tuned  amplifier. 

The  tuned  amplifier  filters  out  a  significant  part  of  the  noise 
which  may  accompany  the  test  signal.   It  also  amplifies  the  signal 
and  shifts  it  90  degrees  to  provide  an  output  in  quadrature  and  an 
output  in  phase  with  the  signal  voltage. 

The  in-phase  and  quadrature  signals  are  applied  to  the  detection 
systems  where  they  are  multiplied  by  the  output  signal  of  the  phase 
shifter  and  averaged  to  give  the  d.c.  voltages,  Ey  and  Eq,  as  indicated. 
The  detection  systems  comprise  a  phase  sensitive  detector  followed  by 
a  d.c.  amplifier- low  pass  filter  combination,  which  provide  additional 
noise  rejection. 

The  usefulness  of  the  lock-in  amplifier  is  contained  in  its  d.c. 
output  voltages  Ey  and  Eq  which  may  be  conveniently  displayed  on 
separate  meters.   Hence,  the  amplitudes  of  Ej  and  Eq  are  directly 
proportional  to  the  signal  amplitude  while  the  phase  of  the  signal  with 
respect  to  the  shifted  reference  signal  is 
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To  illustrate  how  the  lock-in  amplifier  may  be  used,  let  a 
known  resistive  network  be  connected  between  the  reference  and  signal 
channels.   In  this  case,  9  is  zero  and  if  the  phase  shifter  is  adjusted 
until  E„  is  zero,  6  will  also  be  zero.   If  now,  the  resistive  network 
is  replaced  by  an  unknown  system,  Equation  (A4.1)  would  give  the  phase 
shift  of  the  system  transfer  impedance.   To  find  the  magnitude  of  the 
transfer  impedance,  the  resistive  network  would  be  reconnected  and 
adjusted  to  give  the  same  Ej  as  the  unknown  system.   The  value  of  the 
resistive  network  is  then  the  resistive  part  of  the  transfer  impedance 
which  with  the  phase  determines  the  magnitude  and  reactive  part. 

The  foregoing  procedure  may  also  be  used  when  6  is  not  known 
provided  that  it  occurs  as  a  common  error  factor  in  the  calibrate  system 
and  in  the  unknown  system.   This  assumption  is  the  basis  for  the 
substitution  method  used  in  the  transfer  impedance  measurements  of 
Chapter  5 . 
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